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SUMMARY 

A  suite  of  FORTRAN-IV  computer  programs  is  described  which  may  be  used  to 
assist  in  evaluating  relative  combat  aircraft  performance,  using  energy  manoeuvrability 
theory .  The  programs  are  described  in  detail  using  flowcharts,  and  full  operating  instruc¬ 
tions  are  given.  A  selection  of  outputs  illustrates  the  graphical  and  printed  capabilities  of 
the  suite. 
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NOTE  ON  UNIT  SYSTEMS 


The  SI  system  of  units  is  used  as  the  major  unit  system  in  this  report.  Whenever  appropriate, 
Imperial  equivalents  are  given  in  parentheses  ( ), 

The  programs  described  are  capable  of  processing  data  in  both  the  SI  and  Imperial  systems, 
since  aeronautical  practice  stili  makes  use  of  the  Imperial  system. 
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Definition  Unit  or  Value 

Sonic  speed  . m/s  (ft/s) 

Drag  coefficient . — 

Induced  drag  coefficient . — 

Minimum  drag  coefficient  . — 

Store  drag  coefficient  . .  .  . .  . .  — 

Lift  coefficient  . — 

Lift  curve  slope  . .  . .  . .  . .  rad-1 

Maximum  lift  coefficient  . .  . .  . .  — 

Aerodynamic  drag  . N  (lb) 

Total  energy  . J  (ft  lb) 

Specific  energy,  energy  height,  energy  state  .  m  (ft) 

Engine  net  thrust . N  (lb) 

Non-dimensionaiizing  force  . N  (lb) 

Acceleration  due  to  gravity  .  9-80665  m/s2  (32-17405  ft/s2) 

Geopotential  altitude . m  (ft) 

Pressure  altitude . m  (ft) 

Pressure  altitude  increment  . m  (fi) 

Aerodynamic  lift . N  (lb) 

Aircraft  mass  . kg  (slug) 

Mach  number  .  . .  — 

Load  factor  normal  to  aircraft  in  plane  of  sym¬ 
metry  . .  . .  . .  , .  . .  . .  — 

Number  of  turns  at  optimum  conditions  . .  — 

Co-ordinate  origin  . — 

Ambient  pressure .  Pa  (lb/ft2) 

Specific  excess  power,  energy  rate  . .  . .  m/s  (ft/s) 

Turn  radius  . m  (ft) 

Range  . km  (n.m.) 

Range  increment . km  (n.m.) 

Gas  constant  for  air  .  287-053  J/kg.K 

(3089-78  ft2/Ks2) 

Aircraft  reference  area . m2  (ft2) 
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Time  . 

. .  s 

At 

Time  increment . 

. .  s 
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Temperature  . 

..  K 
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Flowrate  scaling  factor . 

1/3600  Ib.hi 

V 

Velocity . 

. .  m/s  (ft/s) 

W 

Aircraft  weight . 

. .  N  (lb) 

AW 

Weight  loss  due  to  fuel  usage 

. .  N  (lb) 

us 

Fuel  flowrate  . 

. .  kg/s  (lb/hr) 

x,  y,  z 

Co-ordinate  axes . 

. .  — 
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Body  incidence  to  flight  path 

rad  (deg) 
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CL 

Effective  angle  of  attack  (--*  \-c) 

. .  rad  (deg) 
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Flight  path  elevation  angle 

. .  rad  (deg) 

y»ir 

Specific  heat  ratio  for  air 

..  1-4 
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Error  function  to  be  solved  for  hp 

. .  m  (ft) 
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Temperature  lapse  rate 

. .  K/m  (K/ft) 

A 

Wing  sweep  angle 

. .  rad  (deg) 

a 

Thrust  incidence  to  body  axes 

. .  rad  (deg) 

X 

Energy  rate  function  . 

. .  — 

<t> 

Roll  angle . 

rad  (deg) 

<P 

Azimuth  angle 

. .  rad  (deg) 

to 

Turn  rate  . 

. .  deg/s 
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Definition 

a 

Ambient  conditions 

av 

Average  value  between  two  energy  states 

b 

Base  or  reference  level 

B 

Body-axis  co-ordinate 

E 

Earth-axis  co-ordinate 
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Iteration  count 

SL 

Sea  level  conditions 
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Trimmed  airciaft 

W 

Wind-axis  co-ordinates 

1,2 

Intermediate  co-ordinate  transformation 

Math 

symbols 

Definition 
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Differentiation  with  respect  to  time 

x’ 

Differentiation  with  respect  to  altitude 

ln(x) 

Natural  logarithm 

abbreviations  and  acronyms 

CAS  Calibrated  airspeed 

Cig.  Centre  of  gravity 

CPU  Central  processing  unit 

deg  Degree 

kt  Knot 

MAC  Mean  aerodynamic  chord 

MMD  Maximum  manoeuvre  diagram 

n.m.  Nautical  mile 

SI  Systeme  International  d’ Unites 

TAS  True  airspeed 
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Program  Names 

AIRCRAFT  Generic  name  of  main  programs  for  each  aircraft 
ANY  Member  of  AIRCRAFT  family  for  data  reprocessing 

BATCON  ARL  system  for  batch  program  operations 

OPTFIT  B-spline  curve-fitting  program 

P2  Plotting  and  data  reorganisation  program 

P4  Contour  plotting  program 

PLOTQ  System  program  for  queuing  plotter  files 

SURFM  B-spline  data  verification  program 


System  Filename  Extension 

.CTL  BATCON  control  file 

.EXE  Absolute  binary  core  image  of  program 

.FOR  FORTRAN  program  file 

.LOG  BATCON  progress  file  written  in  user’s  disk  area 

.LST  Line  printer  file  deleted  when  processed 

.REL  Relocatable  binary  file 


User  Filenames 

P2  P2A  P2DIFF  DUM  P4 


User  Filename  Extensions 
.OPT  .CON  .PLT 


User  Libraries 

PI  P1LIB  P24LIB  GRAF1C  EXTRAS 


Routines  in  Program  AIRCRAFT 

AERO  BLOCK  DATA  THRUST  TR1MCL 

Routines  in  Library  PI 

ALTIT  B1NIN  IDENT  MAXMAN  MONSEP  PARAMS 

PIOUT  PIOUTA  ROMIN  SEP  TABLE 

Routines  in  Library  PI  LIB 
SURF 


PI1N 


Routines  in  Program  P2 

GRID  INMMD  PSCON  PSDIFF  P2IN  RATE1  RATE2  R2DIFF 

Routines  in  Program  P4 

OUTXT  P  PLOTD  P4MAIN 

Routines  in  Library  P24L1B 
LINE  PLOT 

Routines  in  Library  GRAF1C 
CONI  SMOOTH 

Variables  Common  to  PI,  P2  and  P4 


B,  C 

Reserved  storage  arrays  of  equivalenced  variables 

DELES,  DELXM,  DELGN 

Increments  of  energy  state,  Mach  number  and  load 

ES0,  XM0,  GN0 

Initial  grid  values  of  these  variables 

NES,  NXM,  NGN 

Number  of  grid  values  of  these  variables 

I  UN  ITS,  IPSTYP 

Flags  for  types  of  units  and  energy  rate  variables 

PLA 

Power  setting 

Variables  in  Library  PI 

ALB 

Mass  scale  factor  (2-2046225  lb/kg) 

\M 

Length  scale  factor  (0-3048  m/ft) 

AWF 

Fuel  flow  rate  scale  factor  (3600  ALB  lb.s/kg.h) 

HP 

Pressure  altitude 

1ERR 

Error  flag 

IGRID 

Height  variable  flag 

IOPT 

Optimized  grid  flag 

IOUT 

File  output  flag 

OMTAB 

Turn  rate  output  vector 

PSTAB 

Energy  rate  output  vector 

Variables  in  Program  P2 

ICREAT 

On-line  data  generation  flag 

IONDSK. 

Differential  MMD  disk  file  flag 

IOPT 

Processing  option  flag 

NAMOUT 

Output  file  name 

OMEGA 

Turn  rate 

PS 

Energy  rate 

PS0,  DELPS,  NPS 

Energy  rate  grid  definition 

WORK 

Working  storage  vector 

Variable s  In  Program  P4 

CON  LAB,  LABFLG 

Contour  label  flag 

IC 

Running  counter  of  differential  MMD  boundary  points 

1COUNT 

Final  count  of  points  in  each  boundary 

iCHK 

Order  flag  for  x-co-ordinates  in  boundary  vector 

IDATA 

Data  set  leap  counter 

1DEF 

Default  contour  texture  flag 

I  DOT 

Current  contour  ‘exture  flag 

1DOT0 

Texture  flag  for  zero  contour 

1GES 

Energy  state  contour  flag 

IMMD 

Turn  rate  boundary  flag 

NDATA 

Number  of  data  sets 

NLEV0 

Level  number  of  zero  contour 

XPT,  YPT 

Plotter  co-ordinates  transmitted  by  CONT 

XTAB,  YTAB 

Vectors  of  differential  MMD  boundary  co-ordinates 

YLM1N 

Minimum  y  co-ordinate  on  MMD  boundaries 

Z,  ZG 

Working  storage  vectots 
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1.  INTRODUCTION 

With  the  increasing  complexity  of  modern  combat  aircraft  and  the  resulting  extension  of 
their  capabilities,  the  need  has  arisen  for  techniques  to  evaluate  the  relative  combat  performance 
of  new  designs,  und  to  compare  these  designs  against  existing  aircraft  in  service. 

This  report  and  two  companion  reports1**  dsecribe  techniques  developed  at  ARL,  which 
were  considered  to  meet  the  current  needs  of  the  Royal  Australian  Air  Force.  This  report 
describes  a  suite  of  FORTRAN  IV  computer  programs  developed  for  the  computation  and 
presentation  of  data  used  in  evaluating  combat  aircraft  performance  using  the  energy  manoeuvr¬ 
ability  theories3-8  which  have  been  developed  in  the  past  two  decades. 

The  suite  of  programs  produces,  for  a  variety  of  possible  aircraft  configurations,  tabulated 
and  plotted  data  in  either  metric  or  Imperial  units,  which  describe  aircraft  combat  performance 
in  energy  manouvrability  terms.  For  each  aircraft  to  be  evaluated,  the  user  must  supply  sub¬ 
routines  which  calculate  simplified  propulsion  and  aerodynamic  characteristics. 

The  mathematical  basis  for  the  calculations  is  described  in  Chapter  2.  Chapter  3  presents 
the  requirements  and  capabilities  of  the  suite  as  a  whole. 

Chapter  4  describes  the  main  calculation  program,  given  the  generic  name  AIRCRAFT 
in  subsequent  pages.  Chapter  5  describes  the  output  program  P2,  which  provides  plotted  output, 
as  well  as  input  for  the  contour  plotting  program  P4  described  in  Chapter  6. 

Chapters  7,  8  and  9  present  a  self-contained  user's  guide  with  full  instructions  on  how  to 
operate  the  three  programs.  Typical  outputs  are  described  in  Chapter  10. 

Chapter  1 1  concludes  with  some  suggestions  of  possible  extensions  to  the  suite’s  capabilities. 


2.  MATHEMATICAL  PRELIMINARIES 

A  complete  description  of  the  theoretical  background  for  this  suite  of  programs  is  given  in 
Reference  I,  supplemented  by  the  published  literature  of  References  3-6.  The  basic  equations 
are  given  below  without  derivation. 

Where  numerical  methods  are  available  in  the  published  literature,  they  are  quoted  without 
derivation.  An  iterative  technique  developed  by  the  author  is  the  subject  of  an  appendix,  and 
the  results  are  quoted  in  the  body  of  the  text. 


2.1  Equations  of  Motion 


The  equations  of  motion  ir.  flight  path  axes  as  shown  in  Figures  1  and  2,  for  an  aircraft 
in  general  turning  flight  are:1 


mV  —  Fd  cos  W sin  y  (2.1) 

mV  =  (— y  cos  <f>— ificos  y  sin  </>)  —  —  Fn  sin  a'—  L+W'cos  y  cos  <f>,  (2.2) 

mV  =  (  —  y  sin  cos  y  cos  <f>)  —  W cos  y  sin  <£.  (2.3) 

where  m  =  aircraft  mass,  kg  (slug), 

W  =  aircraft  weight,  N  (lb), 

V  -  true  airspeed  along  the  flight  path,  m/s  (ft/s), 
g  =  gravity  constant,  9*80665  m/s2  (32*17405  ft/s2), 

Fn  =  engine  net  thrust,  N  (lb), 

*'  ->  »+q  -  effective  angle  of  attack,  rad. 
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Fig.  1  Earth-and  wind-axes  systems 


Fig.  2  Force  system 
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D  "  aerodynamic  drag,  N  (lb), 

Z.  •  aerodynamic  lift,  N  (lb), 

$  <m  azimuth  angle,  rad, 
y  «  elevation  angle,  rad, 

<f>  «  roll  angle,  rad. 

In  deriving  these  equations,  it  is  assumed  that 

(a)  the  earth  is  flat; 

(b)  acceleration  due  io  gravity  is  constant; 

(c)  the  aircraft  is  considered  as  a  point  mass; 

(d)  atmospheric  parameters  follow  standard  laws; 

(e)  the  velocity  vector,  aerodynamic  forces  and  net  thrust  are  co-planar  in  the  aircraft  frame 
of  symmetry. 


2.2  Energy  Relationships 

Using  the  fundamental  relation  that  energy  state  (is),  or  specific  energy  is  given  by 

Ek-EIH'-.k+V*l2g,  (2.4) 

the  equations  of  motion  for  High*  in  a  horizontal  plane  (y  =  y  -  h  -  0)  can  be  expressed  as 


Ps  «  v~  m  V(Fn  COS  a~D)IW, 
g  dt 

and 

(2.5) 

n  =  1/cos  <f>  -  (Fn  sin  *'+£)/ H-', 

(2.6) 

where  Equation  (2.5)  is  a  restatement  of  Equation  (2.1)  and  Equation  (2.6) 
of  Equations  (2.2)  and  (2.3).  In  these  equations 

Ps  -  dEs/dt  =  energy  rate,  m/s  (ft/s) 

h  ■=  geopotential  height,  m  (ft), 

n  —  load  factor  normal  to  the  aircraft  in  the  plane  of  symmetry. 

In  addition  rate  of  turn  (^)  is  denoted  by  <u,  and  it  is  readily  shown  that 

is  a  combination 

«*  =  g(n*-l)a*IV, 

and 

(2.7) 

•n 

1! 

e 

(2.8) 

where  r  is  the  turn  radius,  m  (ft). 

For  flight  in  a  vertical  plane,  Equations  (2.1)  and  (2.2)  reduce  to 

Pa  -  (~  ■ +  ~j  =  V(Fn  cos «'  - D )/ fV,  (2.9) 


Vy/g-\- cos  y  =  (Fn  sin  *'+L)/IV  =  n. 

(2.10) 

Energy  state  methods  enable  climb  schedules  to  be  estimated  using  the  calculus  of  vari- 

ations7  based  on  Pa  or  some  function  of  Pa : 

<i  *  Pa, 

(2.11) 

X2  =  Pa/wt, 

(2.12) 

ys  =  FgF/wf, 

(2.13) 

3 


where  wr  =*  fuel  flow  rate,  kg/s  (Ib/hr). 

Integration  of  the  maxima  of  these  quantities,  viz: 


t  =a 


m  = 


m/R  = 


I 

I 

I 


Eh 

(1/Xl.m»x)<f£s, 

E»i 

(2.14) 

Eh 

(l/X2,max)</jEs, 

£si 

(2.15) 

E$t 

(1/X3,max)d.£s> 

Eh 

(2.15) 

results  in  approximate  profiles  for  minimum  time,  minimum  fuel  or  maximum  range  respectively. 

Apart  from  P$,  the  energy  functions  X‘  have  no  unique  names.  Hence,  in  the  program 
described  in  later  chapters,  an  output  variable  xt  is  referred  to  either  as  the  “energy  parameter”, 
or  simply  as  the  energy  rate,  with  the  implication  that  any  of  the  three  energy  functions  may  be 
referred  to. 


2.3  Atmospheric  Relationships 

The  algorithms  used  to  determine  atmospheric  pressure,  ambient  temperature  and  geo¬ 
potential  height  have  been  presented  fully  elsewhere,8  and  only  the  essential  details  are  given  here. 

Two  atmospheres,  the  ICAO  Standard  Atmosphere  and  the  AR.DU  Tropical  Atmosphere 
are  provided  as  atmosphere  models.  They  are  approximated  to  by  atmospheric  layers  with  con¬ 
stant  temperature  lapse  rates,  and  are  defined  by  values  of  ambient  temperature  and  pressure 
height  at  the  points  where  these  lapse  rates  change.  The  atmospheric  pressure  is  defined  by  the 
standard  atmosphere  relationships 

P  =  Wb/(7b  f  A  *  0) 

or  (2.17) 

P  =  P*  exp{  — gA/jp/(/?air7h)},  A  .  Oj 

where  P  -  atmospheric  pressure,  Pa  (lb/ft2) 

Mp  =  pressure  height  in  linear  segment  =  Ap — /rpt»,  m  (ft) 
ftp  —  pressure  height,  m  (ft), 

hp\>  —  pressure  height  of  base  of  linear  segment,  m  (ft), 

Pb  =  atmospheric  pressure  at  base  of  linear  segment,  Pa  (lb/ft2), 

7t,  =  ambient  temperature  at  base  of  segment,  K, 

A  =  temperature  lapse  rate,  K/m  (K/ft), 

R* ir  =  gas  constant  for  air,  287  055  J/kg.K  (3089  •  78  ft2/K.s2). 

Atmospheric  temperature  “7”’  is  given  by 

T  =  Tb+AA/ip,  K.  (2.18) 

In  the  ARDU  tropical  atmosphere,  geopotential  height,  required  for  performance  calcula¬ 
tion,  is  obtained  by  integration  over  all  layers  up  to  the  given  pressure  height 

h  =  "{.'hi,  (2.19) 

i - 1 

where 


Ai  =  A{A*/tp+(r&b -7bA./A)ln(l  +Mp/7b)},  A  *  0 


hi  —  fip(7'»b+A*AP/2)/7’b, 


A  =  0 
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and  h  »  geopotential  height,  m  (ft), 

hi  <*  geopotential  height  increment  in  ith  layer,  m  (ft), 

A*  »  lapse  rate  in  ambient  atmosphere,  K/m  (K/fl), 

r»b  ~  temperature  at  base  of  layer  in  ambient  atmosphere,  K, 

n  «  total  no.  of  layers  to  the  given  pressure  height  in  both  standard  and  non-standard 
atmospheies. 

In  the  ICAO  Standard  Atmosphere  geopotential  height  and  pressure  altitude  are  identical. 
Given  pressure,  altitude  and  Mach  number,  and  the  speed  relations 

V  =  Ma  (2.21) 


a  -  (y.ir/l.ifr.)0-6. 


(2.22) 


where  M  =  Mach  number, 

a  «■  sonic  speed,  m/s  (ft/s), 

y»ir  *=  specific  heat  ratio  for  air,  =1-4. 

Equation  (2.4)  may  be  used  to  determine  energy  state. 

However,  to  determine  pressure  altitude,  given  energy  state  and  Mach  number,  the  inter¬ 
dependence  of  geopotential  height,  pressure  height,  temperature  and  Mach  number  requires 
the  iterative  solution  of  Equation  (2.4),  which  is  non-linear  in  Ap  and  hence  in  h.  Appendix  1 
develops  the  solution  using  Newton's  method,  the  result  being  given  below: 


where 


Ap,r+i  —  Ap,i-|-5Ap,<, 


• .  _  h+kT%— £s~| 

p,<  ~  L.  ra/r-f*Aa  Jbp-Ap/ 

k  «  yRnrM2l(2g), 

Tft  =  ambient  temperature,  K 


(2.23) 

(2.24) 

(2.25) 


and  i  denotes  the  r'th  iteration. 

Iteration  continues  until  is  less  than  0-75  m  (0-25  ft).  A  suitable  initial  estimate  of 
Ap  is : 


AP0  =  £s-a8L2A*2/(2g),  (2.26) 

where  asi.  =  sonic  speed  at  sea  level,  340  30  m/s  (1 1 16  46  ft/s). 


The  final  atmosphere-related  quantity  required  is  the  non-dimensionalising  force  “End”, 
which  is  defined  as : 

End  =  O-Sy.irEA^S,  N  (lb),  (2.27) 

where  S  =  aircraft  reference  area,  m2  (ft2). 


3.  REQUIREMENTS  AND  CAPABILITIES 

The  object  of  the  suite  of  programs  for  air  combat  performance  estimation  is  to  produce 
tabular  and  plotted  data  for  accurate  assessment  of  the  performance  and  manoeuvrability  of 
of  any  aircraft  operating  in  a  wide  range  of  weapon/store  configurations. 
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In  particular,  the  programs  produce: 

(a)  specific  excess  power  ( Pa )  plots  on  a  height-Mach  number  grid  for  a  range  of  sustained 
normal  loadings; 

( b )  turn  rate  as  a  function  of  Pa  at  various  height  and  Mach  number  combinations; 

(c)  optimum  turn  rate  plots  on  a  specific  excess  power— energy  state  (Pa— £«)  grid; 

(d)  differential  plots  of  Ps  for  two  aircraft  on  a  height-Mach  number  grid  for  a  range  of 
sustained  normal  loadings;  and 

(<*)  differential  plots  of  optimum  turn  rate  on  a  Ps—E&  grid. 

The  tabular  and  plotted  output  may  also  be  used  to  assist  in  the  production  of  fuel/distance 
diagrams  used  in  a  variety  of  range  calculations.  Data  which  may  be  obtained  from  the  outputs 
include: 

(a)  estimated  fuel  used  and  range  achieved  when  flying  profiles  for  minimum  time,  minimum 
fuel  or  maximum  range; 

( b )  estimated  fuel  used  and  combat  conditions  for  optimum  sustained  manoeuvres; 

(c)  estimated  fuel  used  and  range  achieved  during  accelerations  using  maximum  power; 

(d)  estimated  maximum  speeds  at  constant  altitude; 

(e)  estimated  “corner  velocity’  ’and  sustained  turn  boundaries. 

Finally,  in  order  to  provide  a  concise  picture  of  combat  performance  throughout  the  air¬ 
craft’s  operating  regime,  plotted  output  can  be  produced  on  a  turn  rate/Ps  versus  Mach  number 
grid  for  a  range  of  load  factors  for  any  choice  of  altitudes  or  energy  states. 

The  suite  consists  of  three  main  programs,  whose  interrelations  are  shown  in  Figure  3. 

AIRCRAFT  is  the  generic  name  given  to  the  family  of  programs  which  represents  par¬ 
ticular  aircraft.  An  AIRCRAFT  program  includes  user-defined  routines  which  calculate  propul¬ 
sion  and  aerodynamic  parameters  using  thrust  and  drag  data  files.  The  program  uses  input  data 
describing  the  aircraft  configuration,  and  produces  tabular  output  for  line-printer  listing,  as  well 
as  unformatted  data  files  for  input  to  program  P2.  Program  ANY  is  a  particular  member  of 
the  AIRCRAFT  family  used  to  reprocess  published  performance  graphs  or  unformatted  data 
files,  producing  listings  and  further  input  files  for  program  P2. 

P2,  the  second  member  of  the  suite  processes  the  data  provided  by  program  AIRCRAFT. 
For  some  requirements,  disk  files  are  produced  which  may  be  immediately  plotted  off-line.  For 
other  requirements  program  P2  is  a  pre-processor  for  the  contour  plotting  program  P4,  re¬ 
arranging  data  arrays  and  performing  elementary  operations  such  as  are  required  for  comparing 
two  aircraft. 

Program  P4  is  a  contour  plotting  program  producing  only  files  for  off-line  plotting. 

Associated  with  these  three  programs  are  several  subroutine  libraries  used  in  the  loading 
phase  of  program  operation,  containing  input-output,  curve-fitting,  optimisation  and  plotting 
routines. 

Concise  descriptions  of  the  three  programs  and  the  required  subroutine  libraries  are  given  in 
Table  1.  The  programs  and  libraries  are  described  fully  in  Chapters  4  to  9,  together  with  instruc¬ 
tions  for  using  the  subroutine  libraries  in  the  loading  process. 


4.  PROGRAM  “AIRCRAFT”  DESCRIPTION 
4.1  Program  Structure 

Program  AIRCRAFT  is  the  starting  point  for  all  combat  performance  calculations.  It 
includes  aerodynamic  and  propulsion  routines  peculiar  to  each  aircraft  under  assessment,  and 
calculates  combat  performance  parameters,  such  as  energy  rate  and  turn  rate,  for  given  values  of 
energy  state,  Mach  number  and  load  factor  and  for  specific  aircraft  configurations. 

Output  files  may  be  produced  for  off-line  printing,  containing  full  details  of  all  relevant 
parameters,  or  for  input  to  subsequent  programs  which  prepare  graphic  presentations  of  the 
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Fig.  3  Suite  information  flow 


TABLE  1 

Library  Descriptions 


Name 

Function 

Brief  description 

AIRCRAFT 

Program 

library 

Contains  calling  program,  thrust,  drag  and  data  subpro¬ 
grams;  different  for  each  aircraft. 

P2 

Program 

library 

Processes  output  produced  by  “AIRCRAFT”,  and  pro¬ 
duces  files  for  plotting,  or  data  files  for  input  to  program 
P4. 

P4 

Program 

library 

Processes  data  produced  by  program  P2  and  produces  con¬ 
tour  plots. 

PI 

Subprogram 

library 

Processes  all  input  options  and  produces  text  output  and 
data  output  for  input  to  program  P2;  library  for  program 
“AIRCRAFT”. 

P1LIB 

Subprogram 

library 

Contains  curve-fitting  and  optimisation  routines;  library  for 
program  “AIRCRAFT”. 

P24LIB 

Subprogram 

library 

Contains  system  routines  for  processing  plot  instructions; 
library  for  programs  P2  and  P4. 

GRAF1C 

EXTRAS 

External 

subprogram 

libraries 

External  libraries  containing  contour  plotting  routines. 

TABLE  2 

Program  and  Subroutine  Libraries 


Library  name 

1  . '  " 

Routines 

AIRCRAFT 

MAIN.  AERO  BLOCK.  DATA  THRUST 
(and  any  routines  called  by  AERO  and  THRUST) 

PI 

PI1N  BADINP  TABLE  ATMOS  1NTRP  HEIGHT  HTRUE 
1DENT  PIOUT  B1NIN  ALTIT  PARAMS  MAXMAN  SEP 
PIOUTA  MONSEP 

P2 

MAIN.  GRID  INLAB  1NMMD  MMP  PLTLAB  PSCON 

PSDIFF  P2IN  RATE1  RATE2  R2DIFF  UNITS  WRLAB 

P4 

MAIN.  OUTXT  P  PLOTD  P4MAIN 

PI  LIB 

SURF  CUBICS  CHECKD  SPDER3  ROM  IN 

P24LIB 

INTRP  BADINP  AXIS  NUMBER  LINE  SYMBOL  INLAB 
PLTLAB  UNITS 

GRAFIC 

Routines  required:  CONT  DIAG  SMOOTH  REAL1N 

EXTRAS 

Routine  required:  PROMPT 

. . .  "un 
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printed  output.  Program  AIRCRAFT  user’s  guide  (Chapter  7)  gives  an  example  terminal  input. 
Sample  output  from  co-ordinated  running  of  programs  AIRCRAFT,  P2  and  P4  is  given  in 
Chapter  10. 

An  AIRCRAFT  program  consists  of  routines  in  program  library  AIRCRAFT,  together 
with  routines  from  subroutine  libraries  PI  and  PI  LIB.  The  routines  which  make  up  these  libraries 
are  listed  in  Table  2,  and  briefly  described  in  Appendices  2  to  4. 

The  division  into  program  and  subroutine  libraries  is  required  because  the  program  library 
is  different  for  each  aircraft,  but  the  service  routines  in  the  subroutine  libraries  are  unique.  The 
further  division  of  subroutine  libraries  into  two  parts.  Pi  and  PI  LIB,  arises  from  the  routines 
in  PILIB  having  their  origins  in  the  curve-fitting  program  OPTFIT. 

Appendix  S.  1  lists  the  calling  program,  designated  MAIN,  in  Table  2,  which,  together  with 
data  subroutine  BLOCK  DATA,  must  be  included  as  part  of  the  AIRCRAFT  program  library 
for  each  aircraft. 

The  two  remaining  subroutines  which  make  up  the  AIRCRAFT  program  library,  THRUST 
and  AERO,  are  aircraft  dependent  and  are  fully  described  in  Section  4.5.  Sample  routines  are 
given  in  Appendices  5.2  and  5.3.  A  standardised  form  for  aerodynamic  and  thrust  data  files  is 
described  in  Section  4.6 

Data  communication  between  subroutines  is  principally  by  means  of  labelled  COMMON 
areas  and  EQUIVALENCE  statements.  Appendix  6  lists  the  COMMON  areas  and  the  storage 
allocation  within  each  area.  As  indicated  in  that  appendix,  this  storage  allocation  scheme  is 
also  used  in  programs  P2  and  P4,  with  slight  modifications. 

A  flowchart  for  the  AIRCRAFT  program  structure  is  given  in  Figure  4,  and  it  can  be  seen 
that  program  operations  fall  logically  into  three  areas: 

(i)  input; 

(ii)  unoptimised  grid  calculation  and  output;  and 

(iii)  optimised  grid  calculation  and  output. 

Each  of  these  areas  is  discussed  below. 


4.2  AIRCRAFT  Program  Input 

The  details  of  the  conversational  input  routine  PI  IN  are  shown  in  Figure  5.  Subroutine 
PUN,  called  by  the  main  program,  performs  an  on-line  dialogue  with  the  user,  thereby  defining 
the  type  of  run.  The  user’s  guide  for  program  AIRCRAFT  given  in  Chapter  7  gives  exact  speci¬ 
fications  of  all  FORTRAN  input  parameters,  together  with  example  inputs. 

The  first  task  of  PUN  is  to  request  the  user  to  specify  his  output  units;  PUN  then  sets  scale 
factors  for  length,  mass  and  fuel  flowrate  accordingly. 

The  next  task  is  to  call  subroutine  1DENT,  which  requests  the  name  of  the  disk  file  contain¬ 
ing  the  following  aircraft  configuration  information: 

(i)  descriptive  name  of  the  aircraft ; 

(ii)  type  of  units  used  for  area  and  weight  in  the  file; 

(iii)  aircraft  wing  reference  area ; 

(iv)  gross  weight  for  the  current  role,  and  centre  of  gravity  position  (if  required); 

(v)  description  of  the  current  role; 

(vi)  file  names  of  the  thrust  and  drag  data  files; 

(vii)  tabular  data  defining  external  store  drag;  it  is  assumed  that  this  drag  may  be  defined 
in  terms  of  drag  dount  (equal  to  ACd.sx  10*)  versus  Mach  number,  from  which  linear 
interpolation  yields  the  external  store  drag  count  increment. 

These  data  are  then  read  from  the  aircraft  configuration  file,  followed  by  reading  of  the 
thrust  and  drag  data  files  with  calls  to  subroutine  BININ.  Control  then  returns  to  subroutine 
PUN  to  complete  the  conversational  dialogue. 
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spwpwpppwm 


Use  unity 
scale  factors 
as  defaults 


FIG.  5(a)  Program  AIRCRAFT  Conversational  Input  Flowchart 


Fig.  5{b)  Program  AIRCRAFT  Conversational  Input  Flowchart 
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The  following  run  data  are  next  supplied  in  response  to  program  prompts  from  routine 
PUN: 

(a)  The  type  of  calculation  required.  The  energy  parameter  may  be  one  of: 

(i)  Pa,  IPSTYP  «  1, 

(ii)  Pa/wt,  1PSTYP  *  2, 

(iii)  PaV/wt,  IPSTYP  -  3. 

( b )  The  height  variable  to  be  used  as  an  independent  variable. 

(c)  Grid  definitions  for  the  independent  variable  height,  Mach  number  and  load  factor. 

(</)  Engine  power  setting. 

(e)  Wing  sweep. 

(/)  Atmosphere  profile. 

(g)  Output  files  required. 


Grids  for  height,  Mach  number  and  load  factor  are  supplied  as  the  triplet;  initial  value, 
increment  and  number  of  points.  Two  height  grids  are  available — energy  state  (IGRID  -  1) 
and  pressure  altitude  (IGRID  =  2),  In  the  latter  case,  the  preset  grid  indicated  in  Figure  5  is 
available  as  a  default  (IV  <=  1),  giving  a  reasonably  fine  grid  for  most  purposes.  The  maximum 
value  for  the  Mach  number  grid  is  then  either  I  -  5  or  2  0  depending  on  military  or  afterburner 
throttle  setting. 

When  energy  state  is  requested  as  the  height  variable,  data  may  be  produced  on  an  optimum 
Ps-energy  state  grid  by  selecting  the  maximum  manoeuvre  option.  In  this  case,  a  Mach  number 
grid  is  not  needed,  since  Mach  number  will  be  chosen  by  the  program  to  optimise  the  energy 
parameter. 

Throttle  setting  is  coded  in  the  range  0  to  200  with  ranges  (0,  100)  and  (100,  200)  denoting 
partial  military  or  partial  afterburner  settings.  Note  that  availability  of  thrust  and  fuel  flow  data 
for  these  settings  depends  on  data  supplied  for  each  aircraft.  Wing  sweep  may  be  specified  if 
required  as  an  input  by  the  aerodynamic  routines. 

Where  possible,  each  conversational  input  line  is  checked  against  valid  limits,  and  the  input 
prompt  repeated  if  invalid  input  is  detected.  Valid  ranges  for  the  various  inputs  are  given  in 
Table  3. 


TABLE  3 


Program  AIRCRAFT  Input  Data  Velocity 


Paraineter 

Valid  condition 

Calculation  type 

Height  and  Mach  grid 
Load  factor  grid 

Throttle 

Wing  sweep 

Output  level 

I  «  IPSTYP  «  3 

(initial  >  0,  step  >  0,  0  <  number  50) 
(initial  >  1,  step  >  0, 0  <  number  <  50) 

0  <  PLA  «:  200 

0  ^  wing  sweep  <  80 

1  ^  10UT  ^  3 

4.3  Unoptimised  Grid  Calculations 

Unoptimised  grid  calculation  encompasses  all  calculation  other  than  that  of  the  maximum 
manoeuvre  diagram,  which  optimises  energy  rate  (Pa)  as  a  function  of  Mach  number.  Hence 
unoptimised  calculations  include  those  for  energy  rate  vs.  turn  rate,  energy  rate  contour  plots, 
and  turn  rate/energy  rate  vs.  Mach  number. 

Figure  6  presents  a  flowchart  for  these  unoptimised  calculations,  controlled  by  sub¬ 
routine  TABLE. 
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CALL  TABLE 


IOPT  -  0 

? 


Optimised  grid 
CALL  MAXMAN 


RETURN 


IGRID-  1 


Use  energy  state  as 
control  variable 
ES  -  ALT 


Use  pressure  altitude 
as  control  variable 
HP  -  ALT 


Mach  number  loop 
DO  910  IXM  -  1,  NXM 


Fig.  8(a)  Calculation  Flowchart  for  Unoptimised  Grids 
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height  — 
pressure  j 
altitude  f 


Calculate  tempei  atute 
and  geopotential  height 
explicitly  “ 

CALI  ATMOS 
CALL  HEIGHT 


height  = 
energy 
f  state 


Iterate  for  pressure  altitude 
temperature  and 
geopotential  height 

CALL  ALTIT 


Calculate  speed  and 
pressure  parameters 
initialise  a  — 
CALL  PAR  AMS 


IGRID  »  2 

7 


Yev, 


No  ' 

t _ 

Find  thrust  and 
fuel  flow  - 
CALL  THRUST 

wf  available  > 
if  needed? 

(WF  *  0  end  IPSTYP  *  1 


Es  =  h  +  v2/2g 

1 

1 

to 

908 


9999 


Error  text: 
'NO  WF  DATA' 


~~r~ 

c  st°°  ) 


Fig.  6(b)  Calculation  Flowchart  for  Unoptimised  Grids  (Com.) 
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Output  data  for  on*  M*ch 
numb*r  CALLPIOUT 


910 
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Height  and  Mach  number  loops  are  set  up  according  to  the  grid  data  supplied  by  subroutine 
Pi  IN.  Vectors  of  energy  parameter  ( * )  and  turn  rate  (a>)  arc  then  calculated  for  the  load  factor 
grid  specified.  At  the  end  of  each  load  factor  loop  these  vectors  are  transmitted  to  the  output 
files  before  passing  on  to  the  next  Mach  number  in  the  grid. 

Calculation  is  straightforward,  except  that  when  energy  state  is  the  height  variable,  an 
iterative  calculation  is  required  to  determine  pressure  altitude  and  hence  geopotential  height. 
Subroutine  ALTIT  performs  this  iteration  using  the  algorithm  developed  in  Appendix  1 . 

Output  for  printing  is  arranged  as  a  one  page  summary  for  each  height  requested.  On  this 
page,  each  call  to  PIOUT  produced  a  performance  summary  for  the  current  Mach  number, 
comprising 


Mach  number. 

True  airspeed, 

Energy  state  or  pressure  altitude, 
Geopotential  height 
Fuel  flow  rate, 

Energy  rate 


Turn  rate 


for  the  requested  load  factor  grid. 


These  data  are  also  written  on  the  alternative  output  file  as  input  to  program  P2.  In  this  output 
however,  a  minimum  of  explanatory  text  is  included. 

Provision  is  made  for  user  subroutines  THRUST  and  AERO  to  flag  error  conditions,  by 
returning  nonzero  values  for  the  parameter  IERR.  Error  conditions  encompass  thrust  or  aero¬ 
dynamic  calculations  outside  the  defined  data  envelopes,  inability  to  trim  the  aircraft,  and 
altitude  iterations  which  yield  pressure  altitudes  outside  the  range  sea  level  to  32  km  (104,987 
ft).  When  an  error  condition  is  detected,  an  immediate  jump  is  made  to  call  the  output  routine 
with  default  values  set  for  *i  and  o»i.  These  default  values  are  used  by  PIOUT  to  avoid  printing 
excessive  output  when  error  conditions  occur;  however,  default  values  are  included  on  the  P2 
input  file  and  subsequently  used  by  program  P2.  If  wr  is  required  for  calculating  x  and  no  data 
are  available,  an  error  text  is  typed  on  the  terminal  and  execution  ceases. 

Routine  PIOUT  converts  output  quantities  to  the  required  output  units;  the  only  exception 
to  this  conversion  rule  is  that  energy  parameter  x  is  calculated  directly  in  output  units  in  routine 
TABLE. 


4.4  Optimised  Grid  Calculation 

Optimised  grid  calculation  includes  calculation  of  data  files  to  produce  maximum  manoeuvre 
diagrams  (turn  rate  contours  plotted  on  a  maximum  energy  rate/energy  state  grid)  and  for 
producing  approximate  schedules  for  minimum  time,  minimum  fuel  or  maximum  range  climbs. 
In  addition,  conditions  for  optimum  sustained  turn  rate  are  obtained. 

Figure  7  presents  a  calculation  flowchart  for  the  controlling  routine,  subroutine  MAXMAN. 

Energy  state  and  load  factor  loops  are  set  up  according  to  the  grid  data  supplied  by  sub¬ 
routine  P11N.  The  Mach  number  is  then  determined,  for  each  combination  of  energy  state  and 
load  factor,  which  optimises  the  energy  parameter  x ■  The  optimisation  is  performed,  using  the 
direct  search  method  proposed  by  Rosenbrook"  and  implemented  by  Machura  and  Mulawa,,0iU 
by  subroutine  ROM1N.  The  energy  parameter  to  be  optimised  is  calculated  by  subroutine  SEP. 
Convergence  to  the  optimum  is  monitored  by  ROM1N  via  calls  to  subroutine  MONSEP.  Further 
details  concerning  the  use  of  ROMIN  and  its  communication  with  SEP  and  MONSEP  are  given 
in  Appendices  3  and  4. 

Subroutine  SEP  calculates  energy  parameter  given  energy  state  and  load  factor,  calling 
subroutines  ALTIT,  PARAMS,  THRUST  and  AERO  in  a  similar  manner  to  that  shown  in 
Figure  3.  Error  conditions  resulting  in  nonzero  values  of  IERR  assign  a  default  value  of  x  °f 
—9999-99,  imposing  an  effective  flight  envelope  constraint  on  the  energy  parameter. 

Two  features  of  the  optimisation  process  peculiar  to  the  current  application  deserve  special 
attention.  Firstly,  the  time  required  to  achieve  an  optimum  for  any  energy  state/load  factor 
combination  is  considerably  reduced  if  the  initial  Mach  number  estimate  is  carefully  chosen. 
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Fig.  7(a)  Calculation  Flowchart  for  Optimised  Grids 
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Fig.  7(b)  Calculation  Flowchart  for  Optimised  Grids  (Cont.) 
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Each  time  a  new  energy  state  is  considered,  the  Mach  number  estimate  is  found  from  an 
empirically-determined  function  of  energy  state  as  follows: 


Mi  =  £s/ 10000,  Ea  «  5000 

Mi  =  0  •  425 + 3£s/200000,  5000  <  £s  <  35000 

M  *  0-  lS+fs/50000,  £s  >  35000 


(3.1) 


where  Mi  =  initial  Mach  number  estimate, 

£s  ==  energy  state  in  feet. 

At  each  energy  state,  for  load  factors  other  than  the  first  of  the  grid,  the  final  Mach  number  from 
the  previous  optimisation  is  used  as  Mi  for  the  next  calculation.  A  limiting  value  of  (A/max — 0-05) 
is  applied  whenever  Mi  exceeds  A/m»x,  where  A/m»x  is  the  maximum  Mach  number  for  which 
data  is  available. 

Secondly,  Pa  variation  with  Mach  number  is  such  that  the  energy  parameter  x  way  have  two 
local  maxima,  corresponding  to  subsonic  and  supersonic  regimes.  At  low  values  of  £s  the  sub¬ 
sonic  peak  dominates,  while  at  higher  values  of  £s  the  supersonic  peak  gives  the  global  maximum. 
The  energy  state  at  which  the  transition  of  this  global  maximum  from  the  subsonic  to  the  super¬ 
sonic  peak  occurs  is  not  known  a  priori ;  hence  both  local  maxima  must  be  investigated.  This  is 
done  in  subroutine  MAXMAN  when  after!  urner  power  is  requested,  by  performing  two 
optimisations  for  each  load  factor  for  energy  states  between  25000  ft  and  50000  ft.  M i  for  the 
supersonic  optimum  is  0- 15  units  above  the  optimum  Mach  number  at  the  subsonic  peak.  The 
Pa  values  for  each  optimum  are  compared  and  the  largest  is  taken  as  the  global  maximum. 

Having  obtained  the  maximum  values  of  x  at  several  values  of  the  load  factor  grid,  the  con¬ 
ditions  for  the  best  sustained  turn  are  found  by  interpolating  load  factor,  turn  rate,  speed  and 
pressure  altitude  at  a  value  of  x  ~  0,  regarding  these  variables  as  a  function  of  x  as  load  factor 
increases. 

The  final  calculation  performed  by  MAXMAN  is  to  integrate  the  values  of  (1/x)  at  Ig  load 
factor  over  the  available  £s  range,  giving  an  approximate  climb  profile.  An  estimate  of  the  time 
taken  to  change  energy  state  (At)  is  given  using  Equation  (2.14): 

At  =  (l//>s,»v)A£s,  (3.2) 

where  A£s  is  the  difference  in  energy  states,  and  the  subscript  “av”  denotes  the  mean  of  values 
at  each  energy  state.  Pa  is  found  from  Equations  (2.1 1),  (2.12)  or  (2.13)  depending  on  the  value 
of  IPSTYP  (1,  2  or  3  respectively).  Then  the  fuel  increment  (AfF)  and  the  range  increment 
(A R)  are  found  from 

A  fP  =  w’t,*vAf,  (3.3) 

ana 

A R  —  V &V  cos  y»vAt,  (3.4) 

where 

sin  y»v  —  ^i).xv/(kr avAt).  (3.5) 


Time,  fuel  and  range  increments  are  summed  to  give  the  required  estimates  in  the  form  ofa  time- 
history. 

Routine  PIOUTA  controls  the  printing  of  optimised  grid  data,  arranged  as  a  block  of  data 
for  each  energy  state.  This  block  includes  optimum  energy  rate,  turn  rate,  Mach  number,  true 
airspeed,  pressure  aitutude  and  fuel  flow  rate,  tabulated  at  points  on  the  requested  load  factor 
grid.  It  also  includes  sustained  turn  rate  and  climb  performance  estimates.  Turn  rate  and  energy 
rate  data  are  also  transmitted  to  the  alternate  output  file  to  be  used  as  input  to  program  P2. 

Conversion  of  all  quantities  to  the  required  outputs  is  performed  within  routine  PIOUTA, 
all  other  internal  calculations  being  performed  in  Imperial  units. 


4.5  User-defined  Subroutines 

Only  subroutines  AERO  and  THRUST  need  to  be  defined  by  the  user,  although  the  standard 
main  program  and  the  BLOCK  DATA  subroutine  must  be  included  when  creating  the 
AIRCRAFT  program  library. 


■Mipil  mint  in.  mi 
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Subroutine  THRUST  returns  values  of  powerplant  thrust  and  fuel  flow  in  Imperial  units, 
using  input  values  of  altitude,  Mach  number,  thrust  setting  and,  if  required,  atmospheric 
quantities.  The  latter  data  is  required  if  thrust  data  is  available  either  for  alternative  atmospheres 
or  in  non-dimensional  form  with  corrections  for  atmospheric  variations.  Usually  only  thrust  for 
the  ICAO  Standard  Atmosphere  will  be  available. 

Thrust  and  fuel  flow  values  are  determined  by  curve  fitting,  interpolation,  or  thermodynamic 
calculation,  depending  on  which  form  of  data  is  available.  Data  required  are  as  follows: 

Thrust  «  fn(M,  hp,  throttle,  atmosphere), 

Fuel  flow  =  fn(M,  hp,  throttle,  atmosphere). 

A  flowchart  for  a  typical  THRUST  subroutine  is  given  in  Figure  8  and  a  sample  FORTRAN 
IV  thrust  routine  is  given  in  Appendix  5.2. 

Subroutine  AERO  returns  values  of  angle  of  attack  (degrees)  and  aerodynamic  drag 
(pounds).  Input  values  required  are  initial  angle  of  attack  estimate,  Mach  number,  thrust,  load 
factor,  c.g.  position  and  store  drag  table.  Altitude  will  be  required  where  aerodynamic  data 
makes  provision  for  any  Reynolds  number  corrections,  and  wing  sweep  angle  will  be  required 
where  the  aircraft  has  variable  geometry  features. 

The  subroutine  would  normally  be  written  in  two  sections.  The  first  section  uses  Newton’s 
method  to  solve  the  implicit  equation,  (2.6),  to  determine  trim  angle  of  attack  {at),  and  hence 
trimmed  lift  coefficient  (Clt).  If  this  iteration  fails  to  converge,  or  trims  at  an  a  which  exceeds 
the  maximum  allowed  value  of  Clt,  an  error  flag  is  set  and  a  return  is  made.  Once  Clt  is  satis¬ 
factory,  curve  fitted  or  interpolated  data  is  used  to  determine  minimum  drag  coefficient  (CD.min). 
coefficient  of  drag-due-to-lift  (Cd.i)  and  store  drag  coefficient  (ACd.s).  The  net  aerodynamic 
drag  is  then  determined  by  summing  the  drag  components  and  multiplying  by  the  dimensional- 
ising  force,  Fnd- 

Data  required  for  subroutine  AERO  (or  its  equivalent)  are  as  follows: 

Clt  =/«(«,  M,  A), 

Cl.«  =  fn{ a.  A/,  A), 

Cl. max  **  fn{M,  hp.  A), 

CD.min  —  fn{M,  hp.  A), 

Cd.i  =/m(Clt,  M,  hp.  A), 

ACd.s  =  fn{M,  A,  store  configuration). 

A  flowchart  for  a  typical  AERO  subroutine  is  given  in  Figure  9,  and  a  sample  FORTRAN 
aerodynamics  routine  is  given  in  Appendix  5.3. 

Note  that  any  further  routines  called  by  THRUST  or  AERO  must  be  provided  in  the  AIRCRAFT 
program  library.  The  example  in  Appendix  5  requires  a  subroutine  TRIM  CL  to  determine  aircrqft 
trim.  The  amount  of  coding  obviously  depends  on  the  complexity  of  the  data  required  to  represent  a 
given  aircraft. 


4.6  Standardised  Data  Storage  Allocation 

The  method  of  representation  of  propulsion  and  aerodynamic  data  depends  on  the  format 
of  available  data.  Where  these  data  are  produced  by  readily  available  computer  programs,  these 
programs  can  be  easily  incorporated  in  program  AIRCRAFT  as  subroutines.  The  data  may  be 
in  the  form  of  equations  which  can  be  readily  coded  into  FORTRAN  statements.  Usually, 
however,  the  data  are  in  the  form  of  a  family  of  graphs,  requiring  either  curve  fitting  or  interpola¬ 
tion.  1  he  former  is  more  economical  in  storage,  and  is  preferred,  providing  efficient  algorithms 
are  available  for  evaluating  the  polynomial  expressions  using  the  curve-fitted  coefficients. 

Two  methods  of  polynomial  curve  fitting  are  readily  available  at  ARL.  The  first1*  uses 
Chebyshev  polynomials,  and  was  developed  for  use  with  aerodynamic  functions  of  up  to  three 
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Fig.  9  Typical  AERO  Subroutine  Flowchart 
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independent  variables.  The  alternative  method  using  B-splines  as  a  basis  was  developed  specific¬ 
ally  for  use  in  preparing  data  for  program  AIRCRAFT,  and  experience  has  shown  it  to  provide 
fast,  accurate  evaluation  using  an  economical  storage  allocation. 

The  method  based  on  B-splines  uses  the  algorithms  developed  by  Coxln  14  and  de  Boor15 
incorporated  into  a  program  called  OPTF1T,  which  produces  B-spline  coefficients  with  optimum 
placement  of  knots.  The  coefficients  are  written  on  disk  files  which  are  then  read  by  an  auxiliary 
program,  SURFM.  This  program  enables  checking  of  the  curve  fits  over  any  desired  domain  of 
the  independent  variables,  and  combines  several  files  cf  ASCII  data  into  a  single  binary  file  of 
propulsion  or  aerodynamic  data  for  use  with  program  AIRCRAFT'.  Program  library  PI  LIB 
contains  the  routine  SURF  and  associated  routines  needed  to  evaluate  the  B-spline  coefficients 
and  produce  the  desired  propulsion  or  aerodynamic  data. 

The  sample  THRUST  and  AF.RO  routines  presented  in  Appendix  5  use  data  stored  in  this 
manner.  The  binary  data  file  structure  used  to  represent  aerodynamic  and  thrust  data  in  arrays 
E  and  F  is  presented  in  Table  4. 


TABLE  4 

Binary  Data  File  Structure 

The  data  file  consists  of  a  sequence  of  groups  of  data,  each  representing  a  family  of  one  or 
more  curves  of  an  aerodynamic  or  propulsion  parameter  (military  power  thrust,  induced  drag 

coefficient,  etc.) 


Item 

Word 

Variable 

Description 

1 

1 

AK 

Number  of  words  in  the  current  group  of  data 

2 

2 

TITLE 

One  word  code  for  group  of  data 

3 

3 

NCURV 

Number  of  curves  in  family  (floating  point 

format) 

4 

4-> 

CURV(l)-» 

NCURV  values  of  independent  variable  >>  for 

(NCURV+3) 

CURV(NCURV) 

each  curve 

5 

(NCURV  (-4) 

K0 

Lower  bound  on  independent  variable  x 

6 

(NCURV+5) 

KM 

Upper  bound  on  independent  variable  x 

7 

(NCURV+6) 

NCAP, 

No.  of  arcs  representing yth  curve 

9 

(NCURV  +  7)-> 

K(l)> 

(NCAP— 1)  interior  knots  representing  y'th  curve 

(NCURV+6  + 
NCAP) 

K(NCAPr-l) 

9 

(NCURV+7+ 

Cl(l)-> 

(NCAP-i-3)  coefficients  representing /th  curve 

NCAP,)-> 
(NCURV+6  t- 

CI(NCAP,  +  3) 

2  x  NCAPy) 

Items  7,  8  and  9  are  repeated  for  each  of  the  curves,  giving  a  total  storage  requirement  for 

NCURV 

each  group  of  data  of  (4x  NCURV+2  £  NCAP(./)f5). 

j- i 

4.7  Program  ANY 

If  basic  aeropropulsion  data  for  an  aircraft  is  not  available,  but  comparative  combat  per¬ 
formance  is  still  required,  it  is  possible  to  use  published  performance  curves  of  energy  rate  against 
turn  rate  or  energy  rate  contours  on  a  height/Mach  number  grid. 

For  this  purpose  a  special  member  of  the  class  of  AIRCRAFT  programs  called  ANY  has 
been  compiled.  This  makes  use  of  the  unoptimised  grid  calculations  (Section  4.3),  replacing  the 
calculation  of  energy  parameter  using  thrust  and  drag  data  by  a  conversational  input.  No  input 
files  are  required,  all  data  being  supplied  via  the  user's  terminal.  Output  files  are  produced  as 
before,  with  one  file  intended  for  printing  and  the  other  as  input  for  program  P2.  An  additional 
feature  included  in  program  ANY  is  that  it  will  re-read  the  alternate  output  file  (possibly  with 
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corrections  required  by  on-line  typing  errors)  to  reproduce  the  printed  output  hie. 

Such  a  program  has  been  found  to  be  useful  for  a  variety  of  purposes: 

(a)  Published  turn  rate  curves  or  energy  rate  contours  may  be  replotted  on  a  more  suitable 
scale. 

(/>)  Data  files  can  be  produced  for  differential  contour  plots.  Data  from  aeropropulsion 
sources  can  be  readily  compared  with  that  from  published  graphs,  as  well  as  comparing 
one  set  oi  published  graphs  with  another. 

(c)  The  re-reading  feature  can  be  used  to  reproduce  formatted  printer  files  (provided  P2 
input  is  retained)  without  the  need  to  recalculate  the  complete  grid. 

( d )  Airspeed-altitude  data  can  be  produced  in  printed  or  plotted  form,  to  complement 
energy  rate  data  over  the  required  energy  state/Mach  number/load  factor  grid. 

A  similar  feature  for  replotting  and  comparing  optimised  grids  (maximum  manoeuvre 
diagrams)  is  incorporated  as  an  option  in  program  P2. 


5.  PROGRAM  P2  DESCRIPTION 

5.1  Program  Structure 

Program  P2  is  essentially  a  dual  purpose  program.  Or.  the  one  hand,  it  prepares  plotter 
output  for  turn  rate  diagrams  and  overview  plots  of  the  data  grid.  On  the  other  hand,  it  prepares 
data  files  for  the  contour  plotting  program  P4.  In  addition  assistance  is  given  in  calculating 
turning  endurance  for  given  fuel  quantities. 

The  input  to  the  program  consists  of  data  files  prepared  for  particular  aircraft  by  program 
AIRCRAFT,  and  additionally  in  the  case  of  maximum  manoeuvre  diagrams,  comparative 
aircraft  data  based  on  existing  diagrams  may  be  supplied  in  a  conversational  mode. 

The  modular  st.ucture  of  program  P2  may  be  seen  with  reference  to  Figure  10.  Subroutine 
P21N  opens  data  files  and  reads  identification  headers  and  then  converses  with  the  user  to 
determine  what  operations  to  perform  on  the  data.  This  is  indicated  by  parameter  IOPT,  whose 
range  of  seven  values  determines  the  subroutine  call  for  subsequent  processing.  The  modular 
structure  allows  for  easy  future  addition  of  other  data  processing  operations. 

The  actions  taken  by  the  input  routine  P2IN  and  the  various  subroutine  options  are  con¬ 
sidered  in  the  remainder  of  this  chapter. 

Brief  descriptions  of  all  user  subroutines  required  by  program  P2  are  given  in  Appendices 
7  and  8.  Examples  of  terminal  input  are  given  in  the  user’s  guide  for  program  P2  in  Chapter  8, 
and  examples  of  output  produced  by  co-ordinated  running  of  programs  AIRCRAFT.  P2  and 
P4  are  given  in  Chapter  10. 

5.2  Input  Operations  with  Subroutine  P2IN 

Subroutine  P21N  determines  the  option  parameter  IOPT  and  performs  various  preliminary 
data  input  operations,  depending  o.i  the  value  of  IOPT.  Figure  1 1  provides  a  flowchart  for  this 
subroutine. 

After  initial  settings  of  the  default  output  filename  and  recording  the  current  CPU  time,  an 
option  code  is  requested  to  determine  the  type  of  run.  The  seven  codes,  their  meanings,  and  output 
files  produced,  are  indicated  by  the  following  terminal  prompt: 

OPTIONS  ARE  (OUTPUT  FILENAME  IN  BRACKETS)  :  IOPT 

4A  PS  CONTOUR  PLOT.  (P2.CQN)  1 

4B  PS  VS  TURN  RATE  FOR  GIVEN  HEIGHT.  (P2.PLT)  2 

4C  HAXIHUH  MANEUVER  DIAGRAM  -  MHD.  (P2.0PT)  3 

4D  PS  DIFFERENTIAL  PLOT.  <P2.CON,P2A.CON,P2DIFF.CQN>  4 

4E  HMD  DIFFERENTIAL  PLOT.  (P2.0PT,P2A.0PT,P2DIFF.OPT)  5 

4F  MNP  CALCULATION  ASSISTANCE.  6 

47.  TURN  RATE,  PS  VS  MACH  FOR  GIVEN  HEIGHT.  (P2.PLT )  7 

**  ANY  OTHER  REPLY  PRODUCES  THIS  HELP  TEXT 
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if  MMP  calculation  assistance  is  requested,  an  immediate  return  to  the  main  program  is 
made,  since  all  input/output  operations  for  that  option  are  conversational,  on  the  user's  terminal. 

If  a  differential  MMD  is  requested  (IOPT  «  5),  a  further  test  is  made  if  the  two  files  of 
required  comparative  data  (filenames  P2.0PT  and  P2A.OPT)  already  exist  on  disk.  If  not, 
subroutine  P2IN  continues  in  order  to  produce  these  files.  If  “.OPT”  files  do  already  exist,  the 
flag  IONDSK  is  set,  and  comparison  filename  P2A.QPT  stored  in  NAMOUT  before  returning 
to  the  calling  program. 

If  any  type  of  MMD  is  requested  (IOPT  »  3  or  5)  and  on-line  data  is  to  be  supplied,  flag 
ICREAT  is  set,  and  routine  INMMD  is  called  to  accept  identifying  text  before  returning  to  the 
calling  program. 

For  remaining  options,  a  filename  of  output  from  program  AIRCRAFT  is  supplied  and 
header  data  identifying  the  aircraft  and  the  data  grids  is  read.  The  storage  requirement  is  checked 
against  the  8000  words  available  in  array  WORK  ,  and,  if  it  is  excessive,  art  error  text  is  typed  on 
the  user’s  terminal  and  execution  ceases.  If  sufficient  storage  is  available,  execution  continues  by 
assigning  text  headings  and  echoing  data  units  and  the  type  of  energy  parameter  on  the  terminal. 

If  plotter  output  is  to  be  produced  (IOPT  =»  2  or  7),  plot  file  P2.PLT  is  opened  before  return¬ 
ing  to  the  main  program. 


5.3  Subroutine  PSCON  tc  Produce  Energy  Rate  Contour  Data 

The  purpose  of  subroutine  PSCON  is  to  rearrange  the  grid  data  produced  by  program 
AIRCRAFT  and  produce  energy  rate  contour  output  compatible  with  the  requirements  of  the 
contour  plotting  program  P4. 

The  operation  of  the  subroutine  is  quite  straightforward,  and  can  be  described  without 
reference  to  flowcharts.  The  file  indicated  by  NAMOUT  is  opened  and  the  header  text  written. 
An  outer  loop  for  the  height  variable  and  an  inner  loop  for  Mach  number  are  set  up  and  Mach 
number  and  energy  state  are  read  from  the  input  file,  followed  by  energy  rate  data  over  the  full 
load  factor  grid.  The  matrix  of  energy  stated  on  the  Mach  number/altitude  grid  is  written  on 
the  output  file,  in  case  energy  state  contours  will  be  required. 

The  energy  rate  matrix  on  the  Mach  number/altitudc  grid  is  then  written  on  the  output  file 
for  each  of  the  requested  load  factors.  A  conversational  input  allows  the  user  to  select  the  load 
factors  fot  which  contour  plots  are  desired,  from  the  available  data.  , 

Load  factor  requests  are  stored,  so  that  if  PSCON  is  called  by  subroutine  PSD1FF  to  pro¬ 
duce  data  for  differential  contour  plots,  the  load  factor  grid  will  be  identical. 

At  the  end  of  the  load  factor  loop,  control  returns  to  the  calling  routine,  either  the  main 
program  or  subroutine  PSDIFF. 


5.4  Subroutine  RATE1  to  Produce  Turn  Rate  Data 

Subroutine  RATE1  is  called  when  IOPT  has  the  value  2  (option  4B).  Its  purpose  is  to  read 
the  data  file  produced  by  program  AIRCRAFT  and  produce  plots  of  energy  rate  against  turn 
rate.  The  plots  are  produced  by  joining  the  (turn  rate,  energy  rate)  co-ordinates  of  points  calcu¬ 
lated  at  small  increments  of  load  factor. 

A  flowchart  for  the  subroutine  is  given  in  Figure  12.  The  first  action  of  the  subroutine  is  to 
plot  a  header  text  on  the  plot  file  on  logical  unit  1,  followed  by  setting  of  storage  pointers.  Axis 
scales  and  energy  rate  axis  definition  are  provided  on  the  user’s  terminal.  An  “invalid  input”  text 
is  given  if  either  scale  is  zero  or  if  the  y-axis  length  is  not  in  the  range  0  to  26-7  cm  (10-5  in.). 

One  page  of  plots  is  produced  for  each  requested  height  variable  on  an  A4  size  area.  From 
the  sets  of  height  and  Mach  number  fata  on  the  data  input  file  (logical  unit  5),  particular  com¬ 
binations  of  height  and  Mach  number  may  be  selected  by  replying  to  the  prompts  as  follows: 


Reply 

Meaning 

YES 

. .  Plot  data  for  this  value 

NO 

. .  Skip  to  the  next  value 

ALL 

. .  Plot  data  for  all  remaining  curves  at  this  height 

END 

Terminal  execution  now. 
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Fig.  12(a)  Subroutine  RATE1  Flowchart 
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If  data  are  not  required  for  any  height,  the  input  Die  is  advanced  to  consider  the  next  height. 
Otherwise  axes,  A4  trim  edges  and  height  and  load  factor  labels  are  plotted  before  reading 
turn  rate  data  for  the  first  Mach  number.  If  ail  Mach  numbers  are  requested,  ptotting  is  continued 
without  further  Mach  number  prompts. 

If  data  is  not  requested  at  any  Mach  number,  the  program  reads  in  data  for  the  next  Mach 
number  before  giving  the  request  prompt.  For  each  Mach  number  requested,  scales  and  minimum 
values  of  the  x -  and  y-axes  are  stored  in  locations  required  by  subroutine  LINE  before  calling 
that  subroutine  to  join  the  (turn  rate,  energy  rate)  co-ordinates  with  straight  lines.  These  co¬ 
ordinates  are  at  fixed  increments  of  load  factor,  and  it  is  evident  that  reasonably  small  increments 
are  required  when  running  program  AIRCRAFT  (approximately  0  25 g)  to  produce  a  smooth 
curve.  Ig  increments  of  load  factor  are  indicated  on  the  curve,  together  with  a  label  identifying 
the  Mach  number  and  alternate  height  for  the  curve.  (“Alternate  height”  refers  to  whichever 
of  the  two  height  variables,  energy  state  or  pressure  altitude,  is  not  being  used  as  the  control¬ 
ling  grid.) 

After  all  Mach  numbers  for  the  given  height  have  been  considered,  the  pen  is  moved  to  a 
new  origin  for  the  next  value  of  height.  After  all  heights  have  been  considered,  control  returns 
to  the  main  program  to  terminate  execution. 

5.5  Subroutine  RATE2  to  Produce  MMD  Data 

The  purpose  of  subroutine  RATE2  is  to  produce  turn  rate  contour  data  output  on  an 
energy  state/optimum  energy  rate  grid  prior  to  processing  by  the  contour  plotting  program  P4. 
Usually  the  input  will  be  a  data  file  produced  by  program  AIRCRAFT  as  optimum  energy  rate 
data  on  an  energy  state/load  factor  grid;  the  data  is  rearranged  to  the  desired  output  format. 
Alternatively,  input  may  be  supplied  in  conversational  mode,  reproducing  published  MMD 
plots  on  the  required  MMD  grid. 

Figure  13  presents  a  flowchart  for  this  subroutine.  The  energy  rate  grid  for  the  MMD 
contour  plot  is  supplied  via  the  user’s  terminal  as  the  first  operation.  Definition  of  this  grid  is  not 
required  if  on-line  data  is  being  supplied  or  if  a  file  of  comparative  data  is  being  produced, 
since  the  grid  has  already  been  defined,  either  in  a  call  to  subroutine  INMMD,  or  in  a  previous 
call  to  the  subroutine  RATE2. 

Within  the  loop  for  each  energy  state  two  main  operations  occur.  These  are,  firstly,  to  con¬ 
struct  two  vectors  of  optimum  energy  rate  and  corresponding  turn  rate  data  as  load  factor 
increases  (if  on-line  data  is  being  supplied).  Secondly,  using  either  this  on-line  data  or  data  read 
from  a  disk  file,  to  interpolate  or  extrapolate  turn  rate  values  at  fixed  increments  of  optimum 
energy  rate. 

In  supplying  on-line  MMD  data,  pairs  of  (turn  rate,  optimum  energy  rate)  are  requested. 
With  turn  rate  increasing.  If  a  non-increasing  turn  rate  value  is  detected,  all  data  for  the  current 
energy  state  is  rejected  and  the  data  input  restarted.  Up  to  50  data  points  may  be  supplied,  the 
end  of  the  curve  being  indicated  by  a  negative  value  of  turn  rate. 

In  using  the  PSTAB  and  OMTAB  vectors  to  find  turn  rate  values  at  fixed  increments  of 
optimum  energy  rate  (Ps),  three  regions  are  detected.  If  data  limits  are  exceeded  (e.g.  lift  limit) 
at  large  negative  values  of  Ps,  default  values  for  turn  rate  (-10  0)  and  Ps  (—9999-99)  are 
stored  in  the  grid  matrix  WORK.  If  the  grid  value  of  Ps  is  above  the  zero  turn  rate  boundary, 
nominal  values  of  turn  rate  are  calculated  by  extrapolation.  These  (negative)  values  of  turn  rate 
are  physically  unrealisable,  but  are  useful  in  providing  a  smooth  zero  turn  rate  contour  on 
MMD  plots.  The  third  region  is  that  where  valid  turn  rate  data  is  calculated  by  linear  interpola¬ 
tion  in  the  vectors  OMTAB  and  PSTAB  at  the  grid  value  of  Ps. 

After  all  input  energy  states  have  been  considered,  the  turn  rate  matrix  is  written  on  the 
output  file  on  an  energy  state/optimum  energy  rate  grid,  for  use  as  input  to  the  contour  plotting 
program  P4. 

5.6  Subroutine  PSDIFF  to  Produce  Differential  Energy  Rate  Con.our  Data 

Subroutine  PSDIFF  is  called  by  the  main  program  when  IOPT  has  the  value  4  (option  4D). 
Data  files  produced  by  program  AIRCRAFT  are  read  in  order  to  produce  differential  energy 
rate  contour  data  for  input  to  program  P4. 
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Fig.  13{b)  Subroutine  RATE 2  Flowchart  (Cont) 
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A  flowchart  for  the  subroutine  is  given  in  Figure  14.  The  logic  is  straightforward,  but  input 
and  output  file  operations  need  some  comment.  Using  the  file  indicated  as  primary  input,  sub¬ 
routine  PSCON  is  called  to  produce  base  data  on  file  “P2.CON”.  These  input  and  output  files 
are  closed  and  the  comparison  file  is  opened  and  grid  data  read.  If  the  base  and  comparison 
grids  for  height,  Mach  number  or  load  factor,  or  flags  for  unit  systems  or  types  of  energy  para¬ 
meter  differ  in  any  way,  an  error  text  is  typed  and  execution  terminated. 

If  the  grids  are  identical,  comparison  contour  plot  data  is  written  on  file  P2A.CON  with  a 
second  call  to  subroutine  PSCON.  input  and  output  files  are  again  closed;  base  and  comparison 
data  files  opened  as  input  files,  and  file  “P2DIFF.CON”  opened  to  receive  the  differential 
contour  data.  Grid  data  is  transferred  to  the  output  file,  together  with  identifying  texts  for  both 
aircraft.  Energy  state  contour  data  is  written  for  the  required  Mach  number/height  grid  and 
input  files  positioned  to  read  energy  rate  data. 

Load  factor  and  height  loops  are  set  up,  to  read  vectors  of  energy  rate  data  at  each  height. 
These  vectors  are  subtracted,  and  the  difference  vector  written  on  the  output  file.  The  load 
factor  loop  continues  until  all  requested  load  factors  have  been  processed.  If  any  load  factors 
were  not  requested  for  plotting  during  calls  to  PSCON,  unit  4  will  run  out  of  data  before  the 
loop  terminates,  and  an  information  text  to  that  effect  is  typed  on  the  terminal. 

It  is  seen  that  using  the  option  for  differential  plots  results  in  three  data  files — one  for  each 
aircraft  and  the  differential  data,  any  of  which  may  then  be  processed  with  program  P4. 


5.7  Subroutine  R2DIFF  to  Produce  Differential  MMD  Data 

Subroutine  R2DIFF  is  called  by  the  main  program  when  IOPT  has  the  value  5  (option  4E). 
Data  files  produced  previously  may  be  read,  or  data  may  be  supplied  in  a  conversational  mode, 
in  order  to  produce  differential  maximum  manoeuvre  diagram  data  as  input  to  program  P4. 

A  flowchart  for  the  subroutine  is  given  in  Figure  15.  The  structure  is  similar  to  that  for 
subroutine  PSDIFF,  but  is  made  more  complex  by  the  provision  of  options  for  generating 
data  on-line  and  for  using  existing  “.OPT”  files  to  generate  a  differential  data  file.  The  sub¬ 
routine  is  roughly  in  two  phases:  generation  of  base  and  comparison  files  (Fig.  1 5a)  and 
generation  of  the  differential  file  (Fig.  1 5b). 

If  the  files  for  comparison  are  on  disk  when  R2D1FF  is  called  (lONDSK.  -  1),  the  first 
phase  is  skipped.  Otherwise  the  base  and  comparison  files  are  generated  with  calls  to  subroutine 
RATE2,  using  either  on-line  data  or  data  files  generated  by  program  AIRCRAFT.  Units  5  and 
6  are  used  as  input  and  output  units  respectively,  in  this  phase.  Energy  state  grids  of  base  and 
comparison  data  and  flags  for  unit  systems  and  types  of  energy  parameter  are  compared,  and,  if 
different,  an  error  text  is  typed  before  terminating  execution.  In  the  first  phase  there  is  no  need 
to  compare  energy  rate  grids  because  the  same  input  reply  has  been  used  for  both  sets  of  data. 
Base  and  comparison  data  are  written  on  files  "P2.0PT”  and  “P2A.OPT”  respectively. 

In  the  second  phase,  a  preliminary  file  of  differential  turn  rate  is  written  on  file  “DUM”, 
using  files  “P2.0PT”  and  "P2A.OPT”  as  logical  units  4  and  5  respectively.  Identifying  text 
headers  are  not  needed  on  this  preliminary  output,  and  input  files  are  positioned  accordingly  by 
the  loop  “Skip  file  headers".  If  independently  generated  “.OPT”  flies  are  used,  it  is  necessary 
to  compare  data  grid  definitions,  if  any  difference  appears,  an  error  text  is  typed  and  execution 
terminated. 

Turn  rate  vectors  for  each  energy  rate  level  are  then  read  from  the  input  files,  subtracted, 
and  the  result  written  on  the  output  file.  The  system  program  PIP  is  run  to  combine  the  input 
files  and  the  preliminary  output  file  into  the  final  output  file  “P2D1FF.OPT”.  This  file  now 
contains  contour  data  for  both  aircraft  and  their  comparison;  the  three  sets  of  data  are  required 
for  program  P4  to  generate  zero  turn  rate  (!g)  boundaries  as  well  as  differential  turn  rate 
contours. 

CPU  time  used  is  calculated  by  subroutine  R2DIFF  because  the  running  of  program  PIP 
bypasses  CPU  time  output  generated  by  the  FORTRAN  operating  system. 

5.8  Subroutine  MMP  to  Aid  in  Maximum  Manoeuvre  Persistence  Calculation 

Subroutine  MMP  is  called  by  the  main  program  when  IOPT  has  the  value  6  (option  4F). 
No  input  or  output  files  are  involved  as  all  communication  is  via  the  user’s  terminal.  When 
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Fig.  14  (a)  Subroutine  PSDIFF  Flowchart 
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d (hi  Subroutine  PS01FF  Flowchert  (Cont.) 


fuel/distance  diagrams  are  constructed  for  outward  and  return  legs  of  a  mission,®  the  resulting 
wedge  represents  the  diminishing  fuel  available  for  manoeuvres  as  range  from  base  increases. 
For  any  given  energy  state,  using  maximum  manoeuvre  output  from  program  AIRCRAFT 
the  user  can  determine  optimum  sustained  turn  rate  and  the  fuel  flow  rate  at  the  resulting  Mach 
number  and  altitude.  Subroutine  MMP  aids  in  calculating  the  number  of  turns  («t)  possible  at 
optimum  turn  rate  and  fuel  flow  rate,  given  the  available  fuel  quantity,  according  to 


/it 


l(V  .  W 

VVi'f 


(5.1) 


where  to  =  optimum  sustained  turn  rate  (deg/s), 

>t  f  =  fuel  flow  rate,  kg/s  (Ib/hr), 

W  =  available  fuel,  kg  (lb), 

v  -  unit  scale  factor  =  1/3600  (1  for  Imperial  units). 

The  routine  first  accepts  the  number  of  ranges  and  the  available  fuel  at  each  of  these  ranges. 
Then  for  a  given  energy  state,  optimum  turn  rate  and  fuel  flow  rate  are  provided,  and  Equation 
(5.1)  is  calculated  for  each  range.  New  energy  state  data  are  requested  and  calculation  repeated 
until  a  zero  value  of  energy  state  is  supplied,  when  execution  terminates. 

n  is  termed  the  maximum  manoeuvre  persistence,  and  may  be  plotted  either  against  energy 
state  at  fixed  ranges,  or  against  range  at  fixed  energy  states.  The  former  is  more  informative. 
Example  plots  are  given  in  Reference  2. 


5.9  Subroutine  GRID  to  Produce  Raw  Data  Overview  Plots 


Subroutine  GRID  is  called  when  IOPT  has  the  value  7  (option  4Z).  Its  purpose  is  to  read 
the  data  file  produced  by  program  AIRCRAFT,  and  produce  plots  of  turn  rale  and  energy  rate 
against  Mach  number  for  given  values  of  the  height  variable.  The  plots  are  produced  by  joining 
the  required  points,  calculated  at  small  increments  of  Mach  number.  These  plots  are  intended 
as  an  overview  of  the  data  grid,  and  consequently  are  produced  with  predetermined  scales.  The 
variable  ranges  obtained  using  these  scales  are  shown  in  Table  5. 

A  flowchart  for  the  subroutine  is  given  in  Figure  16.  The  first  action  is  to  plot  £.  header  text 
on  logical  unit  1,  followed  by  allocation  of  storage  constants.  From  the  grid  data  supplied  on  the 
data  input  file  (logical  urit  5),  particular  combinations  of  height  and  load  factor  may  be  selected 
by  replying  “YES”,  “NO”,  “ALL”  or  “END”,  as  indicated  in  Section  6.4.  If  data  is  not  to  be 


TABLE  5 

Ranges  of  Variables  Plotted  by  Routine  GRID 


3 


-75 

-150 


50 

100 


25 

50 


plotted  for  any  height,  the  input  file  is  advanced  to  the  next  height.  Otherwise,  axes  for  turn  rate 
and  energy  rate  are  plotted  in  the  y  direction,  and  Mach  number  axes  are  plotted  in  the  x  direc¬ 
tion,  together  with  a  height  label  for  each  pair  of  graphs.  For  each  Mach  number,  energy  rate 
and  turn  rate  vectors  are  read  into  the  array  WORK. 

At  this  point  all  the  data  for  the  plot  at  the  current  height  has  been  read  in.  The  remainder 
of  the  subroutine  is  concerned  with  plotting  the  data  in  a  different  grid  order  from  that  on  the 
input  file.  A  load  factor  loop  is  set  up  and  requests  for  plots  of  curves  at  selected  load  factors 
arc  supplied  via  the  user's  terminal.  For  each  load  factor  selected  pointers  K1  and  K2  are  cal¬ 
culated  to  indicate  storage  locations  corresponding  to  the  lowest  valid  Mach  number  and  plot 
scale  parameters.  These  parameters  (minimum  value  and  axis  scale)  and  then  stored  for  each  of 
energy  rate,  turn  rate  and  Mach  number,  together  with  a  copy  of  the  Mach  number  data  vector. 
This  storage  in  array  WORK  is  carefully  allocated  to  make  use  of  the  “repeat”  cycle  concept18 
used  in  subroutine  LINE  to  plot  multidimensional  data. 

Subroutine  LINE  is  then  called  twice  to  plot  turn  rate  and  energy  rate  curves  for  the  selected 
load  factor.  The  plotter  origin  is  adjusted  after  each  curve,  ana  again  between  each  set  of  graphs 
for  the  selected  height.  Control  returns  to  the  main  program  to  terminate  execution  after  all 
heights  have  been  considered. 


6.  PROGRAM  P4  DESCRIPTION 
6.1  Program  Summary 

Program  P4  is  a  multi-purpose  contour  plotting  program:  input  consists  of  data  files  pro¬ 
duced  by  program  P2.  The  output  is  a  single  file,  P4.PLT,  to  be  submitted  to  the  system  program 
PLOTQ  for  off-line  plotting. 

The  main  program  is  very  short,  and  simply  accepts  an  input  data  filename,  reads  the  grid 
definitions  and  header  text  and  calls  the  major  control  subroutine  P4MAIN.  This  routine  con¬ 
trols  all  logic  functions  for  the  various  type  of  data,  processes  scale  and  contour  level  requests, 
draws  axes  and  controls  pen  positioning  for  each  page  of  plotter  output. 

Subroutine  CONT,  called  from  subroutine  P4MA1N,  processes  the  data  grid  and  searches 
along  the  requested  contours,  plotting  points  wht-re  the  linear  approximation  changes  direction, 
with  calls  to  the  auxiliary  subroutine  P. 

The  following  sections  discuss  the  types  of  plots  which  may  be  produced  and  describe  the 
routines  involved  with  the  control  logic,  subroutines  P4MA1N  and  P.  Subroutine  CONT  and 
several  utility  routines  are  portions  of  a  larger  contouring  package  used  at  ARL. 

Brief  descriptions  of  all  routines  used  by  program  P4  are  given  in  Appendices  8,  9  and  10. 
Examples  of  terminal  input  are  given  in  the  user’s  guide  for  program  P4  in  Chapter  9,  and 
examples  of  contour  plots  produced  by  co-ordinated  running  of  programs  AIRCRAFT,  P2 
and  P4  are  given  in  Chapter  10. 


6.2  Classification  of  Contour  Plots 

The  control  logic  of  subroutine  P4MAIN  is  capable  of  differentiating  four  types  of  contour 
plot,  based  on  input  data  filenames,  as  shown  in  Table  6. 

Each  data  file  contains  grid  definitions  and  text  descriptions  of  the  aircraft  configurations. 
The  energy  rate  data  files  contain,  in  addition,  data  for  plotting  energy  state  contours  if  desired. 
One  page  of  plots  is  produced  for  each  load  factor,  which  may  be  selected  from  the  available 
load  factor  grid. 

The  differential  MMD  data  files  contain,  in  addition  to  the  differential  data,  turn  rate  data 
for  both  aircraft  from  which  zero  turn  rate  (Ig  load  factor)  boundaries  are  obtained.  Contours 
of  differential  turn  rate  are  then  plotted  only  when  they  are  bounded  by  both  of  these  boundaries. 

For  all  types  of  contours,  a  terminal  dialogue  is  used  to  define  the  range  of  contours  required 
and  the  scales  to  be  used  on  the  x  and  y  axes. 
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TABLE  6 


Contour  Plot  Classification 


Plot  |  input  data 
type  i  filenames 

Brief 

description 

1  1  P2.CON, 

i  P2A.CON 

! 

I 

Energy  rate  contour  plots. 

Input  data  contains  energy  state  contour  data,  followed  by  energy 
rate  contour  data  for  selected  load  fuctors. 

2  P2DIFF.CON 

! 

i 

Differential  energy  rate  contour  plots. 

Input  data  contains  descriptive  texts  for  both  aircraft,  followed  by 
energy  state  and  differential  energy  rate  contour  data  as  above. 

3  |  P2.0PT, 
j  P2A.OPT 

Maximum  manoeuvre  diagrams  (MMD). 

Input  data  contains  only  turn  rate  contour  data. 

4  j  P2DIFF.OPT 

1 

i  ! 

i  1 

Differential  MMD. 

Input  data  contains  copies  of  files  P2.0PT  and  P2A.OPT  for  base 
and  comparison  aircraft,  respectively,  followed  by  differential  turn 
rate  contour  data. 

6.3  Subroutine  P4MAIN 

Subroutine  P4MAIN  controls  all  of  the  logic  functions  for  the  four  types  of  contour  plot, 
and  performs  a  dialogue  with  the  user  to  supply  scale  and  contour  level  data.  It  also  controls 
reading  of  the  input  tile,  plots  axes  on  the  output  file,  positions  the  plotter  pen  for  each  page  of 
plots,  and  calls  subroutine  CONT  to  trace  the  requested  contours  through  the  grid  data.  A 
flowchart  for  the  subroutine  is  given  in  Figure  17. 

Several  important  switches  are  used  to  control  the  logic  operations.  The  first,  IMMD  is 
set  to  I  (otherwise  0)  for  MMD  plots,  based  on  input  data  filename.  If  the  MMD  plot  is  also  a 
differential  plot,  the  data  set  loop  counter  NDATA  is  set  to  3  (otherwise  1).  For  the  first  two 
passes  through  the  data  set  loop  (1DATA  =  1  or  2),  switch  IGRID  is  set  to  1  (otherwise  0)  to 
indicate  that  only  the  zero  turn  rate  boundary  is  to  be  plotted.  For  each  of  these  boundaries 
counter  1C  records  the  number  of  points  used  to  define  the  boundary  (limit  of  400  allowed),  and 
switch  LABFLG  is  set  to  .F.  (otherwise  .T.)  to  indicate  that  contours  are  not  be  labelled. 

For  contour  plots,  switch  IGES  is  set  to  1  (otherse  0)  if  energy  state  contours  are  to  be 
plotted;  these  are  plotted  with  a  dotted  line,  and  do  not  require  contour  labels,  being  readily 
identified  by  their  intersection  with  the  height  axis. 

The  subroutine  proceeds  by  plotting  identifying  labels,  accepting  axis  scales  and  positioning 
the  plotter  pen  for  each  set  of  contour  data  on  the  input  file. 

Two  loops  are  set  up  in  P4MAIN.  The  first  of  these  is  for  the  number  of  data  sets  on  the 
data  file,  NDATA.  After  checking  for  a  differential  MMD  label,  the  first  set  of  contour  data  is 
read  into  array  Z.  Pointers  IGRID  and  1ES  are  given  default  values  of  zero,  and  if  a  differential 
MMD  zero  turn  rate  boundary  is  being  plotted,  the  first  word  of  the  contour  request  vector  is 
set  to  zero. 

MMD  plots  then  proceed  immediately  to  the  load  factor  loop  (executed  only  once);  other¬ 
wise  the  energy  state  contour  request  is  processed.  If  these  contours  are  requested,  matrix  Z  is 
copied  into  matrix  ZG,  switch  IGES  is  set  and  the  energy  state  contour  vector  constructed,  using 
the  user’s  terminal  input. 

The  load  factor  loop  is  then  entered  and,  if  energy  rate  contours  are  being  plotted,  data  is 
read  into  matrix  Z  and  the  load  factor  request  processed.  If  the  data  is  to  be  skipped,  control 
jumps  to  the  end  of  the  load  factor  loop.  Energy  rate/MMD  contours  are  then  specified  and  the 
plotter  origin  reset,  x  and  y  axes  are  drawn  (these  are  skipped  if  the  MMD  zero  turn  rate  boundary 
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Fig.  17(a)  Subroutine P4M AIN  Flowchart 
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Fig.  17(c)  Subroutine  P4MAIN  Flowchart  (Cent.) 
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Next  data  set 


RETURN 


9000  j 
^VPe  end  of  data'I 


Fig.  1 7{d)  Subroutine  P4MAIN  Flowchart  (Cont.) 


is  being  plotted),  and  either  a  load  factor  label  or  a  zero  energy  rate  axis  is  drawn,  depending  on 
the  type  of  plot. 

Data  scales  are  set  to  unity  for  plotting  the  contours  themselves,  as  the  contouring  routines 
plot  data  in  real  inches  of  plot.  Scaling  ;0  these  units  is  performed  in  subroutine  P. 

If  energy  state  contours  have  been  requested  on  energy  rate  plots,  these  contours  are  plotted 
first  as  a  broken  >ine  with  a  call  to  CONT,  after  setting  parameters  as  follows: 

Parameter  Value 

CONLAB . T.  to  label  contours,  .F.  to  omit, 

1DOT,  IDEF. .  . .  1  to  draw  contours  with  a  broken  line  (dashes), 

0  fo'  a  solid  line,  and  --3  for  an  interrupted 
line  (dot-dash). 

NLEV0  . .  Level  number  of  zero  contour 

A  second  call  to  CONT  plots  energy  rate/MMD  contours  with  a  solid  line,  except  for  the  zero 
contour  which  is  plotted  using  a  dot-dash  format  to  distinguish  it  from  other  contours.  Sharp 
corners  of  the  contours  are  smoothed  using  a  smoothing  interval  of  0  05  in.;  calls  to  subroutine 
SMOOTH  with  a  pen  command  of  4  are  required  whenever  the  type  of  contour  line  changes,  to 
finish  the  current  line. 

For  the  remainder  of  the  load  factor  loop,  if  MMD  boundaries  are  being  plotted,  scales  are 
reset  to  unity  and  the  pen  moved  to  the  origin  before  jumping  to  the  end  of  the  loop.  Otherw  ise, 
the  list  of  contour  levels  is  plotted  before  resetting  plotting  scales  to  unity  and  positioning  the 
pen  for  the  next  plot. 

When  all  load  factors  and  data  sets  on  the  data  file  have  been  considered,  control  returns 
to  the  main  program  to  terminate  execution. 

Errors  when  replying  to  terminal  prompts  result  in  the  prompts  being  repeated.  Encountering 
an  eud-of-file  on  the  input  file  before  the  end  of  the  data  set  or  load  factor  loops  results  in  an 
“end  of  file”  text  being  typed  on  the  terminal  before  ceasing  execution. 


6.4  Subroutine  P 

Subroutine  P  is  an  interface  routine  between  subroutine  CONT  and  the  plotter  software, 
which  renders  CONT  independent  of  plotter  conventions.  It  modifies  pen-up  and  pen-down 
movements  commanded  by  CONT,  taking  into  account  the  zero  turn  rate  boundary  when  pro¬ 
ducing  differential  MMD  plots.  It  also  controls  the  parameter  IDOT  specifying  the  type  of  line 
used  in  plotting  a  contour.  A  flowchart  for  the  subroutine  is  given  in  Figure  18. 

P  is  called  with  co-ordinate  and  pen  command  arguments  supplied  by  CONT.  (XPT,  YPT) 
are  the  co-ordinates  of  the  point  to  be  plotted  in  mesh  co-ordinates,  which  have  ranges  (1,  NXM) 
and  (1,  NES)  for  XPT  and  YPT  respectively.  The  co-ordinates  (x,  y)  are  adjusted  to  a  zero  origin, 
and  a  local  copy  of  the  pen  command  is  made  in  NSAV. 

A  test  is  next  made  of  IDATA,  the  counter  for  the  loop  in  P4MAIN  which  has  called  CONT. 
If  a  differential  MMD  boundary  is  being  plotted,  P  also  stores  copies  of  the  co-ordinates  in 
vectors  XTAB  and  YTAB  for  future  reference.  Counter  IC  is  incremented  each  time  a  point  is 
stored,  and  the  two-word  vector  1COUNT  records  the  value  of  IC  at  which  data  for  the  two 
boundaries  finishes.  If  the  maximum  values  of  IC  of  400  is  exceeded,  execution  ceases  after 
typing  a  failure  text  on  the  terminal.  The  minimum  value  of  y  on  either  boundary  is  recorded  in 
YLMIN,  and  flag  1CHK  set  to  zero.  This  records  that  the  order  cf  x  values  in  XTAB  is  undefined 
(ascending  or  descending).  A  jump  is  then  made  to  statement  1000  to  plot  the  boundary  point. 

For  all  other  types  of  contour  plot,  y  is  first  tested  against  YLMIN.  Unless  a  differential 
MMD  is  being  plotted,  this  has  a  default  value  equal  to  the  maximum  y  mesh  co-ordinate;  hence 
the  test  is  always  satisfied  and  a  jump  always  made  to  statement  1000  to  plot  the  point.  If  the 
test  is  not  satisfied,  the  point  (x,jO  must  be  further  tested  against  both  of  the  stored  MMD 
boundaries. 

The  first  time  this  area  of  code  is  entered,  ICHK  is  zero,  and  the  data  in  vectors  XTAB 
and  YTAB  are  rearranged  if  necessary  to  guarantee  each  boundary  is  stored  in  ascending  order 
of  x  values.  This  is  done  in  the  loop  ending  at  statement  580.  ICHK  is  set  to  1  so  that  future  calls 
to  P  skip  this  rearrangement. 
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CALL  P  (XPT,  YPT,  N) 


Um  local  pan  command: 
HSAV-N 

Convert  to  zero  origin: 

X  ■  XPT  —  1 Y-  YPT  -  1. 


Type  'array 
size  exceeded' 


<IDATA>tNDATA?< 
(Differential  MMD, 
’'^•boundary?)1^^ 


5oqjL^ 

^.Y<YLMIN?Vw 
<(Y  below  min.  bounder 


,^ICL<400?>> 
Boundary  storage 
available?**''* 


Count  b'dry  points:  1C  ■  1C  +  1 
!  Store:  ICOUNT  (IDATA)  -  1C 
XTAB  (IC>  -  X;  YTAB  (1C)  -  Y; 

Set  min.  Y  on  boundaries: 
YLMIN  ■  AMIN1  (Y.YLMIN); 
Order  flag:  ICHK-0 


Ensure  boundary  tables 
are  in  the  order  of 
increasing  X. 


-^Tchk-  ir^ 

(order  checked)? 


Set  flags: 
ICHK-  1;  IL  -  1 


Boundary  Ipop 
DO  680  J-  1.2 


Set  table  pointers 
I  LI  -  IL;  IL2  -  IL  +  1 


Fig.  18(a)  Subroutine  P  Flowchart 
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<  XTAB  (KR-S^ 
X  Co-ord.  bounded^ 


Interpolate  in  XTAB 
YTAB  for  YUM 


Y  <  YUM?  . 
(X,  Y)  under 
^boundary?-'* 


Lift  Pen: 
NSAV  -  5 


© 


1 
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Fig.  18(d)  Subroutine  P  Flowchart  (Cont.) 
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Subsequent  differential  MMD  calls  jump  to  statement  599  to  compare  the  .v  co-ordinate 
with  the  YTAB  co-ordinates  at  the  current  value  of  x.  The  XTAB  values  for  each  boundary 
ere  bounded,  and  interpolation  yields  two  values  against  which  y  must  be  compared.  If  y  is 
greater  than  either  of  these  two  values,  the  pen  is  lifted  (NSAV  ™  5)  before  proceeding. 

At  statement  1000,  P  is  now  ready  to  plot  the  point  and  the  co-ordinates  are  converted  to 
inches  of  plot,  if  contour  smoothing  is  not  requested,  the  point  is  plotted  by  a  call  to  PLOTD. 
Otherwise  it  is  plotted  with  a  call  to  SMOOTH.  In  the  latter  case,  the  texture  of  xero  contours 
is  varied  by  changing  the  parameter  1DOT.  However,  each  time  the  texture  is  changed,  an 
additional  call  to  SMOOTH  is  required  to  terminate  the  previous  contour.  Flag  1DOTC  is  used 
to  indicate  a  change  in  contour  texture. 

The  beginning  of  contours  are  flagged  with  their  contour  level  number  via  a  call  to  routine 
OUTXT,  unless  the  pen  has  been  lifted  by  P  or  energy  state  contours  are  being  plotted  (CONLAB 
is  then  .FALSE.).  Having  plotted  the  point  control  returns  to  the  calling  routine,  CONT. 


6.5  Communication  with  System  Routine  PLOT 

The  basic  unit  for  plotting  a  point  using  the  CALCOMP  plotting  software16  is  a  call  to 
PLOT  of  the  form 


CALL  PLOT  (LU,  X,  Y,  N), 

where  LU  is  the  logical  unit  number  of  the  output  channel; 

(X,  Y)  are  the  co-ordinates  of  the  point ; 

N  is  an  integer  specifying  the  pen  command — 

1  use  (X,  Y)  as  axis  scaling  factors, 

2  establish  new  origin  with  current  pen  position  at  (X,  Y), 

3  lift  pen,  move  to  (X,  Y),  drop  pen, 

4  drop  pen,  draw  straight  line  to  (X,  Y), 

5  lift  pen,  move  to  (X,  Y). 

The  pen  commands  used  by  the  contouring  routines  are  the  same  as  these,  except  that  the 
meaning  for  N  =  4  has  been  transformed  to  N  =  2.  Hence  the  sole  task  of  routine  PLOTD, 
whose  call  statement  is  identical  to  that  above,  is  to  replace  a  value  of  N  =  4  by  N  =  2,  before 
calling  PLOT. 

Routine  SMOOTH  also  plots  by  calls  to  routine  PLOTD,  but  the  call  includes  the  additional 
calling  parameter  1DOT,  used  to  indicate  the  mark  :  space  ratio  for  broken  lines.  A  smoothing 
interval  of  0  05  is  transmitted  via  the  COMMON  area  GDMDOT. 

In  operation,  SMOOTH  has  a  plotting  lag  of  one  point,  required  by  the  splining  technique 
used.  Successive  points  are  joined  by  cubic  arcs  with  tangential  coincidence  at  the  end  points. 
The  plotting  lag  results  in  an  additional  call  to  SMOOTH  being  required,  to  finish  the  previous 
contour,  whenever  a  change  in  contour  texture  is  made.  The  pen  command  N  =  4  is  used  to 
request  this  updating  procedure. 


7.  PROGRAM  “AIRCRAFT”  USER’S  GUIDE 
7.1  Loading  and  Saving  of  Absolute  Files 

Given  a  file  A1RCRA.FT  containing  a  main  program,  the  BLOCK  DATA  subroutine  and 
thrust  and  drag  routines,  the  commands  for  loading  the  file  and  producing  an  absolute  copy  of 
the  program  on  disk,  together  with  a  symbol  map,  are : 
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.LJAP/MAP  AIRCRA. H./5EARCH  PI. P !UF 

FORTRAN:  AIRCRA 

MAIN. 

.BLOCK 

TRINCL 

AERO 

THRUST 

LINK:  Loading 

EXIT 

.SAVE 

AIRCRA  saved 


In  this  example,  and  in  subsequent  examples  of  terminal  operations,  user  responses  are 
underlined. 

The  above  operations  produce  a  file  AIRCRA.EXE  on  disk  which  is  then  executed  with  the 
RUN  command. 


7.2  Program  Execution 

"Tie  conversational  input  to  the  program  is  best  described  with  reference  to  a  typical  run  to 
produce  data  for  plots  of  energy  rate  against  turn  rate  at  several  altitude  /Mach  number  com¬ 
binations.  Numbers  in  the  left-hand  margin  indicate  subsections  in  which  the  various  responses 
are  discussed. 


7.2.1  Imperial  (0)  or  SI  (1)  Units 

Replies  of  0  or  1  indicate  that  output  quantities  are  to  be  in  Imperial  or  SI  units  respectively. 
Any  other  reply  results  in  the  text 

INVALID  INPUT!  TRY  AGAIN 

being  typed  on  the  user's  terminal,  and  the  prompt  is  repeated. 
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.RUN  A1RCRA 


MANEUVERABILITY  ARID  CALCULATION 

1  IMPERIAL  (0)  OR  8.1.  (1)  UNITS  T  1 

2  AIRCRAFT  DATA  FILENAME  T  ASTORE 

AIRCRAFT  EXAMPLE  BATE  1 -Apr-82  TIME  12x17  37.1 
AIR  SUPERIORITY  -  2  AAM  ♦  SUN  +  502  FUEL 

3  PS,  P8/UF  OR  PS*V/UF  (1,  2  OR  3)  T  J. 

4  IS  HP  THE  HEIGHT  VARIABLE  ?  X 

5  PRESET  GRID  (Y,  N  OR  C/R)  ?_ 

PRESET  GRID  : 

hP  =  0(4000)68000  (IMP.)  OR  0(1000)20000  (S.I.) 

(6)  MN  ■  0(0.05)1 .5  (MIL)  OR  0(0.05)2.0  (MAX) 

G  =  1(2)9 

PRESET  GRID  (Y,  N  OR  C/R)  ?  N 

7  HPO  ,  HP  STEP  (M  ),  NO.  OF  STEPS  ?  0,5000,4 

8  HACHO,  MACH  STEP  ,  NO.  OF  STEPS  ?  .8,.V 

INVALID  INPUT!  TRY  AGAIN  - 

MACHO,  MACH  STEP  ,  NO.  OF  STEPS  ?  ,8, .1,2 

9  GO  ,  G  STEP  ,  NO  OF  STEPS  ?  1 ,.25,24 

10  POUER  <HIL=100,MAX=200)  ?  200 

11  UING  SUEEP  (IF  VARIABLE)  ? _ 

12  ATH05PHERE  ,  DEVIATION  ? _ 

ATHOSPHERE  IS  INVALID;  ICAO  ASSUMED 

13  OUTPUT  (1*TEXT, 2-NOS. ,3*B0TH>  ?  3. 

14  0/P  FILENAME  FOR  UNIT  6  ?  TEST.LST 

0/P  FILENAME  FOR  UNIT  8  ?  TEST.NUH 

CALCULATION  : 

ALTITUDE  HP  *  0.0  M 

ALTITUDE  HP  =  5000.0  M 

ALTITUDE  HP  a  10000.0  M 

ALTITUDE  HP  *  15000.0  M 

STOP 

END  OF  EXECUTION 

CPU  TIME:  6.96  ELAPSED  TIME:  3:13.76 
EXIT 
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7.2.2  Aircraft  Data  File  Name 

Supply  up  to  10  characters  giving  the  name  of  a  file  or.  disk  containing  data  for  the  given 
aircraft  for  the  current  run.  In  the  example,  file  ASTORE  is  as  follows: 


EXAHPLE 

350. . 0 

20442. . 30. 

AIR  SUPERIORITY  -  2  AAM  +  SUN  +  50Z  FUEL 
CF.BIN 
CA.BIN 
8 

.2,. 92, .96, .98, 1.01, I. 1,1. 2, 2.0 
8,8,10,13,17,15,14,14 


Items  in  this  file  are  as  follows 


Line 

Format 

Description 

1 

2A5 

Aircraft  name — up  to  10  characters 

G 

Wing  reference  area  (ft2  or  m2),  file  units  flag  (0  or  1  respectively) 

2G 

Gross  weight  (lb  or  kg),  c.g.  position  (%MAC) 

14A5 

Role  description — up  to  70  characters 

2A5 

Thrust  data  filename — up  to  10  characters 

2A5 

Aerodynamic  data  filename — up  to  10  characters 

G 

No.  of  points  in  store  drag  table 

8 

10G 

Mach  no.  list  in  store  drag  table  L 

9 

10G 

Drag  count  list  in  store  drag  table  j 

*  10  items  per  line;  carry  on  to  next  line  if  more  than  10  entries  on  any  line.  Not 
needed  if  0  points  indicated  by  item  7. 


As  shown,  reference  area  and  weight  are  in  units  indicated  by  the  file  units  flag;  0  indicates 
Imperial  units  and  1  indicate  SI  units. 

Action  in  the  case  of  input  error: 

(a)  If  the  specified  file  does  not  exist,  an  erro  message  is  typed  and  a  new  filename  is  requested. 
Either  type  the  correct  filename  or  return  to  the  monitor  (fC)  and  generate  the  required 
file. 

( b )  If  items  in  the  file  are  incorrect  or  out  of  order,  no  checking  is  possible  and  execution 
errors  will  result. 
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7.2.3  PS,  PS/WF  or  PS.V/WF  (1, 2  or  3) 

Replies  of  1,  2  or  3  indicate  as  follows: 

Reply  Meaning 

1  Calculate  P&  as  dependent  variable.  This  is  the 

most  common  case.  Use  for  minimum  time 
climb  schedule  (1;  load  factor). 

2  ..  ..  Calculate  Pt/wt  as  dependent  variable.  Use 

for  minimum  fuel  climb  schedule  (1;  load 
factor). 

3  Calculate  PaV/lOOOwi  as  dependent  variable. 

Use  for  maximum  range  climb  schedule  (Ig 
load  factor). 

Values  1,  2  and  3  are  valid  as  input.  Any  other  value  produces  the  text 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed,  and  the  prompt  is  repeated. 

The  type  of  energy  parameter  indicated  by  IPSTYP  and  units  indicated  by  IUNITS  are 
evident  from  the  printed  listing;  this  information  is  also  transmitted  to  input  files  to  programs 
P2  and  P4  as  well,  so  that  all  outputs  from  all  programs  carry  correct  identification  of  output 
variables  and  units. 

7.2.4  Is  HP  the  Height  Variable? 

A  reply  “Y”  (for  “YES”)  indicates  that  pressure  altitude  is  to  be  used  as  the  height  grid 
variable.  Any  other  reply  will  result  in  energy  state  being  used  as  the  height  variable. 

7.2.5  Preset  Grid  (Y,  N  or  C/R) 

For  some  contour  plots,  such  as  differential  energy  rate  plots,  a  preset  or  standardised  grid 
may  be  useful.  A  carriage  return  reply  results  in  the  preset  grid  values  being  typed  on  the  uset’s 
terminal,  as  shown.  The  values  indicated  for  each  grid  are  in  the  form  “initial  value  (increment) 
final  value”.  A  reply  “Y”  results  in  this  grid  being  used  for  the  pressure  altitude  (no  preset  grid 
is  available  for  energy  state);  grid  requests  7,  8  and  9  are  then  omitted. 

Any  other  reply  will  result  in  requests  for  grid  data  for  pressure  altitude,  Mach  number  and 
load  factor. 


7.2.6  Maximum  Manoeuvre 

This  prompt  (not  shown  in  the  example  given)  is  made  if  energy  state  is  selected  as  the  height 
variable.  Any  reply  other  than  “Y”  will  result  in  an  unoptimised  grid,  and  a  Mach  number  grid 
will  be  requested. 

Reply  “Y”  will  cause  an  optimisation  of  energy  rate  as  a  function  of  Mach  number,  prepar¬ 
ing  data  for  a  maximum  manoeuvre  diagram  (MMD),  climb  schedules  or  optimum  sustained 
turn  rates. 


7.2.7  HP0,  HP  Step,  Number  of  Steps 

This  prompt  requests  parameters  defining  the  height  grid  (“HP”  is  replaced  by  “ES”  in 
the  prompt  if  energy  state  is  the  height  variable).  Three  items  are  requested,  separated  by  commas. 
These  are 

•  initial  value, 

•  a  positive  increment, 

•  number  of  points,  including  first  and  last. 

Units  for  initial  value  and  increment  (FT  or  M)  are  indicated  in  the  prompt. 
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A  valid  combination  of  the  three  parameters  is: 

(!)  initial  value  >  0, 

(ii)  increment  value  >  0, 

(iii)  0  <  number  of  values  <  SO. 

Any  other  combination  causes  the  text 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed,  and  the  prompt  is  repeated. 

7.2.8  MACHO,  MACH  Step,  Number  of  Steps 

The  Mach  number  grid  request  is  similar  to  the  height  grid  request  together  with  its  error 
test.  This  prompt  is  not  given  if  an  optimal  grid  has  been  requested. 

7.2.9  GO,  G  Step,  Number  of  Steps 

The  load  factor  grid  request  is  also  similar  to  the  height  grid  request. 

7.2.10  Power  (MIL  =  100,  MAX  =  200) 

Military  power  is  indicated  by  a  value  in  the  range  0  to  100,  and  maximum  power  in  the 
range  100  to  200.  Tens  and  units  digits  indicate  a  percentage  of  the  maximum  of  that  range. 

A  valid  reply  is  0  <  Power  <  200.  Any  other  reply  causes 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed,  and  the  prompt  is  repeated. 


7.2.11  Wing  Sweep  (if  Variable) 

Supply  wing  sweep  in  degrees  if  aerodynamic  data  requires  it.  Otherwise  reply  with  a  carriage 
return. 

A  nominal  valid  reply  is  0  <  Wing  sweep  <  80.  Any  other  reply  causes 
INVALID  INPUT!  TRY  AGAIN 
to  be  typed,  and  the  prompt  repeated. 

If  any  other  wing  sweep  limit  is  desired,  it  should  be  included  in  the  aerodynamic  routines 
for  the  aircraft. 


7.2.12  Atmosphere,  Deviation 

This  is  an  input  provided  if  a  THRUST  subroutine  is  loaded  for  a  given  aircraft  which  can 
handle  alternative  atmospheres.  Valid  replies  are 

ARDuS  )F0RMAT  <A4'°> 

where  X  is  a  deviation  from  the  nominated  atmosphere  in  degrees  Celsius.  It  may  be  omitted  if 
X  is  zero. 

Action  in  Case  of  Input  Errors 

If  the  first  four  characters  are  neither  ICAO  nor  ARDU,  the  text 

ATMOSPHERE  (name)  IS  INVALID;  ICAO  ASSUMED 
is  typed,  and  calculation  proceeds  assuming  an  ICAO  atmosphere. 
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7.1.13  Output  (1  m  Tixt,  2  *1  N«k,  3  *»  Both) 

Replies  1, 2  or  3  indices  u  follows: 

(1)  formatted  output  with  fell  ASCII  text  is  requested; 

(2)  numerical  (data)  output  is  required,  for  input  to  program  P2; 

(3)  both  types  of  output  are  required. 

A  reply  K  it  valid  if  1  <  K  <  3.  Any  other  reply  causes 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed,  and  the  prompt  is  repeated. 

7.2.14  O/P  Name  for  Unit  re 

Supply  a  \0- character  filename  for  output  as  requested,  m  -  6  is  the  logical  unit  number  for 
foil  format  output,  m  —  8  is  the  logical  unit  number  for  numerical  output. 

Formatted  output  should  be  printed  with  /P/B  switches  to  produce  listings  suitable  for 
permanent  retention. 

Numerical  output  may  be  printed  if  required. 

7.3  Core  Storage  Requirement*  and  Execution  Speed 
Storage  requirements  and  execution  speed  are  both  very  much  dependent  on  user-defined 
routines  for  (he  propulsion  and  aerodynamic  cakrlations. 

A  typical  example  requires  a  total  coresire  of  S3  pages  (about  2/Kg  words)  of  core  ard  pro¬ 
cesses  an  energy  state  contour  plot  request  using  the  default  grid  in  3  min.  18  sec.  of  CPU  time. 
This  example  uses  a  thermodynamic  model  of  the  propulsion  system  giving  thrust  and  fuel 
flow  data  at  ali  power  settings,  and  B-spline  representations  of  maximum  Ci„  Ct  versus  a  and 
full  drag  polars  for  one  representative  c.g.  location. 

Progress  of  five  calculation  is  indicated  on  the  user’s  terminal  by  typing  out  each  value  of  the 
height  grid  as  it  is  processed. 

8.  PROGRAM  P2  USER’S  GUIDE 
8.1  Loading  and  Saving  of  Abaolute  Files 

The  commands  for  loading  program  P2  and  saving  an  absolute  copy  on  disk,  together  with 
a  symbol  map  are: 

.LQA/NAP  P2,P?4LIB/SEARCH 

FORTRAN:  P2 

MAIN. 

GRID 

INLAB 

INMMD 

NHP 

PLTLAB 

PSCON 

P8DIFF 

P2IN 

RATE! 

RATE2 

R2DIFF 

UNITS 

URLAB 

LINK:  Loading 

EXIT 

•ML 

P2  savtd 
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The  absolute  copy  of  the  program  is  saved  on  the  disk  file  P2.EXE,  which  is  then  executed 
with  the  RUN  command.  Core  storage  requirement  is  34  pages  (17Kg  words).  If  the  storage 
allocation  for  array  WORK  (8000  words)  is  found  to  be  insufficient  for  large  grids,  there  is  thus 
ample  core  storage  in  reserve,  should  a  program  change  be  necessary. 

8.2  Program  Execution 

The  basic  inputs  to  program  P2  consist  of  data  files  created  in  a  standard  format  by  program 
AIRCRAFT.  The  data  are  then  manipulated  according  to  commands  supplied  at  the  user’s 
terminal  during  the  execution  of  the  program.  Some  options  allow  the  user  to  supply  on-line 
data,  in  which  case  the  data  format  is  specified  during  the  terminal  dialogue.  In  either  case, 
run-time  input  file  preparation  is  not  required  for  program  P2,  all  data  files  being  generated  by 
program  AIRCRAFT. 

Output  files  are  of  two  types.  Options  4B  and  4Z  produce  files  for  off-line  plotting  using 
the  system  program  PLOTQ.  Other  options  produce  data  files  for  input  to  the  contour  plotting 
program  P4. 

The  remaining  sections  of  this  chapter  describe  the  terminal  operations  required  for  each 
option,  with  reference  to  example  dialogue.  User  input  to  terminal  prompt  is  underlined. 


8.3  Option  4A — Energy  Rate  Contour  Data 


.RUN  P2 


COMBAT  PERFORMANCE  PROCESSING 


OPTION  OR  (CR)  FOR  KELP  :  4A 

4A  PS  CONTOUR  PLOT.  (P2.C0N) 

DATA  BASE  FILENAHE  :  A.NUH 

***  DATA  ARE  IN  METRIC  UNITS,  ENERGY  PARAMETER  IS  PS/UF  *t* 

REPLY  "YES",  "NO”,  "ALL"  OR  "END"  : 

OUTPUT  DATA  FOR  N  *  1.00  ?  1 

OUTPUT  DATA  FOR  N  *  3.00  ? _Y 

OUTPUT  DATA  FOR  N  =  5,00  ?  _Y 

OUTPUT  DATA  FOR  N  »  7.00  ?  N. 

OUTPUT  DATA  FOR  N  =  9.00  ? 

END  OF  EXECUTION 

CPU  TIME!  0.88  ELAPSED  TIME:  1 : t 1 .50 
EXIT 


This  option  produces  energy  rate  contour  data  on  a  Mach  number  (x  axis)  versus  altitude 
(y  axis)  grid,  at  selected  values  of  load  factor. 


63 


i'1  "i'»»«'>  v^tTjvTvn-iT** 


l'*'T7fJ"V?J^>f,Vl(T*’,T’ 


’"  T’-^-V  T  -n»w 


84.1  Option,  or  (CR)  for  Kelp 

Reply  with  the  appropriate  option  code.  A  carriage  return  (CR)  is  interpreted  as  a  request 
for  help,  and  a  one  line  description  of  each  option,  together  with  the  output  files  produced,  is 
typed  on  the  terminal.  An  example  of  this  help  text  is  given  in  Section  S.2. 

84.2  Data  Bate  Filename 

Reply  with  the  name  of  the  appropriate  data  file  produced  by  program  AIRCR  AFT  as 
logical  unit  8.  After  trading  the  header  on  the  files,  P2  echoes  the  unit  system  and  type  of  energy 
parameter  on  the  user’s  terminal. 

MJ  Ontpot  Data  for  N  —  bums 

Reply  with  “Y”  or  "N",  depending  on  whether  a  contour  plot  is  required  for  the  load 
factor  specified  or  not.  Output  is  written  on  file  P2.CON,  for  input  to  program  P4.  This  file  is  in 
ASCII  and  may  be  printed,  but  it  is  not  formatted  for  output  listing. 

8.3.4  Input  Errors 

(a)  If  the  specified  input  file  does  not  exist,  an  error  message  is  typed  and  a  new  filename  is 
requested.  Supply  the  correct  name  to  continue  or  (fC)  to  abort  the  job. 

(b)  Replies  other  than  Y  to  load  factor  requests  are  interpreted  as  N. 

8.4  Option  4B — Turn  Rate  Plots 

OPTION  OR  < CR )  FOR  HELP  :  4B 

4B  PS  VS  TURN  RATE  FOR  GIVEN  HEIGHT.  (P2.PLT) 

DATA  BASE  FILENAHE  :  B.NUH 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

SCALES  IN  UNITS/INCH  OF  PLOT 
QHEGA  /  ENERGY  VARIABLE  ;  4,400 

MAXIMUM  FOR  ENERGY  AXIS  :  .1200 

ENERGY  AXIS  LENGTH  (IN)  s  2 

REPLY  "YES",  "NO",  "ALL"  OR  "END"  ! 


HEIGHT 

= 

0,0  FT 

? 

x 

MACH 

S 

.800  ?  X 

MACH 

3 

.900  ?  N. 

HEIGHT 

8 

10000.(5  FT 

1 

x 

MACH 

3 

.800  ?  _N 

MACH 

S 

.900  ?  X 

HEIGHT 

Z 

20000.0  FT 

? 

A 

HEIGHT 

3 

30000.0  FT 

T 

E 

END  OF  EXECUTION 

CPU  TINE:  6.80  ELAPSED  TINE:  2:6.50 
EXIT 
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This  option  produces  plots  of  energy  rate  against  turn  rate  for  selected  values  of  Mach 
number.  One  page  of  plots  is  produced  for  each  requested  altitude. 

8.4.1  Data  Base  Filename 

Reply  with  the  appropriate  filename,  as  in  Section  8.3.2. 

8.4.2  Scales  in  Unlts/lnch  of  Plot 

Reply  with  the  amount  of  the  physical  quantity  which  one  inch  of  plot  would  represent  on 
each  axis.  The  scales  presented  in  Table  5  may  be  used  as  a  guide  in  selecting  suitable  values. 

OMEGA,  the  turn  rate  variable,  will  vary  from  0  to  20  deg/s  or  more,  and  4  deg/s  per  inch 
is  a  representative  value.  OMEGA  is  plotted  as  the  x  axis. 

PS  will  usually  represent  specific  excess  power  and  will  have  significance  between  values  of 
(approxima.ely)  —600  and  +400  m/s  (  —  1800  and  +1200ft./s).  A  value  of  150  m/s 
(400  ft/s)  per  inch  is  representative.  When  PS  represents  other  energy  parameters,  the 
scales  of  Table  5  should  be  used  as  a  guide.  PS  is  plotted  as  the  y  axis. 

8.4.3  Maximum  for  PS  Axis 

Reply  with  the  approximate  maximum  PS  value  of  interest.  T  his  value,  together  with  the 
PS  scale,  determines  the  range  of  values  on  the  y  axis. 

8.4.4  PS  Axis  Length  (in.) 

Reply  with  the  length  in  inches  of  the  vertical  axis.  A  value  of  7  in.  is  representative  for  easy 
trimming  to  A4  size. 


8.4.5  Height  =  nnmum.n  M  (ft) 

One  page  of  curves  is  plotted  for  each  height  selected.  The  replies  and  their  significance  are 
indicated  below. 


Reply 
Y(ES)  . . 

N(O)  . . 

A(LL)  . . 

£(ND)  .. 


8.4.6  Mach  =  n.nnn 


Meaning 

Plotting  is  required  for  that  height  value.  Mach  numbers 
will  be  listed  from  which  a  choice  is  to  be  made. 

Plotting  is  not  required  for  that  height  value.  Go  to  the 
next  set  of  data,  for  which  a  further  height  prompt  will 
be  given. 

Plot  all  the  data  given  for  that  height  value.  Use  this  reply 
when  the  Mach  number  set  is  known,  to  avoid  repe¬ 
tition  of  Mach  number  requests. 

Finish  considering  height  values.  This  in  effect  terminates 
program  execution  before  all  the  data  has  been  pro¬ 
cessed. 


The  replies  and  their  meanings  are  indicated  as  follows. 
Reply  Meaning 


Y(ES) 

N(O) 

A(LL) 

E(ND) 


Plot  a  curve  for  this  Mach  number. 

Skip  this  Mach  number. 

Plot  curves  for  all  remaining  Mach  numbers  in  the  data, 
for  the  current  height. 

Skip  the  remaining  Mach  numbers  in  the  data,  and  go  on 
to  the  next  height  value. 
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The  output  file  P2.PLT  contains  the  requested  plots,  and  is  submitted  to  the  off-line  plot  queue 
by  running  the  system  program  PLOTQ.  As  a  rough  guide,  each  height  value  requires  about 
30  cm  (12  in.)  of  plot,  depending  on  the  x  axis  scale  used. 


8.4.7  Input  Errors 

(a)  Incorrect  filename— as  in  Section  8.3.4. 

(b)  If  either  of  the  axis  scales  is  zero,  the  message 

INVALID  INPUT!  TRY  AGAIN 
is  typed,  and  new  scales  are  requested. 

(c)  Any  reply  to  the  HEIGHT  or  MACH  prompts,  other  than  the  four  listed,  is  interpreted 
as“NO”. 

8.5  Option  4C — Maximum  Manoeuvre  Diagram  (MMD) 


OPTION  OR  <CR)  FOR  HELP  :  4C 

4C  MAXIMUM  MANEUVER  DIAGRAM  -  MMD.  (P2.0PT) 

ON-LINE  DATA  ?  N^ 

DATA  BASE  FILENAME  s  C.MMD 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

DEFINE  ENERGY  GRID  -  (POSITIVE  INCREMENT) 

START,  STEP,  NO  OF  STEPS;  -1200,100,25 

END  OF  EXECUTION 

CPU  TIME;  0.64  ELAPSED  TIME;  49.46 
EXIT 


This  option  produces  turn  rate  contour  data  on  an  energy  state  (x  axis)  versus  optimum 
energy  rate  (y  axis)  grid. 


8.5.1  On-Une  Data? 

Reply  “Y”or  "N”  depending  on  whether  data  from  a  prior  MMD  is  to  be  supplied  on-line 
or  not.  If  the  reply  is  "N",  data  will  be  expected  on  a  file  produced  by  piogiam  AIRCRAFT. 
(An  example  of  on-line  data  is  given  below.) 


8.5.2  Data  Base  Filename 

Reply  as  in  Section  8.3.2. 
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8.5.3  Energy  Grid— Start,  Step,  Number  of  Steps 

The  reply  to  this  prompt  defines  ?.n  energy  parameter  grid  which,  together  with  the  energy 
state  grid  defines  the  independent  variable  for  the  MMD. 

Three  systems  are  requested,  separated  by  commas: 

•  minimum  value ; 

•  increment  value— must  be  positive ; 

•  number  of  values,  including  first  and  last. 

Typical  ranges,  together  with  appropriate  units  depending  on  the  values  of  parameters 
IPSTYP  and  1UN1TS,  are  given  in  Table  5. 

8.5.4  Input  Errors 

(a)  Incorrect  filename— as  in  Section  8.3.4. 

\ 

(b)  Energy  rate  grid — a  valid  combination  will  ensure: 

(i)  increment  value  >  0;  and 

(ii)  0  <  number  of  values  <  50. 

An  invalid  combination  causes 

INVALID  INPUT!  TRY  AGAIN 

to  be  typed  on  the  terminal,  and  the  prompt  repeated. 

Output  is  written  in  ASCII  mode  on  file  P2.0PT  for  input  to  program  P4.  Again,  this  file 
may  be  printed  if  desired,  but  it  is  not  formatted  for  output  listing. 

8.5.5  On-line  Data  Example 

A  reply  “Y”  to  “ON-LINE  DATA?"  results  in  a  dialogue  as  shown  over.  Data  defining 
system  of  inputs,  type  of  energy  parameter,  aircraft  name,  role,  weight,  wing  sweep  (optional!, 
power  setting  and  atmosphere  are  supplied  in  response  to  the  appropriate  prompts.  1  he  energy 
state  and  energy  rate  grids,  over  which  turn  rate  data  are  to  be  plotted,  are  also  given. 

8.5.6  nnnnnn.M  (ft) 

For  each  energy  state,  pairs  of  points  on  an  energy  rate  vs.  turn  rate  curve  are  supplied. 
The  turn  rate  values  must  be  in  increasing  order;  a  negative  value  of  turn  rate  indicates  the  end 
of  data  for  that  energy  state.  Linear  interpolation  at  fixed  levels  of  energy  rate  are  performed 
using  these  points,  so  the  accuracy  of  the  final  diagram  depends  largely  on  using  increments  of 
turn  rate  over  which  a  linear  approximation  is  valid.  (In  the  example  given,  a  coarse  grid  is  used 
to  abbreviate  the  text.)  In  practice,  a  grid  interval  of  1500  m  (5000  ft)  in  energy  state,  together 
with  about  ten  pairs  of  points,  will  produce  a  reasonably  smooth  contour  plot. 

8.5.7  On-line  Input  Errors 

(a)  Energy  state  and  energy  rate  grids — valid  combinations  of  parameters  are  such  that 

(i)  increment  value  >  0;  and 

(ii)  0  <  number  of  values  ^  50. 

An  invalid  combination  causes 

INVALID  INPUT!  TRY  AG.  ’N 
to  be  typed  on  the  terminal,  and  the  prompt  repeated. 

( b )  Non-increasing  values  of  turn  rate  cause 

? OMEGA  NOT  INCREASING— START  AGAIN: 
to  be  typed  on  the  terminal,  and  the  loop  for  accepting  data  points  to  be  restarted. 

(c)  Typing  mistakes  in  values  of  turn  rate  or  energy  rate  may  be  corrected  after  running 
program  P2  by  using  the  text  editing  program  TECO. 
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run  P2 


COHi AT  PERFORMANCE  PROCESSING 
OPTION  OR  <CR>  FOR  HELP  t  4£ 

4C  MAXIMUM  MANEUVER  DIAGRAM  •  HMD.  (F2.0PT) 

ON-LINE  DATA  T  V_ 

IMPERIAL  (0)  OR  S.I.  <U  UNITS  *  ± 

DATA  REPRESENTS  PS, 

PS/UF  OR  PS*V/ 1 OOOUF  ((,2,  OR  3)  » 

AIRCRAFT  NAME  t  EXAMP 

ROLE  :  AIR  SUPERIORITY 

WEIGHT  <  NO),  UI  MG  SWEEP  (DEO  t  9000 

ESO,  ES  STEP  (  N  ),  NO.  OF  STEPS  t  3000,3000,5 

PSO,  P3  STEP,  NO.  OF  STEPS  x  -300,25,25 

PCUER  <HIl«100,HAX>200)  l  200 

ATMOSPHERE,  DEVIATION  :  ICAO 

SUPPLY  (OMEGA.PS )  PAIRS,  OMEGA  INCREASING. 

FINISH  EACH  ES  SET  U1TH  NEGATIVE  OMEGA. 

3000.  M 

1  t  0,138 

2  :  4,133 

3  «  7.5,110 

4  x  14.4,0 

5  l  19.7,-270 

4  1  zL 

4000.  M 
1  t  0,145 


1  s  0,50 

2  * 

3  :  3.5,0 

4  l  S./  ,-270 

5  :  zL 

END  OF  EXECUTION 

CPU  TIME:  1.30  ELAPSED  TIME;  4:24.78 
EXIT 
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8.6  Option  4D — Differential  Energy  Rite  Contour  Diti 

This  option  functions  by  processing  two  files,  based  on  identical  altitude,  rnach  number 
and  load  factor  grids,  as  in  option  4A,  and  then  subtracting  the  data  at  the  grid  points.  An 
example  follows : 


OPTION  OR  (CR)  FOR  HELP  :  4D 

43  PS  DIFFERENTIAL  PLOT.  (P2.C0N,P2A.CQN,P2DIFF.CUN) 

DATA  BASE  FILENAME  :  D.NUH 

***  DATA  ARE  IN  METRIC  UNITS,  ENERGY  PARAMETER  IS  FS/UF  **♦ 

REPLY  "YES",  "NO",  "ALL"  OR  "END"  : 

OUTPUT  DATA  FOR  N  1.00  ?  Y_ 

OUTPUT  DATA  FOR  N  =  3.00  ?  N_ 

OUTPUT  DATA  FOR  N  =  5.00  ?_A^ 

COMPARISON  FILENAME  :  A.NUM 

END  OF  DATA  ON  L0G4 

END  OF  EXECUTION 

CPU  TIME:  1.78  ELAPSED  TIME:  47.54 
EXIT 


8.6.1  Data  Base  Fils  Name 


8.6.2  Output  Data  Base  for  N  =  n  r*si 

Replies  to  these  prompts  are  the  «ame  as  for  thos  e  of  option  4A.  Loa-  l actor  selection  then 
applies  both  to  the  base  file  and  to  the  comparison  file  considered. 


8.6.3  Comparison  Filename 

Reply  with  the  filename  of  the  data  file  to  be  compared  with  the  base  file.  Data  files  are 
considered  to  be  valid  for  comparison  if  the  following  parameters,  specified  when  running 
program  AIRCRAFT  to  create  the  data,  are  identical  for  both  files: 

H0,  H  STEP,  NO.  OF  S  '  EPS  (of  HP  or  ES) 

MACHO,  MACH  STEP  NO.  OF  STEPS  (of  M) 

GO,  G  STEP,  NO.  OF  STEPS  (of  G) 

IPSTYP,  IUNITS 

Altitude,  Mach  number  and  load  factor  grid  identitv  is  assured  if  the  PRESET  I.V. 
INCREMENTS  option  is  used  when  running  AIRCRAFT  to  create  both  data  files. 


B.6.4  Input  Errors 

(e.)  If,  for  any  reason,  tht  grid  parameters  are  not  identical,  thv.  message 
?  FILE  ARGUMENTS  ARE  DIFFERENT 
is  typed  together  with  a  listing  of  the  grid  parameters,  and  execution  finishes  immediately. 
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(6)  If  either  of  the  specified  input  files  does  not  exist;  an  error  message  is  typed  and  a  new 
filename  is  requested.  Supply  the  correct  name  to  continue  or  (fC)  to  return  to  the 
monitor. 

(c)  Replies  to  “?”  other  than  “Y”  are  interpreted  as  “N”. 

Notes: 

(1)  When  considering  whether  a  comparison  file  is  valid,  only  the  grid  parameters 
mentioned  above  are  checked.  Other  parameters,  such  as  power  setting,  wing  sweep 
and  atmosphere  type  are  not  checked. 

(2)  If  not  all  of  the  available  load  factors  are  selected  for  output,  the  informative 

END  OF  DATA  ON  LOG4 

is  typed  on  the  terminal.  This  is  not  an  error  condition,  but  simply  echoes  that  only 
selected  load  factors  will  appear  on  the  output  file. 

Output  for  this  option  is  in  ASCII  mode  on  three  files,  viz.  P2.CON,  P2A.CON  and 
P2DIFF.CON. 

P2.CON  is  the  same  as  would  be  produced  by  running  option  4A  with  the  same  replies, 
using  DATA  BASE  FILENAME  as  input. 

P2A.CON  is  the  same  as  would  be  produced  with  option  4A  using  COMPARISON 
FILENAME  as  input. 

P2DIFF.CON  is  the  file  for  the  differential  contour  plot. 

All  three  files  are  intended  as  input  to  program  P4,  but  may  be  printed  if  desired. 


8.7  Option  4E— Differential  MMD 

This  option  functions  by  processing  two  sets  of  data  and  subtracting  the  data  at  the  grid 
points  to  produce  values  of  differential  turn  rate.  An  example  follows: 


OPTION  OR  (CR)  FOR  HELP  :  4E 

4E  MHO  DIFFERENTIAL  PLOT.  ( P2 . OPT , P2A .OPT , P2D IFF . OPT ) 

ARE  BOTH  ".OPT"  FILES  ALREADY  ON  DISK  T  N_ 

ON-LINE  DATA  ?  N_ 

DATA  BASE  FILENAME  :  E.MHP 

*M  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARA*1::;:?  IS  F‘S  **+ 

DEFINE  ENERGY  GRID  -  (POSITIVE  INCREMENT) 

START,  STEP,  NO  OF  STEPS:  -1200,100,25 

ON-LINE  COMPARISON  DATA  ?  £ 

COMPARISON  FILENAME  :  C.MHD 

CPU  TIME  USED  =  0  MINS  1.56  SECS 
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8.7.1  Are  Both  “.OPT”  Files  Already  on  Disk? 

Answers  “Y”  if  comparison  files  P2.0PT  and  P2A.OPT,  produced  by  prior  runs  of  program 
P2,  are  to  be  used  as  input  data.  This  facility  is  included  to  enable  comparison  of  files  generated 
using  on-line  data,  without  the  need  to  repeat  the  on-line  dialogue  to  recreate  the  comparison 
files.  An  answer  “Y”  results  in  completion  of  execution  without  any  further  user  interaction. 
Execution  proceeds  normally,  with  interruption  only  if  energy  state  or  energy  rate  data  grids 
differ. 

Answer  “N”  if  on-line  data  or  program  AIRCRAFT  output  files  are  to  be  used  as  input 
data. 


8.7.2  On-line  Data 

8.7.3  Data  Base  Filename 

8.7.4  PS0,  PS  Step,  Number  of  Steps 

Reply  to  these  three  prompts  as  for  option  4C.  The  energy  rate  grid  chosen  applies  for  both 
the  base  file  and  the  comparison  file. 


8.7.5  On-line  Comparison  Data 

Comparison  data  may  also  be  in  a  disk  file  or  supplied  on-line  from  the  terminal.  Reply 
“Y”  or  “N”  as  required. 

8.7.6  Comparison  Filename 

This  prompt  is  given  if  data  is  not  being  supplied  on-line.  Reply  with  the  name  of  the  file 
of  comparison  data.  Data  files  are  considered  valid  for  comparison  if  the  energy  state  grid  para¬ 
meters  are  valid  for  both  sets  of  data. 

8.7.7.  Input  Errors 

(a)  If  the  grid  definitions  are  not  identical,  the  text 

?  FILE  ARGUMENTS  ARE  DIFFERENT 
is  typed,  followed  by  the  arguments  for  the  two  grids;  execution  then  terminates. 

(b)  If  either  of  the  specified  input  files  cannot  be  found  on  disk,  an  error  message  is  typed, 
and  a  new  filename  is  requested.  Supply  the  correct  name  to  continue  or  fC  to  abort 
the  job  and  return  to  the  monitor. 

( c )  Energy  rate  grid— a  valid  combination  ensures: 

(i)  increment  value  >  0;  and 

(ii)  0  <  number  of  values  ^  50. 

An  invalid  combination  causes  the  text 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed  on  the  terminal,  and  the  prompt  is  repeated. 

(d)  On-line  input  errors — remarks  made  in  Section  8.5.7  also  apply  for  on-line  differential 
data. 

Output  for  option  4E  is  in  ASCII  mode  on  three  files,  viz.,  P2.0PT,  P2A.OPT  and 
P2DIFF.OPT,  containing  data  for  base,  comparison  and  differential  MMD  contour  plots, 
respectively.  Each  is  intended  as  input  to  program  P4,  but  may  be  printed  if  required. 
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8.8  Option  4F —Maximum  Manoeuvre  Persistence  Aid 

This  option  supplies  an  on-line  calculation  aid  in  determining  manoeuvre  capability  at  a 
number  of  ranges  from  base.  An  example  follows: 


OPTION  OR  (CR)  FOR  HELP  I  4F 

4F  NHP  CALCULATION  ASSISTANCE. 
IMPERIAL  <0)  OR  S.I.  (I)  UNITS  ♦  _0 

NO  OF  RANEES 
«  A 


FUEL  VECTOR  -  LB 
■*  4133 

*  3216 

*  2220 
*  1180 


ES  -  FT,  OH  -  DEG/S, 
*  5000,24.014,38710 

UF  - 

LB/HR 

TURNSi  16.99  13.15 
*  10000,22.376,56073 

9.08 

4.83 

TURNSi  16.57  12.83 
*  13000,20.492,52263 

8.86 

4.71 

TURNS:  16.28  12.60 
♦  20000,18.776,49240 

8.70 

4.62 

TURNSi  15.84  12.26 

8.47 

4.50 

«  23000,17.233,46319 

TURNSi  15.40  11.93 

8.23 

4.38 

*  30000,15.523,46133 

TURNS:  13.97  10.82 
■*  33000,12.457,48465 

7.47 

3.97 

TURNSi  10.67  8.27 
«  40000,9.106,37840 

3.71 

3.03 

TURNSi  6.54  3.06 
♦  43000,8.071,49598 

3.30 

1 .86 

TURNSi  6.76  5.23 
«  30000,6.990,46831 

3.61 

1.92 

TURNSi  6.20  4.80 

3.31 

1.76 

,  > 


STOP 

END  OF  EXECUTION 

CPU  TIHEi  0.99  ELAPSED  TIN  5(17.44 
EXIT 
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8.5.1  Imperial  (0)  or  SI  (1)  Units 

Reply  0  or  1  to  indicate  Imperial  or  SI  unit  systems,  respectively. 

8.8.2  Number  of  Ranges 

Reply  with  the  number  of  ranges  for  which  available  fuel  quantities  are  known. 

8.8.3  Fuel  Vector — kg  (lb) 

Reply  with  the  values  of  fuel  availability  each  at  range  from  base,  one  value  per  line. 

The  values  should  be  at  distances  either  increasing  or  decreasing  from  base. 

8.8.4  ES— m  (ft),  OM— deg/s,  WF— kg/s  (Ib/hr) 

Reply  with  values  of  optimum  turn  rate  and  fuel  flow  rate,  in  the  units  indicated,  at  each 
energy  state  for  which  turning  data  is  required,  three  values  per  line,  in  response  to  the  asterisk 
prompt. 

The  program  responds  by  typing  the  numbers  of  turns  possible  using  the  given  fuel  vector. 
Execution  terminates  by  replying  to  the  prompt  with  a  carriage  return. 

8.8.5  Input  Errors 

(a)  If  the  number  of  range/fuel  availability  pairs  is  greater  than  50,  or 

( b )  If  the  value  of  fuel  flow  rate  is  not  positive, 
the  error  text 

INVALID  INPUT!  TRY  AGAIN 
is  typed  on  the  terminal  and  the  prompt  is  repeated. 

8.9  Option  4Z — Overview  of  Data  Grid 

This  option  presents,  in  a  concise  form,  the  fundamental  data  representation  for  the  energy- 
manoeuvrability  method.  Turn  rate  and  energy  parameter  are  plotted  against  Mach  number  for 
several  values  of  load  factor.  One  page  of  graphs  is  produced  for  each  value  of  the  height  varaiable, 
which  may  be  energy  state  or  pressure  altitude.  An  example  follows : 

OPTION  OR  <CR)  FOR  HELP  :  jZ_ 

4Z  TURN  RATE,  PS  VS  HACH  FOR  GIVEN  HEIGHT.  (P2.PLT ) 

DATA  BASE  FILENAME  :  A. NUN 

■■***  DATA  ARE  IN  METRIC  UNITS,  ENERGY  PARAMETER  IS  PS/UF  *** 

REPLY  "YES",  “NO",  "ALL"  OR  "END"  s 
HEIGHT  =  0.0  M  ?  Jf 

G  ■  1 .00  ?  X 

G  =  3.00  ?  X 

G  «  5.00  ?  X 

G  *  7.00  ?  £ 

HEIGHT  =  10000.0  M  ?  A_ 

HEIGHT  *  20000.0  M  ?  _N 

END  OF  EXECUTION 

CPU  TINE:  4.45  ELAPSED  TIME:  42.92 
EXIT 
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8.9.1  Data  Base  Filename 

Reply  with  the  name  of  the  data  file  created  for  this  option.  In  operation,  any  data  file 
created  by  program  AIRCRAFT  will  be  able  to  be  plotted  except  those  created  for  optimised 
grids  (maximum  manoeuvre  diagrams). 

8.9.2  Height  «=  nnannn.n  m  (ft) 

8.9.3  G  -  nn.nn 

For  these  two  prompts,  replies  “YES”,  “NO”,  “ALL"  or  “END”  have  the  same  meanings 
as  comparable  replies  when  using  option  4B  (Sections  8.4.S  and  8.4.6).  This  range  of  replies 
allows  a  selection  to  be  made  from  the  data  available  on  the  input  file. 

The  output  file  P2.PLT  contains  the  requested  plots,  and  is  submitted  to  the  plotter  queue 
in  the  usual  way.  Each  height  value  requested  produces  approximately  20  cm  of  plotted  output. 

Since  this  is  an  additional  plot  provided  for  an  overall  view  of  the  grid  for  an  aircraft,  no 
options  for  variation  of  plotting  scales  are  provided.  Each  group  of  curves  is  plotted  on  a  double 
graph  of  size  25  cm  by  14  cm  (10  in.  by  6  in.)  with  a  common  Mach  number  axis  (x  axis).  The 
scales  and  ranges  for  the  axes  are  shown  in  Table  5. 

8.9.4  Input  Errors 

(a)  In  correct  filename — the  correct  filename  is  requested,  as  in  other  options. 

( b )  Any  reply  to  a  height  or  load  factor  prompt  other  than  the  four  above  is  interpreted  as 
“NO”. 

9.  PROGRAM  P4  USER’S  GUIDE 
9.1  Loading  and  Saving  of  Absolute  Files 

The  commands  for  loading  and  saving  an  absolute  copy  of  program  P4  on  disk,  and  pro¬ 
ducing  a  symbol  map  are : 

■LQA/MAP  6P4LQAD 
FORTRAN:  P4 
MAIN. 

QUTXT 

P 

PLOTD 

P4HAIN 

LINK:  Loading 

EXIT 

.SAV 

P4  saved 


The  indirect  loading  command  is  employed  to  instruct  the  linking  loader  to  use  the 
loading  sequence  contained  in  the  disk  file  P4LOAD.  The  content  of  this  file  is  the  string 
P4, /SEARCH  P24LIB, C1033, 1Q223GRAFIC, EXTRAS, 6RAFIC 

File  P24LIB.REL  and  files  GRAFIC.REL  and  EXTRAS.REL  (both  on  disk  area  [1033,  1022]) 
are  searched  as  user  libraries. 

The  absolute  copy  of  the  program  is  saved  on  the  disk  file  P4.EXE,  which  is  then  executed 
using  the  RUN  command.  Core  storage  requirement  is  37  pages  (approximately  19Ks  words). 
The  storage  allocation  of  array  WORK  is  again  set  at  8000  words,  as  for  program  P2.  This 
storage  limit  will  need  changing  only  if  the  allocation  in  program  P2  is  changed. 
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9.2  Program  Execution 

The  inputs  to  program  P4  consist  of  data  files  created  in  standard  formats  by  program  P2. 
A  selection  is  then  made  using,  terminal  inputs,  of  data  to  be  plotted.  The  sole  output  of  the 
program  is  the  plotter  file  P4.PLT. 

A  classification  of  input  data  files  into  four  types  is  given  in  Section  6.2.  The  four  types  are 
energy  rate  contour  plots,  differential  energy  rate  contour  plots,  maximum  manoeuvre  diagrams, 
and  differential  maximum  manoeuvre  diagrams. 

The  following  sections  of  this  chapter  describe  the  terminal  operations  required  for  each 
type  of  data,  using  sample  dialogue.  User  input  to  terminal  prompts  is  again  underlined. 


9.3  Energy  Rate  Contour  Plots 
Example : 

RUN  P4 


INPUT  FILENAME  :  P2.CQN 

***  DATA  ARE  IN  METRIC  UNITS,  ENERGY  PARAMETER  IS  PS/UF  *** 
CONTOUR  PLOTTING 

SCALES  IN  UNIT5/IN  OF  PLOT  -  X,  Y  :  0.4,4000 
SMOOTHED  AND  TEXTURED  CONTOURS  ? 

PLOT  ES  CONTOURS  ?  Y 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  0,2000,25 
I’LOT  N  =  1.00  ?  _Y 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  0,5,15 
PLOT  N  =  3.00  ?  Y_ 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  -10,5,15 
PLOT  N  =  3.00  ?  N_ 

PLOT  N  *  7.00  ? 

CONTOUR  LEVELS  -  START,  STEP,  NO.  J  -30,5,150 

INVALID  INPUT!  TRY  AGAIN  - 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  -30,5,15 

PLOT  N  *  9.00  ?  £ 

STOP 

END  OF  EXECUTION 

CPU  TIME:  14.81  ELAPSED  TIME:  2:57.08 
EXIT 
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9.3.1  inpat  Filename 

Program  P2  produces  files  with  filenames 

P2.CON  or  P2A.CON 

for  single  aircraft  contour  plots  (the  latter  during  a  comparison  run).  One  of  these  two  names  is 
required  as  a  reply.  F4  responds  with  the  unit  system  for  the  data  and  the  type  of  energy  parameter. 

9.3.2  Scale  in  Unite/hL  of  Plot— x,  y 

Reply  with  the  amount  of  the  physical  quantity  which  one  inch  of  plot  would  represent  on 
each  axis.  For  energy  rate  contour  plots,  the  x-  and  y-axes  represent  Mach  number  and  altitude 
respectively. 

x-axis 

The  Mach  number  covers  the  range 

(Min,  Min+(No.  of  Machs— i)  x  Inc.) 

where  the  values  of  Min,  No.  and  Inc.  are  the  defining  grid  values  specified  when  running  program 
AIRCRAFT.  A  scale  between  0-2  and  0-4  unit/in.  is  usual,  such  that  the  axis  length  is  less  than 
20  cm  (8  in.)  for  adequate  A4  trim. 

y-axis 

Altitude  range  is  determined  by  a  similar  formula  to  that  above.  Usual  range  is  approxi¬ 
mately  0  to  20,000  m  (0  to  60,000  ft),  so  that  a  scale  of  4000  m/in.  (10,000  ft/in.)  provides 
adequate  A4  trim.  The  altitude  variable  will  normally  be  pressure  altitude,  but  it  is  also  possible 
to  plot  energy  state  as  the  altitude  variable,  in  which  case  the  scales  suggested  should  be  halved. 


9.3.3  Smoothed  and  Textured  Contours 

Reply  “Y”  or  “N”  as  required. 

The  smoothing  interval  of  0  05  in.  chosen  for  P4  results  in  rounding  off  the  sharp  corners 
where  the  otherwise  essentially  linear  contour  approximations  meet.  It  also  provides  a  visually 
pleasing  mark  :  space  ratio  for  contour  texturing,  which  cannot  be  chosen  unless  smoothing 
is  requested. 

For  most  plots  except  those  with  very  fine  grid  spacings,  smoothing  results  in  a  more 
attractive  plot  without  loss  of  accuracy,  at  the  expense  of  a  slight  increase  in  computing  time. 


9.3.4  Plot  ES  Contours 

Information  is  carried  in  the  data  file  to  enable  contours  of  constant  energy  state  to  be 
superimposed  on  the  requested  plots.  Reply  “Y”  if  these  are  required. 


9.3.5  Contour  Levels-Start,  Step,  Number 

This  prompt  is  given  to  define  energy  state  contours  (if  specified)  and  to  define  energy 
rate  contours  for  each  load  factor  specified.  For  energy  slate  contours  the  following  are  typical 
values : 

START:  2000  m  (5000  ft), 

STEP  :  2000  m  (5000  ft), 

NO.  :  up  to  30,  depending  on  ranges  for  x~  and  y-axes. 

For  energy  rate  contours ,  values  of  START,  STEP  and  NO.  should  cover  the  complete  range 
of  energy  parameters  expected  for  the  given  aircraft  configuration  and  load  factor.  Table  5 
gives  a  guide  to  the  ranges  for  the  energy  parameters  in  the  two  unit  systems.  Obviously  as  load 
factor  increases,  more  negative  contours  should  be  plotted. 
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9.3.6  Plot  N  -  x.xx 

The  prompt  will  supply  each  load  factor  in  turn  for  which  data  exists  on  the  data  file. 
Reply  “Y”  or  “N”  as  desired.  Contour  level  definition  will  then  be  requested  for  each  load  factor 
to  be  plotted. 

Depending  on  the  x-axis  scales  used.  20  to  35  cm  (8  to  14  in.)  of  plot  is  required  for  each 
load  factor  requested.  Submit  the  output  file  P4.PLT  to  the  plot  queue  in  the  usual  way  for 
plotting. 


9.3.7  Input  Errors 

(a)  Incorrect  filename. — An  error  message  is  typed  on  the  user’s  terminal  if  the  specified 
input  file  cannot  be  found  on  disk,  and  a  new  filename  is  requested.  Supply  the  correct 
name  to  continue  or  ( jC)  to  abort  the  job. 

(b)  Axis  scales. — If  either  scale  is  zero,  the  message 

INVALID  INPUT!  TRY  AGAIN 
is  typed,  and  new  scales  are  requested. 

(c)  PLOT  prompts.— Any  reply  other  than  “Y”  is  interpreted  as  “N”, 

(d)  Contour  level  specification. — A  valid  combination  ensures: 

(i)  increment  value  ^  0,  and 

(ii)  0  <  number  of  steps  <  50. 

An  invalid  combination  causes  the  text 

INVALID  INPUT!  TRY  AGAIN 
to  be  typed  on  the  terminal,  and  the  prompt  is  repeated. 


Notes: 

(1)  If  any  of  the  load  factors  specified  in  the  grid  definition  have  been  bypassed  by  program 
P2,  the  informative  text 


END  OF  DATA  CN  LOG5 

will  be  typed  on  the  terminal.  It  does  not  indicate  an  error  condition,  but  is  informative 
only. 

(2)  File  P4.PLT  must  be  submitted  to  the  plotter  queue  after  each  run  of  P4,  otherwise  the 
plot  may  be  accidentally  overwritten  the  next  time  P4  is  run. 
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9.4  Differential  Energy  Rate  Contour  Plots 


.RUN  P4 


INPUT  FILENAME  :  P2DIFF.C0N 

*»*  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

CONTOUR  PLOTTING 

SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  s  0.4,10000 
SMOOTHED  AND  TEXTURED  CONTOURS  ?  Y_ 

PLOT  ES  CONTOURS  ?  £ 

PLOT  N  =  1.00  ?  _Y_ 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  -500,50,16 
PLOT  N  =  3.00  ?  Y_ 

CONTOUR  LEVELS  -  START,  STEP,  NO.  !  -1000,100,16 
PLOT  N  =  5.00  ?  X 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  -1000,100,16 

PLOT  N  =  7.00  ?  X 

STOP 

END  OF  EXECUTION 

CPU  TIME:  11.44  ELAPSED  TIME:  2:25.24 
EXIT 


I 


Filename  P2DIFF.C0N  indicates  differential  energy  rate  data. 

Replies  for  this  type  of  run  &.re  very  similar  to  those  for  an  energy  rate  contour  plot  for  a 
single  aircraft.  The  only  difference  is  that  requests  for  contour  levels  now  refer  to  the  difference 
between  the  two  aircraft  at  each  grid  point.  The  aircraft  referred  to  by  program  P2  as  the  DATA 
BASE  will  be  indicated  by  positive  contours,  and  that  referred  to  as  the  COMPARISON  wilt  be 
indicated  by  negative  contours.  Full  identification-headers  are  included  on  the  plotted  output 
on  file  P4.PLT. 
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9.5  Maximum  Manoeuvre  Diagrams 

This  type  of  run  produces  a  single  graph  of  turn  rate  contours  on  an  optimum  energy  rate 
(y-axis)  vs.  energy  state  (x-axis)  grid. 


.RUN  P4 


INPUT  FILENAME  :  P2.QPT 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

CONTOUR  PLOTTING 

SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  !  10000,400 
SMOOTHED  AND  TEXTURED  CONTOURS  ?  N 
CONTOUR  LEVELS  -  START,  STEP,  NO.  :  0,2,16 
STOP 

END  OF  EXECUTION 

CPU  TIME:  6.75  ELAPSED  TiHE:  1:17.86 
EXIT 


9.5.1  Input  Filename 

Program  P2  produces  files  with  filenames 

P2.0PT  or  P2A.OPT 

for  a  single  aircraft  MMD  (the  latter  during  a  comparison  run).  One  of  these  two  names  should 
be  typed  as  a  reply.  P4  responds  with  the  unit  system  for  the  data  and  the  type  of  energy 
parameter. 


9.5.2  Scales  in  Units/in.  of  Piet — X,  ¥ 

Reply,  as  before,  with  the  amount  of  the  physical  quantity  which  one  inch  of  plot  would 
represent  on  each  axis.  For  MMD  plots,  the  x-  and  .y-axes  represent  energy  state  and  energy  rate, 
respectively. 

x-axis 

Energy  state  covers  the  range  defined  when  running  program  AIRCRAFT  (see  Section 
9.3.2).  The  usual  range  will  be  approximately  2000  to  24,000  m  (5000  to  75,000  ft),  so  that  a 
scale  of  4000  m/in.  (10,000  ft/in.)  provides  adequate  A4  trim. 

y-axis 

Energy  rate  covers  the  range  defined  in  the  usual  way  when  running  program  P2.  The  scales 
given  in  Table  5  may  be  used  as  a  guide  for  all  combinations  of  type  of  energy  parameter  and 
unit  system. 


9.5.3  Smoothed  and  Textured  Contours 

Reply  “Y”  or  “N”  as  required. 


9.5.4  Contour  Levels— Start,  Step,  Number 

Reply  with  parameters  defining  the  turn  rate  contours.  Typical  values  are: 
START:  0  deg/s; 

STEP  :  2  deg/s ; 

NO.  :  up  to  16. 


9.5.5  Input  Errors 

Input  errors  in  filename,  axis  scales  and  contour  levels  are  the  same  as  those  for  the  energy 
rate  contour  plots  (Section  9.3.7) 

After  execution,  output  file  P4.PLT  is  submitted  to  the  plotter  queue  in  the  usual  way. 
About  40  cm  (16  in.)  of  plot  will  be  produced  for  each  run. 


9.6  Differential  Maximum  Manoeuvre  Diagrams 


.RUN  P4 


INPUT  FILENAME  :  P2DIFF.0PT 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

CONTOUR  PLOTTING 

SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  :  10000,400 
SMOOTHED  AND  TEXTURED  CONTOURS  ?  Y _ 

CONTOUR  LEVELS  -  START,  STEP,  NO.  :  -10,2,11 
STOP 

END  OF  EXECUTION 

CPU  TINEs  4.34  ELAPSED  TIME:  1:13.52 
EXIT 


Filename  P2DIFF.OPT  indicates  a  differential  optimum  turn  rate  (MMD)  plot.  Replies 
for  this  type  of  run  are  similar  to  those  above  for  a  MMD  plot  for  a  single  aircraft.  Again,  the 
only  difference  is  that  contour  levels  now  represent  differences  between  two  aircraft,  and  the 
levels  will  cover  a  range  from  negative  to  positive,  rather  than  being  strictly  positive  as  for  a 
single  aircraft. 

Base  airciaft  superiority  will  be  indicated  by  positive  contours,  and  the  comparison  aircraft 
superiority  will  be  indicated  by  negative  contours.  Full  identification  headers  are  included  on 
the  plotted  output  on  file  P4.PLT. 
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10.  INPUT/OUTPUT  EXAMPLES 

10.1  Introduction 

Previous  chapters  have  discussc  I  fully  the  various  programs  in  the  suite.  This  chapter 
presents  examples  of  inputs  and  cutputs  when  using  the  programs  in  a  co-ordinated  fashion  to 
produce  tabular  summaries  and  contour,  turn  rate  and  maximum  manoeuvre  plots. 

10.2  File  Handling 

For  any  given  aircraft,  a  variety  of  input  and  output  files  will  be  used,  requiring  some  form 
of  standardised  housekeeping  for  efficient  storage  and  retrieval  of  files. 

Each  aircraft  is  represented  by  files  of  propulsion,  aerodynamic  and  configuration  data, 
together  with  source  and  absolute  versions  of  the  FORTRAN  programs.  These  files  are  con¬ 
veniently  stored  on  DECtape,  since  changes  to  one  or  more  files  will  be  needed  as  data  banks 
and  programs  are  developed. 

Output  files  for  line-printer  listing  or  plotting  may  be  preserved  or  deleted  as  required,  but 
it  is  useful  to  maintain  copies  of  files  needed  for  comparative  energy  rate  contour  plots  or 
maximum  manoeuvre  diagrams.  Once  produced,  these  tiles  are  not  changed,  so  it  is  convenient 
to  store  them  on  magnetic  tape.  A  suitable  naming  convention  is  required,  indicating: 

«  type  of  aircraft 

»  power  setting 

»  configuration  or  role 

•  wing  sweep  (if  applicable) 

•  atmosphere 

•  type  of  run  (MMD,  contour  plot,  etc.). 

Other  data  identifying  the  files,  such  as  unit  system,  type  of  energy  parameter  and  grid  parameiers, 
are  carried  within  the  files  and  used  by  program  P2  to  ensure  valid  comparative  data.  Thus 
there  is  no  necessity  to  identify  these  data  in  the  file  name. 

As  an  example,  the  filename  “AMXSUP.MMD”  might  be  used  to  identify  maximum 
manoeuvre  diagram  (MMD)  data  for  aircraft  A  operating  in  the  air  superiority  role  (SUP)  with 
maximum  power  (MX).  ICAO  Standard  Atmosphere  is  implied,  and  no  wing  sweep  variation 
is  needed.  A  standardised  data  grid  is  used  (discussed  below)  and  separate  magnetic  tapes  may  be 
used  for  data  in  the  two  unit  systems. 

10.3  Energy  Rate  Contour  Plots 

Figure  19  shows  a  sample  run  of  a  typical  AIRCRAFT  program  and  the  processing  of 
data  using  programs  P2  and  P4  to  produce  a  set  of  energy  rate  contour  plots. 

The  example  shown  uses  the  preset  data  grid  (see  Fig.  5(b)  or  Section  7.2).  Program 
AIRCRA  produces  a  file  line-printer  listing  (AMXSUP.LST)  and  a  file  for  input  to  progr  m  P2 
(AMXSUP.CON). 

Figure  20  shows  portion  of  file  AMXSUP.LST  and  illustrates  the  main  features  of  the 
listed  output.  The  header  page  identifies  the  run  and  echoes  configuration  data  supplied  hi  the 
configuration  file.  Subsequent  pages  group  energy  rate  data  and  turn  rate  data  for  each  specified 
altitude.  The  lowest  Mach  number  shown  for  each  altitude  is  the  lowest  Mach  number  specified 
by  the  data  grid  within  the  1  g  lift-limited  envelope.  No  structural  limits  are  indicated  at  higher 
Mach  numbers;  these  placard  limits  must  be  apptied  ofi-linc. 

The  file  AMXSUP.CON  contains  all  the  information  of  the  line-printer  listing  in  a  con¬ 
densed  format.  In  addition,  energy  rate  values  outside  the  lift-limited  envelope  are  flagged 
with  the  values  of  -9999-99. 

Figure  i9  shows  that  the  running  of  program  P2  using  option  4A  is  straightforward.  In  this 
case,  ail  load  factors  are  specified,  and  data  are  written  on  file  P2.CON. 

The  scales  used  when  running  program  P4,  and  the  ranges  of  contour  level  arc  discussed  in 
the  User’s  Guide  for  program  P4  (Section  9.3).  It  is  likely  that  the  first  attempt  o  plot  contours 
may  omit  some  conto.rs  of  interest,  and  this  trial  run  may  be  accelerated  by  omitting  the  time- 
consuming  smoothing  process. 
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.RUN  AIRCRA 


maneuveraiility  GRID  CALCULATION 


IMPERIAL  <•)  OR  8*1*  UNIT8 
AIRCRAFT  RATA  FILENAME 


r  • 

T  ASTORE 


AIRCRAFT  EXAMPLE 
AIR  SUPERIORITY  - 


DATE  3l-Apr-79  TINE  14i10  24.9 
2  AAN  ♦  BUN  ♦  3BX  FUEL 


PS,  PS/UF  OR  PS*V/UF  <1,  2  OR  3)  T  1 

IS  HP  THE  HEIBHT  VARIABLE  T  T 

PRESET  GRID  <Y,  N  OR  C/R)  T  1 

POWER  <HIL“100»NAX*20B)  *  5 

UIN8  SHEEP  (IF  VARIABLE)  J  , 

ATMOSPHERE  ,  DEVIATION  J  ] 

OUTPUT. (1-TEXT, 2-NOS. ,3-B0TH>  t  ; 

O/P  FILENAME  FOR  UNIT  A  T 

n/p  FILENAME  FOR  UNIT  8  T  ) 


T  20* 

T 

T  ICAO 
?  3 

T  AHXSUP.LST 
*  AMXSUP.CON 


CALCULATION  I 

ALTITUDE  HP  -  0.0  FT 

ALTITUDE  HP  *  4BBB.B  FT 

ALTITUDE  HP  *  8 001.0  FT 

ALTITUDE  HP  -  12BBB.B  FT 

ALTITUDE  HP  -  14I0I.0  FT 

ALTITUDE  HP  •  2BBOB.P  FT 

ALTITUDE  HP  «  24101 .f  FT 

ALTITUDE  HP  ■  2800B.0  FT 

ALTITUDE  HP  •  32000.#  FT 

ALTITUDE  HP  -  34l#t.»  FT 

ALTITUDE  HP  -  AiBfff.f  FT 

ALTITUDE  HP  *  44B0B.B  FT 

ALTITUDE  HP  »  48B0I.0  FT 

ALTITUDE  HP  *  32#Ol.f  FT 

ALTITUDE  HP  »  3410*. •  FT 

ALTITUDE  HP  ■  ABOOB.O  ft 

ALTITUDE  HP  ■  44000.#  FT 

ALTITUDE  HP  ■  48101.0  FT 

STOP 


END  OF  EXECUTION 
CPU  TIHEi  1 i 19.02 
EXIT 


ELAPSED  TIME*  5i59.42 


Fig.  19(a)  Sample  dialogue  for  energy  rate  contour  plot 
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.RUN  P2 


COMDAT  PERFORMANCE  PROCESSING 


OPTION  OR  <CR>  FOR  HELP  l  4A 

4  A  PS  CONTOUR  PLOT.  (P2.C0N) 

DATA  RASE  FILENAME  s  AMXSUP.CON 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS 

REPLY  "YES",  "NO",  "ALL"  OR  "END"  i 
OUTPUT  DATA  FOR  N  >  Kit  T  A 

END  OF  EXECUTION 

CPU  TINEt  14.39  ELAPSED  TINEl  Ii8.f2 
EXIT 

.RUN  PA 


INPUT  FILENAME  s  P2.C0N 

DATA  ARE  IN  IMPERIAL  UNITS,  ENER6Y  PARAMETER  IS  PS 
CONTOUR  PLOTTING 


SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  I  0.25,10000 
SMOOTHED  AND  TEXTURED  CONTOURS  ?  Y 


PLOT  ES  CONTOURS 

?  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

i 

5000,5000,25 

PLOT  N  -  1.00 

?  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

i 

0,100,12 

PLOT  N  -  3.90 

?  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

i 

-500,100,17 

PLOT  N  -  5.90 

?  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

l 

-1000,100,21 

PLOT  N  «  7.90 

T  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

i 

-1900,100,21 

PLOT  N  «  9.00 

?  Y 

CONTOUR  LEVELS  - 

START, 

STEP, 

NO. 

i 

-1500,100,21 

STOP 

END  OF  EXECUTION 

CPU  TINEl  57. fS  ELAPSED  TINEl  7:13.66 
EXIT 


Fig.  19(b)  Sample  dialogue  for  energy  rate  contour  plot  (cont.) 
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AIRCRAFT  MANEUVERABILITY  TABLE  FOR  AIRCRAFT  EXAMPLE 


Fig,  20{a)  Sample  listing  of  energy  rate  contour  data  -  header  page 
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Fig.  20{b)  Sample  listing  of  energy  rate  contour  data  -  extract 


\p\M\ 


The  \g  load  factor  energy  rate  contour  plot  produced  using  the  dialogue  in  Figure  19  is 
shown  in  Figure  21.  The  size  o?  the  plot  is  a  function  of  the  Mach  number  and  altitude  grids 
and  the  scale  specified.  The  different  contour  textures  available  for  energy  state  contours  and 
the  zero-energy  rate  contour  are  clearly  seen.  Contours  are  labelled  at  their  start  with  a  level 
number.  Contour  values  for  each  level  are  plotted  at  the  right  hand  side  of  the  plot;  units  are 
given  on  the  printed  output  (file  AMXSUP.LST).  Configuration  identification  is  given  at  the 
start  of  each  plotting  run  (not  shown  in  Fig.  21). 

Lift  limits  and  placard  limits  may  be  added  by  hand  (as  shown  on  Fig.  21).  Manual  retouching 
may  be  used  to  delete  contours  falling  outside  the  operating  envelope  or  to  smooth  contours  in 
the  vicinity  of  the  lift  limit  boundary. 


10.4  Differential  Energy  Rate  Contour  Plots 

The  same  sequence  of  running  programs  AIRCRAFT,  P2  and  P4  is  used  to  generate  dif¬ 
ferential  energy  rate  contour  plots,  as  used  in  the  previous  example. 

Data  files  of  comparative  data  for  input  to  program  P2  are  generated  by  running  programs 
for  each  aircraft;  alternatively  published  energy  rate  data  can  be  used  in  conjunction  with 
program  ANY,  as  shown  in  Figure  22.  Energy  rate  values  at  the  specified  increments  of  load 
factor  are  entered  on  each  line.  In  the  example  shown,  only  one  value  is  entered  per  line  since  the 
load  factor  grid  has  specified  only  the  1  g  load  factor  level.  Data  outside  the  lg  lift  or  structural 
limits  is  indicated  by  replying  to  the  prompt  with  a  carriage  return.  (This  requires  that  data 
values  of  zero  be  entered  as  small  non-zero  numbers,  e.g.  0  001). 

Option  4D  is  used  when  running  program  P2,  generating  input  files  to  program  P4  for  each 
aircraft,  and  a  third  file  of  differential  data.  An  example  of  a  differential  energy  rate  plot  produced 
by  running  program  P4  is  shown  in  Figure  23. 

Obviously,  differential  contours  only  have  meaning  within  the  operating  envelope  common 
to  both  aircraft.  Outside  this  envelope  (shown  by  hatching  in  the  figure),  contours  would  be 
deleted  by  hand  in  preparing  a  final  plot;  smoothing  is  also  required  in  the  vicinity  of  the  inner 
lift  limit  boundary. 


10.5  Turn  Rate  Plots 

Figure  24  shows  a  typical  set  of  terminal  dialogue  for  producing  turn  rate  plots  using 
programs  AIRCRAFT  and  P2. 

The  data  grid  shown  in  the  example  provides,  in  the  file  AMXSUP.TRN,  data  for  high 
subsonic  manoeuvres  up  to  9000  m  (approximately  30,000  ft).  A  load  factor  increment  of  0-25 
provides  sufficiently  smooth  curves  when  the  data  is  plotted.  Figure  25  shows  part  of  the  listing 
file  AMXSUP.LST,  and  illustrates  the  tabular  layout  for  the  33  load  factors  specified.  In  this 
example,  SI  units  are  nominated,  and  ali  output  headings  and  numerical  values  vary  accordingly. 

Four  combinations  of  Mach  number  and  altitude  are  selected  for  plotting  when  running 
program  P2  using  option  4B.  The  User’s  Guide  for  program  P2  (Section  8.4)  should  be  con¬ 
sulted  if  uncertain  of  the  scale  and  range  for  the  energy  rate  axis.  Each  turn  rate  curve  plotted  is 
identified  with  Mach  number  and  energy  state  annotation.  The  turn  rate  plot  for  3000  m 
altitude  from  the  above  example  is  shown  in  Figure  26. 


10.6  Maximum  Manoeuvre  Diagrams 

The  usual  sequence  of  running  programs  AIRCRAFT,  P2  and  P4  may  be  used  to  obtain 
maximum  manoeuvre  diagrams.  However,  the  optimization  process  involved  in  program 
AIRCRAFT  is  time-consuming,  and  on-line  running  of  this  program  can  be  replaced  by  batch 
running  using  the  ARL  Computer  Centre's  BATCON  system.17  This  simply  involves  grouping 
the  commands  required  to  run  AIRCRAFT  into  a  single  file  with  the  extension  *‘.CTL”. 

Such  a  file,  named  MMD.CTL,  is  shown  in  Figure  27.  The  commands  in  the  file  run 
AIRCRAFT,  supplying  all  the  conversational  replies  required,  prints  the  listing  file 
AMXMMD.LST,  and  leaves  the  input  data  for  program  P2  in  file  AMXSUP.MMD  on  the 
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Fig.  24  Sample  dialogue  for  turn  rate  plot 
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Fig.  25  Sample  listing  of  turn  rate  data  (extract) 
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Fig.  27  Batcon  commands  in  file  MMD.CTL 
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user’s  disk  area.  The  energy  state/'load  factor  grid  used  in  that  example  (15  states  from  5000  ft 
to  75,000  ft  and  17  load  factors  from  ig  to  9 g)  has  been  used  as  a  suitable  compromise  between 
computing  time  (approximately  6  minutes  CPU  time  for  the  example  shown)  and  interpolation 
accuracy.  Finer  grids  could  be  used  if  required. 

An  extract  of  file  AMXMMD.LST  for  the  above  example  is  shown  in  F'gure  28.  This 
extract  indicates  the  data  calculated  at  each  energy  state,  as  well  as  the  approximate  time 
histories  of  fuel  used  and  range  obtained  by  integration  from  the  initial  energy  state,  using 
Equations  (3.2)  to  (3.5). 

The  lg  load  factor  points  provide  the  Mach  number/altitude  schedule  for  the  optimum 
energy  climb.  The  right-hanu  load  factor  column  provides  optimum  sustained  turn  conditions 
for  use  in  maximum  manoeuvre  persistence  calculations.  k,<  the  example  shown,  cxtiapolaiion 
beyond  the  9g  data  limit  indicates  that  at  the  energy  states  shown,  the  sustained  turn  would  be 
load-factor  limited,  rather  than  thrust  limited.  The  listing  is  of  considerable  use  in  identifying 
limits  applied  to  the  final  plot  for  lift,  load  factor  and  airspeed  limitations. 

Running  programs  P2  (option  4C)  and  P4  is  straightforward,  as  outlined  in  the  User’s 
Guides  (Sections  8.5  and  9.5).  An  example  output  plot  from  program  P4  for  a  typical  modem 
generation  air-superiority  fighter  is  shown  in  Figure  29.  Lift  limits,  structural  limits  and  airspeed 
limits  are  applied  off-line  as  shown.  Contours  below  these  limits  represent  excursions  outside 
the  flight  envelope,  and  are  deleted  by  hand.  Contour  levels  are  labelled,  and  aircraft  con¬ 
figuration  identified  in  the  usual  way  with  text  at  the  right-  and  left-hand  side..-  of  the  plot  (the 
latter  not  shown). 

In  the  example  shown,  a  small  region  in  the  left-hand  corner  of  the  plot  shows  where  the 
grid  is  too  coarse  for  linear  interpolation  to  give  satisfactory  results.  This  situation  could  be 
remedied  by  using  a  finer  grid  for  the  whole  plot,  or  by  calculating  data  for  the  suspect  region 
on  a  very  fine  grid  and  adjusting  the  plot  by  hand. 


10.7  Differeuttal  Maximum  Manoeuvre  Diagrams 

Generation  of  MMD  data  files  for  aircraft  comparisons  is  performed  using  program 
AIRCRAFr  as  in  the  previous  section,  or  using  program  P2  to  supply  data  on-line  (see  Section 
8.5.5). 

Option  4E  is  used  when  running  program  P2.  In  the  example  shown  in  Figure  30,  data  fibs 
AMXSUP.MMD  and  BMXSUP.MMD  were  generated  using  the  BATCON  system  in  separate 
runs.  MMD  data  file1-  for  both  aircraft,  as  well  as  for  their  comparison,  arc  created  by  running 
program  P2  and  the  example  shows  that  program  P4  can  be  used  with  any  of  the  three  data 
files  in  successive  runs,  so  long  as  copies  of  tiic  output  file  P4.PLT  are  made  using  a  different 
filename. 

Running  P4  is  a  simple  operation;  guidance  with  scales  and  contour  levels  is  given  in  the 
User's  Guide  (Sections  9.5  and  9.6).  The  plot  produced  in  the  above  example  is  shown  in 
Figure  31.  The  zero  turn  rate  boundaries  for  the  two  aircraft  are  plotted  by  program  P4;  lift, 
load  factor  and  airspeed  limits  for  the  two  aircraft  must  be  added  by  hand.  Again,  contours 
exterior  to  the  common  flight  envelope  (shown  hatched)  have  no  significance,  and  would  be 
deleted.  Some  touching  up  would  also  be  necessary  in  the  immed:atc  vicinity  of  the  lower  turn 
rate  boundary,  if  the  plot  were  to  be  required  for  publication.  That  the  area  of  comparison  is 
so  small  indicates  the  vast  superiority  of  aircraft  “A”  (1980  generation  fighter)  over  aircraft 
“B”  (1960  generation  fiighter). 


10.8  Overview  Plots  of  Data  Grid 

Programs  AIRCRAFT  and  P2  are  used  to  produce  overview  plots  of  the  data  covering 
any  selected  portions  or  the  altitude  (or  energy  state)/Mach  number/load  factor  grid. 

Input  files  to  program  P2  may  be  generated  for  particular  altitudes  of  interest,  or  program 
P2  may  be  used  to  select  altitudes  from  the  energy  rate  contour  plot  data.  Figure  32  is  an  example 
of  these  grid  plots  at  a  typical  combat  altitude.  Each  plotted  curve  is  identified  by  a  load  factor 
annotation. 
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Fig.  28  Sample  maximum  manoeuvre  diagram  data  (extract) 


II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II 
~Nf9'ttn«>Ncoa>a*-Nro*mtD 


96 


ENERGY  STATE  FT  CX103 


RUN  P2 


CONIAT  PERFORMANCE  PROCESSING 


OPTION  OR  <CR)  FOR  HELP  t  4E 

4E  NMD  DIFFERENTIAL  PLOT.  <P2.OPT,P2A.0?T,P2DIFF.OPT> 

ARE  BOTH  ".OPT"  FILES  ALREADY  ON  DISK  T  N 

ON-LINE  DATA  ?  N 

DATA  BASE  FILENAME  t  AHX3UP.NND 

*M  DATA  ARE  IN  IHPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS 

DEFINE  ENERGY  GRID  -  < POSITIVE  INCREMENT) 

START,  STEP,  NO  OF  STEPSi  -12M,1M,2S 
ON-LINE  COMPARISON  DATA  ?  N 

COMPARISON  FILENAME  t  BMXSUP.HMD 

CPU  TINE  USED  *  f  MINS  5.94  StCS 

.RUN  P4 


INPUT  FILENAME  l  P2.0PT 

***  DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARAMETER  IS  PS  *** 

CONTOUR  PLOTTING 

SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  :  1««ff,4lf 
SMOOTHED  AND  TEXTURED  CONTOURS  ?  Y 
CONTOUR  LEVELS  -  START,  STEP,  NO.  i  #,2,16 
STOP 

END  OF  EXECUTION 

CPU  TIME!  9.21  ELAPSED  TIME*  1*23.32 
EXIT 

.COPY  HMD.PLT-P4.PLT 
.RUN  P4 


INPUT  FILENAHE  *  P2DIFF.0PT 

DATA  ARE  IN  IMPERIAL  UNITS,  ENERGY  PARANETER  IS  PS  •** 

CONTOUR  PLOTTING 

SCALES  IN  UNITS/IN  OF  PLOT  -  X,  Y  t  1«Mf,4f« 

SMOOTHED  AND  TEXTURED  CONTOURS  T  Y 
CONTOUR  LEVELS  -  START,  STEP,  NO.  *  -ff,2,2l 
STOP 

END  OF  EXECUTION 

CPU  TINE*  11.44  ELAPSED  TIME*  1*2f.S4 
EXIT 

.COPY  DIFMMD.PLT»P4.PLT 

Fig.  30  Sample  dialogue  for  differential  MMD  plots 
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Fig,  31  Sample  differential  MMD 
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HT  =  10000.  M 


As  with  other  plots,  lift,  load  factor  and  airspeed  limits  may  be  added  by  hand  to  give  the 
aircraft’s  flight  envelope  at  that  altitude,  The  Ps  =  0  axis  may  be  used  to  plot  the  sustained 
turn  boundary  on  the  upper  graph,  and  the  lift  limit/load  factor  limit  intersection  provides  an 
estimate  of  t'.ie  “corner  velocity”. 

Although  intended  as  an  overview  plot  when  the  suite  of  programs  was  designed,  plots  on 
this  format  have  proven  to  be  of  use  to  operational  pilots,  and  the  format  could  be  further 
developed. 


II.  FUTURE  PROGRAM  DEVELOPMENT 

11.1  Introduction 

The  development  of  the  suite  of  programs  described  has  been  completed,  and  the  range  and 
type  of  outputs  produced  satisfy  the  requirements  of  the  combat  performance  evaluation  task. 
The  quality  of  output  is  superior,  for  example,  to  that  produced  by  the  specific  excess  power 
programs  developed  for  use  with  the  Langley  differential  manoeuvring  simulator.18 

Nevertheless,  it  is  always  possible,  given  the  need  and  the  time,  to  improve  a  working 
product,  and  the  suite  described  is  no  exception.  This  chapter  discusses  several  areas  in  which 
accuracy,  presentation  and  ease  of  operation  could  be  improved  if  the  need  arises.  Any  improved 
accuracy  resulting  from  the  use  of  improved  computational  techniques  should  always  be  weighed 
against  the  errors  inherent  in  the  data  presentation  and  representation. 

11.2  Accuracy  Improvements 

There  are  two  areas  where  variations  in  numerical  methods  could  give  improved  accuracy, 
albeit  at  the  expense  of  greater  complexity  and  reduced  operating  speeds. 

The  first  area  is  in  the  use  of  higher  order  interpolation  schemes,  particularly  in  the  maximum 
manoeuvre  options  of  program  P2.  At  present,  linear  interpolation  is  used  and  a  finer  grid  may 
be  used  to  overcome  irregularities  in  the  plotted  output.  Routines  such  as  SURF  are  already 
included  in  the  subroutine  library  PI  LIB  and  little  extra  programming  effort  would  be  needed. 

The  second  area  where  improved  numerical  methods  may  be  useful  is  in  the  integration  of 
the  energy  function  to  provide  time  histories  of  optimum  climb  profiles. 

The  present  simple  Euler  method  could  be  replaced  by  an  algorithm  incorporating 
Simpson’s  rule,  but  in  this  case  the  improvements  in  calculation  accuracy  are  limited  by  model¬ 
ling  considerations.  The  energy  state  approximations  assume  that  potential  and  kinetic  energy 
interchanges  (i.e.  “zoom”  dives  or  climbs)  at  constant  energy  state  occur  in  zero  time,  and  do 
not  take  into  account  short  period  puli-ups/push-overs  required  in  the  transition  from  accelerated 
flight  to  constant  energy  flight.  These  factors  are  considered  in  more  detail  by  Spillman.7 

11.3  Presentation  Improvements 

Two  aspects  of  presentation  could  be  improved  by  further  development. 

At  present,  operating  envelopes  for  particular  aircraft  are  added  to  plotted  output  off-line. 
Lift  limits  are  obtainable  from  plotted  or  printed  output  by  extrapolation.  Structural  and  airspeed 
limitations  are  obtainable  only  from  data  sources,  All  limits  are  configuration-dependent  and 
thus  conveniently  determined  and  applied  by  hand  after  each  run.  It  is  possible  to  include  these 
limits  with  the  configuration  data  file,  but  this  will  involve  programming  changes  to  plotting 
routines  to  control  pen  operation,  as  well  as  careful  monitoring  of  input  files  as  configurations 
are  changed. 

A  second  presentation  change  which  could  be  needed  as  service  personnel  become  familiar 
with  the  plotted  outputs,  is  to  extend  the  grid  option  of  program  P2  (option  4Z)  by  plotting 
energy  rate  contours  on  the  Mach  number/turn  rate  axes  (i.e.,  combine  the  two  plots  shown  in 
Fig.  32).  The  changes  required  to  do  this  are  not  trivial,  but  the  benefits  in  presentation  may 
justify  the  programming  effort  involved.  The  modular  structure  of  program  P2  minimizes  the 
problems  associated  with  such  a  change,  and  facilitates  the  addition  of  any  future  options  which 
might  be  thought  necessary. 
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11.4  Operation  Improvement* 

A  recent  addition  to  the  PDP-10  system  facilities  at  ARL  is  the  MIC  system  for  on-line 
execution  of  commands  in  a  user-defined  control  file. 

Changes  could  be  made  to  all  the  programs  in  the  ARL  suite  to  take  advantage  of  this 
system  program,  particularly  for  production  running,  where  terminal  responses  vary  little  from 
run  to  run.  The  changes  would  involve  writing  the  user  responses  to  the  programs  in  a  series  of 
control  files,  which  the  MIC  system  is  then  commanded  to  process.  User  interaction  is  then 
minimal,  ensuring  error-free  production  running.  The  amount  of  program  changes  required 
depends  only  on  the  extent  to  which  the  MIC  system  is  involved. 


12.  CONCLUSION 

A  suite  of  FORTRAN  IV  computer  programs  has  been  developed  for  the  computation  and 
presentation  of  data  used  in  evaluating  combat  aircraft  performance  using  energy  manoeuvr¬ 
ability  theory. 

The  arrangement  of  programs  and  subroutine  libraries  allows  for  flexibility  of  data  repre¬ 
sentation  of  comparison  aircraft,  producing  a  variety  of  printed  and  plotted  outputs.  The  pro¬ 
grams  have  been  fully  described,  presenting  self-contained  user’s  guides  and  reference 
documentation. 
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APPENDIX  1 

Newton’s  Method  Iteration  for  Calculating  Geopotential  Height 
Energy  state  (£s)  is  expressed  in  terms  of  geopotential  height  ( h )  and  true  airspeed  (V)  as: 

£s  =  E/W  m  h  +  V*llg  (A.l) 

or,  since  both  h  and  V  are  functions  of  pressure  height  (Ap): 

£s  =A(M+m*p)}a/2g.  (A.2) 

By  definition,  Mach  number  ( M )  is  related  to  true  airspeed  by: 

M  =  Via  =  K/(y*lr£»itra)°»  (A.3) 

where  a  =  sonic  speed,  m/s  (ft/s), 

y»ir  =  specific  heat  ratio  for  air,  1  -4, 

R%\t  =  gas  constant  for  air,  287  055  3/kg.K  (3089-78  fta/K.s2), 

7’a  r-=  ambient  temperature,  K. 

Solving  Equation  (A.3)  for  V  and  substituting  in  Equation  (A.2)  yields  £s  in  terms  of  hp: 

£s  -  /i(/»n)+{yatf«.irA/?/(2g)} .  T>(hv).  (A. 4) 

When  £s  and  M  are  specified,  Equation  (A.4)  represents  a  non-linear  equation  in  /»p,  which 
may  be  solved  using  Newton’s  method  (see  example  tef.  19,  pp.  92,  157-160). 

Let 

k  -  y*(r£alr.W2/(2g)  (A. 5) 

which  is  a  constant  for  a  given  M,  and  rewrite  Equation  (A.4)  as  an  equation  in  /»p,  viz. 

«(/»p)  -  h{hp)+kT^hP)  -Es  -  0.  (A. 6) 

Let  /ip,/  denote  the  t'th  approximation  to  /ip.  Then,  using  Newton's  method, 


,  ,  «(V<) 

Zlp.M  -  «p,<  — 

«  (»P.f) 


(A.7) 


where 


«'(^P.O  —  Lt>e/bAp]Ap-=Ap,(  (A. 8) 

and  '  denotes  differentiation  w.r.t.  hv . 

Differentiating  (A.6)  w.r.i.  hv  gives 

«'(A„)  =  A'(^>P)+ A7'»'(AP).  (A.9) 

For  the  atmospheric  models  used,  it  can  be  shown8  that,  in  a  given  linear  segment  in  a  non¬ 
standard  atmosphere, 


t.e. 


k  (Zip)  —  (A/>P)  —  (7  ab  x  Aa/i/ip)/7b+AA/ip), 
A'(Ap)  =  2a(AP)/7Vfp), 


(A.  10) 
(A.  11) 


l . . 


where  atmospheric  temperature  ( T )  is  given  by 

T  *»  Tb+AAAp. 

Differentiation  of  Equation  (A.  1 2)  w.r.t.  Ap  gives 

r.  (Ap)  -  A.(AAp)'  -  A„ 


since 


(AAp)'  ■»  (hp—hp\})'  “  hp'  =*  i. 

Thus,  substituting  Equations  (A.i  i)  and  (A.  13)  into  Equation  (A. 9)  gives 

-  n(hp)IT(hP)  \ k\>, 

and  substituting  Equations  (A.  15)  and  (A.6)  into  Equation  (A.7)  gives 

ft  p.t+i  =  fip,<"l'8Ap1< 


where 

_h(hp,i)+kTtXhp'i)—E& 
p,<  =  _  r*(hp.i)iT(hp.i)-\-kkp * 

A  suitable  iterative  algorithm  can  now  be  written  as  follows: 


(A.  12) 

(A.  1 3) 
(A.  14) 

(A.  15) 


(A.  16) 


(1)  estimate  V  =  Es--asi,2M*/(2g).{  -ho)', 

(2)  set  i  =  1 ; 

(3)  use  subroutines  ATMOS  and  1NTRP  given  in  Reference  8  to  determine  h,  T»,  T  and  L; 

(4)  solve  (A.  16)  for 

(5)  if  8Ar<  <  tolerance  [say,  0  075  in  (0  25  ft)],  exit; 

(6)  increment  i  by  1  and  return  to  step  3. 

On  exit,  the  required  values  of  /ip,  h  and  7»  are  the  exit  values  of  these  parameters  as  deter¬ 
mined  by  subroutine  ATMOS. 
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APPENDIX  2 

“AIRCRAFT”  Program  Library 

A2.I  AIRCRAFT  Piogram 

Purpose 

This  is  a  small  routine  to  call  for  program  input  and  to  start  calculation  of  aircraft  perfor¬ 
mance  parameters. 

Use 

.RUN  AIRCRAFT  (monitor  command) 

User  routines  called 
PUN  TABLE 


A2.2  AERO  Subroutine  (example) 

Purpose 

AERO  is  m  aircraft-dependent  routine  which  determines  trimmed  aircraft  aerodynamic 
parameters,  using  curve-fitted  data  in  COMMON  EDATA.  An  initial  estimate  of  angle  of 
attack  is  required  as  input,  and  trimmed  angle  of  attack  is  returned  as  output.  Input  and  output 
are  via  COMMON  B. 

Use 

CALL  AERO 
Input 

GN  XM  W  FN  ALPHA  DIM  CGPERC  NDCDS 
XMDCDS  DELCDS 

(and  any  other  parameters  needed) 

Output 

ALPHA  D  IERR 
Routines  Called 

SURF  (and  any  other  aircraft-dependent  routines) 

Calling  Routines 
TABLE  SEP 

A2.3  BLOCK  DATA 

Purpose 

To  define  data  in  COMMON  C,  representing  atmospheric  profiles,  text  constants  and  other 
numeric  constants. 

Use 

Loaded  during  loading  phase. 

A2.4  THRUST  Subroutine  (example) 

Purpose 

THRUST  is  m  aircraft-dependent  routine  which  determines  engine  thrust  and  fuel  flow, 
using  curve-fitted  data  in  COMMON  FDATA.  An  auxiliary  output  is  the  lower  limit  of  Mach 
number  validity. 

Input  and  output  are  via  COMMON  B. 

Use 

CALL  THRUST 
Input 

XM  HP  PLA  DTEMP  MT 
(and  any  other  parameters  needed) 


Output 

FN  WF  XMLIM  IERR 
User  Routines  Called 

SURF  (and  any  other  aircraft-dependent  thrust  routines) 

Calling  Routines 
TABLE  SEP 
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APPENDIX  3 
PI  Subroutine  Library 

A3.1  ALTIT  Subroutine 

Purpose 

ALTIT  performs  a  Newton’s  method  iteration  to  determine  pressure  height  and  other 
atmospheric  quantities,  when  the  independent  variables  are  energy  height  and  Mach  number. 

Input  and  output  are  via  COMMON  B.  Atmosphere  profiles  defined  in  COMMON  C  are 
used. 

Use 

CALL.ALT1T 

lnpu> 

COMMON  B:  ES  XM 

COMMON  C:  TICAO  GAMMA  R  G 

Output 

COMMON  B:  T  TS  HP  H  ALAPSE  Cl  1ERR 

User  Routines  Called 

ATMOS  HEIGHT 

System  Routines  Called 
ABS 

Calling  Routines 
TABLE  SEP 

Comment 

This  subroutine  embodies  the  algorithm  of  Appendix  1 .  A  failure  exit  is  made  with  error 
indicator  IERR  set  to  201  if  more  thin  20  iterations  are  required  for  ES  convergence  to  within 
0-25  ft. 


A3. 2  ATMOS  Subroutine 


Purpose 

ATMOS  calculates  atmospheric  parameters,  giving  temperature  as  a  function  of  pressure 
height  in  either  ICAO  or  AROU  Atmospheres. 

Input  and  output  are  via  COMMON  B.  Atmosphere  profiles  defined  in  COMMON  C 
are  also  used. 

Use 

CALL  ATMOS 


Input 

COMMON  B: 
COMMON  C: 

Output 

COMMON  B: 


HP  MT  DTEMP 

HICAO  TICAO  HARDU  TARDU 

ALAPSE  TS  T  ALPSTD  LEV 


User  Routine  Call 
INTRP 


Calling  Routines 

TABLE  ALTIT 

A3.3  BADENP  Subroutine 

Purpose 

BAD1NP  types  out  an  error  message  on  the  user’s  terminal.  It  is  used  while  checking  input 
validity  during  conversational  input  operations.  If  invalid  input  is  detected,  recovery  is  achieved 
by  repeating  the  input  prompt  after  return  from  BADINP. 

Use 

CALL  BADINP 
Calling  Routine 

pun 

A3.4  BININ  Subroutine 

Purpose 

BININ  opens  and  reads  data  into  a  nominated  storage  area  from  a  specified  disk  file.  This 
file  must  consist  of  a  sequence  of  binary  records,  with  the  first  word  in  each  record  being  the 
number  of  data  items  in  that  record. 

Use 

CALL  BININ  (A,  NA,  LIMA,  FILNAM) 

Input 
A 

NA 
LIMA 
FILNAM 

Output 

A  Array  containing  the  data. 

Calling  Routine 
IDENT 


A3.5  HEIGHT  Subroutine 

Purpose 

HEIGHT  calculates  the  geopotential  height  corresponding  to  a  given  pressure  height  in  the 
nominated  atmosphere. 

Input  and  output  are  via  COMMON  B.  Atmosphere  profiles  defined  in  COMMON  C  are 
used. 

Use 

CALL  HEIGHT 
Input 

COMMON  B:  HP  MT  DTEMP 

COMMON  C:  HICAO  HARDU  TICAO  TARDU 

Output 

COMMONS:  H 

User  Routines  Called 
1NTRP  HTRUE 


The  array  into  which  data  is  to  be  read 
The  index  of  the  first  location  in  A 
The  dimension  of  A 

A  double  precision  text  variable  or  literal  giving  the  name  of  the  disk  file. 


System  Routine  Called 
AMINI 

Calling  Routines 

TABLE  ALT1T 

A3.6  HTRUE  Function 
Purpose 

HTRUE  is  called  by  HEIGHT  to  integrate  over  constant  temperature  lapse  rates  to  obtain  a 
geopotential  height  increment,  using  the  method  of  Relcrence  8. 

Use 

HG  =  HTRUE  (A,  B,  C,  D,  H2,  HI) 

Input 

A  Ambient  lapse  rate 

B  ICAO  lapse  rate 

C  Sea  level  extrapolated  ambient  temperature 
D  Sea  level  extrapolated  ICAO  temperature 

H2  Upper  pressure  height 

H I  Lower  pressure  height 

Output 

HG  Geopotential  height  increment 

System  Routine  Called 
ALOG 

Calling  Routine 
HEIGHT 

A3.7  IDENT  Subroutine 

Purpose 

IDENT  is  an  input  subroutine  which  reads  the  aircraft  identification  data  file  and  calls 
BIN  IN  to  read  thrust  and  drag  data  files. 

Aircraft  identification  data  is  transmitted  via  COMMON  B,  while  aerodynamic  and  thrust 
data  are  transmitted  via  COMMON  EDATA  and  COMMON  FDATA. 

Use 

CALL IDENT 
Output 

COMMON  B:  3REF  W  FRAME  SWPMIN  SWPMAX  NDCDS 

XMDCDS  DELCPS  ROLE  CGPERC 
COMMON  EDATA:  E 

COMMON  FDATA:  F 

User  Routine  Called 
BININ 

Calling  Routine 
PUN 

A3.8  INTRP  Subroutine 
Purpose 

INTRP  is  a  linear  interpolation  routine,  returning  ordinate  and  gradient  information. 

Use 

CALL  INTRP  (TABT,  TABH,  DTDH,  HP,  TEMP,  N) 


Input 

TABT  Table  of  ordinate  values  (x<) 
TABH  Table  of  abscissa  values  (yi) 
HP  Independent  variable  (.v) 

N  Index  of  next  ,v<  >  x 

Output 

DTDH  Slope  (dy/dx) 

TEMP  Dependent  variable  (>>) 


Calling  Routine 

ATMOS  HEIGHT  MAXMAN 


A3.9  MAXMAN  Subroutine 

Purpose 

MAXMAN  controls  the  calculation  and  output  of  the  optimised  manoeuvrability  grid, 
using  energy  height  and  load  factor  as  grid  variables.  Energy  parameter  is  optimised  by  varying 
Mach  number,  and  a  simple  Euler  integration  is  performed  at  each  energy  state  to  estimate  a 
climb  time  history. 

Input  and  output  are  via  COMMON  areas  B,  TABLES,  TABLET,  ZEROPS  and 
CLIMBS.  Data  constants  are  input  via  COMMON  C. 

Use 

CALL  MAXMAN 


Input 

COMMON  B:  PLA  DELES 

1PSTYP 

COMMON  C:  G  RADIAN 


NES  DELGN  NGN  ES0  GN0 
HOUR 


Output 

COMMON  B: 

COMMON  TABLES 
COMMON  TABLET 
COMMON  ZEROPS 
COMMON  CLIMBS 

User  Routines  Called 

SEP  1NTRP  ROMIN  PIOUTA 


HP  IES 
GNTAB 

HPTAB  WFTAB 
XMPS0  VPS0  HPPS0  WFPS0 


MONSEP 


ES  GN  XM 
PSTAB  OMTAB 
XMTAB  VTAB 
GNPS0  OMPS0 
T1  FU  RA 


System  Routines  Called 

SQRT  ASIN  ABS  COS 


Calling  Routine 
TABLE 


A3. 10  MONSEP  Subroutine 
Purpose 

MONSEP  is  a  routine  called  by  ROMIN  to  monitor  the  convergence  of  the  energy  para- 
meter  optimisation.  The  convergence  criterion  to  be  satisfied  at  the  ith  optimisation  step  is 

\fi{XM)-f  ^XM)\  <  |I0  ■*/<(A»/)|. 

If  more  than  50  evaluations  of f{Xiu)  are  required  the  optimisation  is  terminated  by  setting 
the  convergence  flag  and  returning  to  ROMIN. 

All  input  and  output  are  via  the  arguments  of  the  subroutine  call,  except  for  input  parameter 
1FAIL,  which  is  transmitted  via  COMMON  B. 

Use 


CALL  MONSEP  (N,  AX,  F,  R,  BETA,  CON,  NR) 


'''T’T'i'SH  p i)  m _ _._.  _ 
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Input 

N  The  number  of  independent  variables  ( « 1 ,  in  this  program) 

AX  The  vector  of  independent  variables  (Mach  number) 

F  The  current  value  of  the  (negative)  energy  parameter 

R  The  actual  number  of  energy  parameter  evaluations 

BETA  The  value  of  the  Euclidian  norm  of  the  vector  representing  the  total  progress 
since  the  last  axis  rotation 

NR  A  monitor  index — 

0  indicates  initial  function  evaluation 

1  indicates  a  single  function  evaluation 

2  indicates  completion  of  an  iteration  stage 

Output 

CON  A  logical  variable  set  to  .TRUE,  if  the  convergence  criterion  has  been  satisfied, 
otherwise  .FALSE. 

System  Routine  Called 
ABS 

Calling  Routine 
ROM  IN 


A3.ll  PAR  AMS  Subroutine 

Purpose 

PARAMS  calculates  airspeed  ana  pressure  parameters,  and  sets  the  initial  angle  of  attack 
estimate  for  each  aircraft  trim  calculation. 

Input  and  output  are  via  COMMON  B,  and  atmospheric  constants  stored  in  COMMON  C 
are  used. 

Use 

CALL  PARAMS 
Input 

COMMON  B:  XM  T  TS  HP  SREF  LEV  ALPSTD  C2 

COMMON  C:  H1CAO  TICAO  PTAB  GAMMA  R  G 

Output 

COMMON  B:  P  A  V  ALPHA  DIM 

System  Routine  Called 

SQRT  EXP  ABS 

Calling  Routines 
TABLE  SEP 


A3.12  PIIN  Subroutine 
Purpose 

PUN  controls  all  input  data.  All  aircraft-independent  parameters  defining  the  type  of 
calculations  required  are  accepted  in  conversational  mode,  with  validity  checks  where  possible. 
Aircraft-dependent  data  is  read  from  disk  via  a  call  to  1DENT.  Output  files  are  opened  before 
returning  to  the  calling  program. 

Output  variables  are  transmitted  via  COMMON  B,  ard  constants  stored  in  COMMON  C 
are  used. 

Use 

CALL  P11N 
Input 

COMMON  C.  NIO  AMET  R  G  HOUR  GRID1 


WWW**  wmpnpqminipp 


I’^mrp^*  w* 


Output 

COMMON  B: 


r  r,’*T^rr^'H7r»':rTT\  r  y?-v»f  ,■*»*•  ••« 
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SWP  PJLA  DTEMP  MT  C2  DELES  NES 
DELXM  NXM  DEIGN  NGN  ES0  XM0 
GN0  IOFT  FRAME  DAY  CLOCK  IOUT 
IGRID  ROLE  UNITS  AM  ALB  AWF 


User  Routines  C ailed 
IDENT  BADINP 


System  Routines  Called 
DATE  TIME 


Calling  Routine 

AIRCRAFT  program 


AT  13  JPIOUT  Subroutine 
Purpose 

PI  OUT  provides  output  control  for  unoptimized  performance  calculations.  A  text-formatted 
output  file  is  produced  for  printing,  an  alphanumeric  output  file  is  produced  as  input  to  program 
P2,  and  altitude  variable  is  output  to  the  user’s  terminal  to  indicate  progress  of  the  calculation. 
Output  units  may  be  Imperial  or  SI,  as  selected  by  the  user  during  the  input  dialogue. 

Input  of  data  to  be  printed  is  via  COMMON  B  and  COMMON  TABLES.  Data  constants 
used  are  transmitted  via  COMMON  C. 

Use 

CALL  Pi  OUT 
Input 

COMMON  B:  ES  XM  HP  V  H  SWP  PLA  SREF  W 

DTEMP  MT  WF  DELES  NES  DELXM  NXM 

DELGN  NGN  ES0  XM0  GN0  FRAME  1ES 

IXM  DAY  CLOCK  IOUT  IGRID  NDCDS 

XMDCDS  DELCDS  ROLE  IPSTYP  CGPERC 

UNITS  IPSTYP  AM  ALB  AWF 
COMMON  C  NIO  AMET  GRID1  GR1D2 
COMMON  TABLES:  PSTAB  OMTAB  GNTAB 

System  Routines  Called 
FLOAT  MIN0 

Calling  Routine 
TABLE 

A3.I4  PIOUTa  Subroutine 
Purpose 

PIOUTA  provides  output  control  for  optimized  performance  calculations,  in  a  like  manner 
to  subroutine  PIOUT. 

Input  of  data  to  be  printed  is  via  COMMON  areas  B,  TABLES,  TABLET,  ZEROPS  and 
CLIMBS.  Data  constants  used  are  transmitted  via  COMMON  C. 

Use 

CALL  PIOUTA 
Input 

COMMON  3:  As  for  subroutine  PIOUT 

COMMON  C:  NIO  AMET 


COMMON  TABLES: 

As  for  subroutine  PIOUT 

COMMON  ZEROPS: 

GNPS0 

WFPS0 

OMPS0 

SMPS0 

VPS0  HPPS0 

COMMON  CLIMBS: 

TI  FU 

RA 

COMMON  TABLET: 

XMTAB 

VTAB 

HPTAB 

WFTAB 

«wpi mr 


. . 


System  Routine  Called 
FLOAT  MINO 

Calling  Routine 
MAXMAN 


A3. 15  SEP  Subroutine 

Purpose 

SEP  calculates  the  energy  parameter  at  a  given  grid  point,  for  use  with  the  optimisation 
routine  ROM1N. 

Input  and  output  are  via  the  subroutine  arguments  listed  below.  During  calculation  variable 
data  is  transmitted  between  subroutines  via  COMMON  B.  Data  constants  in  COMMON  C  are 
also  used. 

Use 

CALL  SEP  (N,  X,  F) 

Input 

N  The  number  of  independent  variables  ( =  1 ,  in  this  program) 

X  The  vector  of  independent  variables  (Mach  number) 

Output 

F  The  current  value  of  the  (negative)  energy  parameter 
COMMON  Input 

COMMON  B:  W  IERR  IPSTYP 

COMMON  C:  HICAO  HOUR 

COMMON  Area  Output 

COMMON  B:  XM  HP  V  FN  WF  ALPHA  D  IERR 

User  Routines  Called 

ALTIT  PARAMS  THRUST  AERO 

System  Routine  Called 
COSD 

Calling  Routines 

MAXMAN  ROM  IN 


A3.16  TABLE  Subroutine 

Purpose 

TABLE  controls  the  calculation  and  output  of  the  unoptimised  manoeuvrability  grid,  using 
altitude,  Mach  number  and  load  factor  as  grid  variables.  (The  altitude  may  be  energy  height  and 
pressure  height,  depending  on  the  input  request.)  The  principle  outputs  are  energy  parameter  and 
turn  rate. 

Input  and  output  are  via  COMMON  areas  B  and  TABLES.  Data  constants  stored  in 
COMMON  C  are  also  used. 

Use 

CALL  TABLE 
Input 

COMMON  B:  DELES  NES  DELXM  NXM  DELGN  NGN 

ES0  XM0  GN0  IOPT  IGR1D  IPSTYP  AM 
ALB 
HICAO 


COMMON  C: 


Output 

COMMON  B:  ES  GN 

ALPHA  D 
COMMON  TABLES:  PSTAB 


XM  HP  V  H 
IES  IXM  IERR 
OMTAB  GNTAB 


W 


FN  WF 


User  Routines  Called 

MAXMAN  ATMOS  HEIGHT  ALTIT 
AERO  PIOUT 


PARAMS  THRUST 


System  Routines  Called 
SQRT  COSD 


Calling  Routine 

AIRCRAFT  program 


APPENDIX  4 


P1LIB  Subroutine  Library 

The  curve-fitting  techniques  used  in  preparing  thrust  and  drag  data  are  fully  described  in 
Reference  9;  a  brief  description  of  the  polynomial  evaluation  routines  is  given  below. 

A4.1  CHECKD  Subroutine 

Purpose 

CHECKD  checks  the  validity  of  data  supplied  to  subroutine  CUBICS  for  curve-fitting.  In 
the  present  application,  CUBICS  is  called  to  fit  a  simple  cubic  to  four  points  by  subroutine 
SURF,  and  so  the  parameter  I  FAIL  should  always  return  a  zero  value. 

Use 

CALL  CHECKD  (M,  NCAP,  X,  W,  K,  1FAIL,  NCAP3,  B) 

Input 

M  The  number  of  data  points  (here,  4) 

NCAP  The  number  of  intervals,  n  (here,  1) 

X  The  data  abscissae  (x(j),  j  —  1 ,  M) 

W  The  weight  vector  (w(j),j  =  1 ,  M)  (here  w(j)  -  1,  all  j) 

AK  The  vector  of  knots  ( ak(j ),  j  =  —3,  4  +  3).  (Here  there  are  no  external  knots, 

so  only  the  storage  space  for  internal  knots  is  required) 

NCAP3  An  integer  =  n+3 
B  A  working  vector  ( b(j),j  =  1,  n) 

Output 

IFAIL  The  error  indicator — 

=  0  for  valid  data 

=  1  to  6  for  invalid  data  (see  subroutine  CUBICS) 

Calling  Routine 
CUBICS 

A4.2  CUBICS  Subroutine 
Purpose 

CUBICS  applies  a  least-squares  cubic  spline  fit  of  B-splines  to  weighted  data  points  with 
selected  knots.  In  the  present  task,  CUBICS  is  only  required  to  fit  a  simple  cubic  to  four  points; 
hence  the  fit  is  exact,  all  weights  are  unity,  and  no  external  knots  are  required. 

Use 

CALL  CUBICS  (M,  NCAP,  X,  Y,  W,  AK,  C,  SS,  IFAIL,  NCAP3,  A,  DIAG,  B) 


Input 

M 

NCAP 

X 

Y 

W 

AK 

NCAP3 

A 

DIAG 

B 


The  number  of  data  points  (here,  4) 

The  number  of  intervals  (here,  1) 

The  data  abscissae  (x(j),  j  *>  1,  M) 

The  data  ordinates  (>'0'),  j  =  1 ,  M ) 

The  weight  vector  (vt(y ),  j  =  1,  M)  (here  w(J)  =  1,  ally) 
The  vector  of  knots  ( ak(j),j  —  —3,  n+3) 

An  integer  =  n+3 

A  working  array  (( a{i,j ),  »  =  1,  =  2,  4) 

A  working  vector  (diag  (j),j  =  1,  «+3) 

A  working  vector  (b(j),  j  -  1 ,  n) 


Output 

C 

ss 

IFAIL 


The  vector  of  B-spline  coefficients  (c(j),  j  -  1,  u-j-3) 

The  residual  sum  of  squares  (here,  ~0) 

The  error  indicator — 

=  0  for  a  successful  call 
=  1  for  un-ordered  knots 
=  2  for  non-positive  weights 
=  3  for  un-ordered  abscissae 

=  4  if  the  no.  of  distinct  data  abscissae  does  not  exceed  n  by  at  least  three 
=  5  if  there  are  too  many  knots  for  the  number  of  data  points 
=  6ifNCAP3  *  NCAP +3 

(The  application  here  is  such  that  IFAIL  -  0  should  be  the  only  output) 


User  Subroutine  Called 
CHECKD 

Calling  Routine 
SURF 


A4.3  ROMIN  Subroutine 
Purpose 

ROMIN  finds  the  local  minimum  of  an  unconstrained  function  of  n  variables  using  the 
method  of  Rosenbrock.10  12  In  this  application  ROMIN  is  called  by  MAXMAN  with  n  ■=  1 
(Mach  number  is  the  only  independent  variable). 

Use 

CALL  ROMIN  (N,  X,  FUNCT,  STEP,  MON1TR) 
with  input 

N  The  number  of  independent  variables,  n 

X  A  vector  giving  an  initial  estimate  of  the  solution  (x(i),  i  =  1,  n) 

FUNCT  (N,  X,  F)  An  EXTERNAL  subroutine  to  calculate  the  function /(*)  to  be 
minimised 
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STEP  An  initial  step  length  for  all  co-ordinate  directions 

MONITR  (N,  X,  F,  R,  B,  CON,  NR)  EXTERNAL  subroutine  to  monitor 

convergence  and  provide  diagnostics 

and  subroutine  arguments 

F  The  value  of  the  function  f(x)  to  be  minimised 

R  The  total  number  of  function  evaluations 

B  The  Euclidian  norm  of  the  vector  representing  progress  since  the  last  axis 

rotation 

CON  A  logical  variable  set  to  .FALSE,  by  ROM1N  initially,  and  set  to  .TRUE,  by 
by  MONITR  to  stop  the  process 
NR  A  monitor  index  supplied  by  ROM1N — 

-  0  for  initial  calculation 
=  1  for  function  evaluation 

-  2  for  new  axis  rotation 

Output  is  via  subroutine  MONITR  and  ti  c  vector  X. 

Note:  In  the  current  application  ROMIN  is  called  by  MAXMAN  as 
CALL  ROMIN  (I,  X,  SEP,  STEP,  MONSEP) 

with  STEP  =01. 

User  Subroutines  Called 
SEP  MONSEP 

System  Routine  Called 
SQRT 

Calling  Routine 
MAXMAN 


1 


A4.4  SPDER3  Subroutine 


Purpose 

SPDER3  evaluates  the  cubic  spline/(c)  and  its  first  derivative  df(x)/dx,  based  on  normalised 
B-spline  coefficients  and  associated  knot  positions. 


Use 


CALL  SPDER3  (NCAP,  AK,  C,  X,  1FA1L,  F,  DFDX,  NCAP3) 
with  input 

NCAP  The  number  of  intervals,  n 
AK  The  vector  of  knots  ( ak(j),j  =  —3,  n+3) 

C  The  B-spline  coefficients  (c(j),j  -  1,  n+3) 

X  The  value  x  at  which /(a)  is  required 

NCAP3  An  integer  =  n -f  3 
and  output 

F  The  value  of  f(x) 

DFDX  The  derivative  df{x)jdx 

IFAIL  An  error  indicator — 

=  0  if /(*),  df(x)dx  are  successfully  calculated 
=  1  if  x  is  outside  the  valid  range  ak{0)  <  x  ak(n) 
-  2  if  NCAP3  /  NCAP-f  3 


Calling  Routine 
SURF 
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A4.5  SURF  Subroutine 

Purpose 

SURF  takes  spline-fit  information  representing  a  surface  stored  in  vector  form  and  provides, 
at  the  grid  point  ( x ,  y),  the  function  value  z  =  /(x,  y),  and,  if  requested,  the  derivatives  2 >//<>x, 
iifldy  and  iy]lxty. 

The  spline  data  represent  knots  and  B-spline  coefficients  for  least-square  curve-fits  of  z 
against  x  for  discrete  values  of  y.  /(x,  y)  is  evaluated  by  bounding  the  requested  y  with  four 
appropriate  values  (>■(/),  i  =  1,4),  evaluating  the  four  values  (z(»)  =  /(x,  y(i)),  and  performing 
a  final  cubic  fit  to  these  four  points  to  obtain  z  =  /(x,  y).  The  calculation  includes  checking  x 
and  y  against  valid  ranges  allowed  hj  the  data. 

SURF  also  allows  for  the  degenerate  case  where  the  “surface"  is  a  single  curve. 


Use 


CALL  SURF  (TAB,  X,  Y,  Z,  DZDX,  DZDY,  D2ZDXY,  1ERR) 


Input 


TAB  7  he  array  of  spline  coefficients  and  knots  in  “standard"  form  (see  Section  4.6) 
X  The  first  independent  variable 

Y  The  second  independent  variable 

1ERR  An  integer  set  to  —1  if  derivatives  are  required 


Output 

Z 

DZDX 

DZDY 

D2ZDXY 

IERR 


The  function  value  z  —  /(x,  y ) 

The  first  derivative 
The  first  derivative  i »//2>y 
The  second  derivative  32//2 )x5y 
An  error  indicator  =  10M  FlFAIL,  where 
M  =  J  if  number  of  curves  in  TAB  is  2,3  or  more  than  30 
=  2  if  y  is  outside  the  valid  range 
=  3  if  number  of  intervals  for  any  curve  excees  19 
=  4  if  an  error  occurs  evaluating  (z(i),  /  -  1,4) 

=5  5  if  an  error  occurs  fitting  a  cubic  to  z(/) 

=  6  if  an  error  occurs  evaluating  z  or  Dz/i>y 
-  7  if  an  error  occurs  fitting  a  cubic  to  iz/ix 
=  8  if  an  error  occurs  evaluating  ’izj'ix  or  S2z/3x3y 
and  1FA1L  is  defined  as  for  CUBICS  ( M  —  5,  7)  or  for  SPDER3  (Af  =  4,  6,  8) 

IERR  -  0  ;ndicates  a  successful  calculation 

(In  use  with  the  present  program  !ERP„  values  of  20  or  41,  indicating  y  or  x  values  out  of 
range  respectively,  will  normally  occur  when  invalid  Mach  number  or  altitude  requests  are  made.) 

User  Routines  Called 

SPDER3  CUBICS 


Callint;  Routines 
THRUST 


AERO  (if  needed) 
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APPENDIX  5 

Sample  AIRCRAFT  Program  Library  Listing 


A5.1  Main  Program 


C  PROGRAM  AIRCRAFT  EXANPLE 

C  G.  KIPP  ARL  FEB  1976 
C 

CALL  P1IN 
CALL  TABLE 
STOP 
END 

BLOCK  DATA 
C 

C  COMMON  DATA  FOR  PROGRAMS  PI,  P2  AND  P3 
C 

COMMON  /C/C(200) 

DIMENSION  H ICAO < 4 ) ,  TICAO(4),  HARDU < 6 > ,  TARDU ( A ) ,  AMET ( 3 ) , 

+  N 1 0 < 10) ,  PT AB < 3 ) ,  GRID1 (2) ,  GRID2 < 6 ) 

EQUIVALENCE  (  HICA0,C<  1)),(  TICA0,C(  5)>,(  HARDU, C(  9)), 

+(  TARDU,C(  15)), <  NI0,C<  21)), (  PTAB,C<  31)), (  AMET ,C(  34)), 

♦  <  GAMMA, C(  37)), <  R,C<  38)), (  G,C(  39) ), (RADIAN, C(  40)), 

♦  <  HOUR , C (  41)), (  GRID1 ,C< 121 ) ) , 

+  <  G R I D 2 , C ( 1 23) ) 

DATA  HICAO  /O.O,  36089.24,  65616.8,  104986.88/ 

DATA  TICAO  /288.15,  216.65,  216.65,  228.65/ 

DATA  HARDU  /O.O,  25000.,  45000.,  54000.,  70000.,  104986.88/ 
DATA  TARDU  /301 . ,  253.,  205.,  193.,  2137,  238.912/ 

DATA  AMET  /"ICAO  ARDU  V 

DATA  GAMMA,  R,  G,  RADIAN,  HOUR 
+  /I. 4,  3089.78,  32.17405,  57.2957795,  3600./ 

DATA  NIO  /I,  2,  3,  4,  5,  6,  7,  8,  9,  10/ 

DATA  PTAB  /2116.22,  472.680,  114.345/ 

DATA  GRID1  /"ES" ,"HP"/ 

DATA  GRID2  /"ENERGY  STATE",'  ALTITUDE'/ 

END 


non  o  o  n  o 
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A5.2  Sample  THRUST  Subroutine 


SUBROUTINE  THRUST 
C 

C  THRUST  ROUTINE  FINDS  NET  THRUST  AND  FUEL  FLOW  FOR  COBRA  AIRCRAFT, 

C  BASED  ON  STANDARDISED  DATA  FORMAT  FOR  ARRAY  F 
C  REV  GUK  FEB  77. 

C 

COMMON  /B7  B < 2 00 ) 

COMMON  /FBATA/  F(2000) 

EQUIVALENCE  (  XM,8<  3)),(  HP,B(  6)),(  PLA  ,B (  15)), 

+  (  DTEMP,B<  19)), (  MT,B(  19)), (  FN,B<  41)), (  UF,B(  42)), 

+  (  XMLIM,B(  61)), (  IERR,B<  75)) 

C 

C  IDENTIFY  POWER  LEVEL  (MIL:  1=1,  MAX:  1=2) 

C 

DO  50  1=1,2 

AI  =  FLOAT ( I ) 

IF  <ABS(PLA-100.*A1)  .LE.  0.001)  GOTO  100 

50  CONTINUE 
70  TYPE  75,  PLA 

75  FORMAT  <//  CANNOT  PROCESS  ',F6.2/X  THRUST  REQUEST  YET  ) 

STOP 

100  NJ  =  1 

IF  (I.EQ.2)  GOTO  200 

C 

C  IF  MAX  THRUST.  SKIP  TO  CALCULATE;  OTHERWISE  JUMP  TO  MIL  DATA 
C 

DO  150  1=1,2 

L  =  F ( N  J )  +  0.001 
150  NJ  =  L  +  NJ 
200  NJ  =  NJ  +  1 

SET  NOMINAL  MACH  NO  LIMIT 

*****  N.B.  SUPPLY  NEXT  TWO  LINES  FOR  ALL  AIRCRAFT 

LX  =  F (N J+1 )  +  0.001 
XML1H  =  F (NJ+LX+3) 

PROCESS  SPLINE  DATA  FOR  FN  (1=1)  AND  UF  (1=2) 

CO  250  1=1,2 

CALL  SURF  (F(NJ),  XM,  HP,  D1,  DFDH,  DFDH,  D2FDHh,  IERR) 
IF  ( IERR.NE.O)  GOTO  300 
L  =  F(NJ-I)  +  0.001 
NJ  =  L  ♦  NJ 
GOTO  (215,230),  I 
215  FN  =  D1 

GOTO  250 

230  WF  =  D1 

250  CONTINUE 
300  CONTINUE 
RETURN 
END 


1  .  ASJ  Sample  Aerodynamics  Routines 


SUDRUUUNE  AERO 
C 

C  AERO  POGIINE  FINPS  DRAG,  I,  ANB  ANGLE  OF  ATTACK,  ALPHA  (DEGREES) 
C  liMX  ARL  MAR  74 
C 

COMMON  /■/  1(200) 

COMMON  /C/  C(20V) 

COMMON  /EDATA/  E(2000) 

DIMENSION  XNDCBS(20),  BELCDSI20) 

EQUIVALENCE  (RADIAN, C (  40) > 


♦  i 

EQUIVALENCE  ( 
FN,»<  419), 

6N,I( 

2) ) ,  ( 

XN,I(  3) ) ,  ( 

U,t( 

17)), 

♦  ( 

♦  ( 

ALPHA, l(  43)), ( 
CLALF,t(  47)), ( 

DIM, 11 
1ERR,B( 

44)), ( 
73)), 

D,B(  45)), ( 

CLT ,D( 

441), 

♦  ( 

NDCBS.K  80) ) ,  (XH0CDS,B< 

01 ) ) ,  (DELCDS.Dl 1 01 ) ) 

FINTRP(X,XI,VI,X2,Y2)  ■  Tl  ♦  (X-XI)  »  (T2-Y))  /  (X2-XI) 

C  • 

C  ITERATE  FOR  ALPHA  AND  CLT .  ASSUME  INITIAL  ALPHA  IS  GIVEN. 

C 

DALF  *  O.S 
NIT  =  0 
100  CONTINUE 

CALL  TRIHCL 

IF  UERR.NE.O)  GOTO  900 

IF  (AIS(  DALF )  .LT.O.Ot )  GOTO  MO 

DALF  *  (FN*S1ND( ALPHA)  ♦  CLY*DIN  -  OH*U)  / 

♦  (FN»C0S0<ALPHA>/RAD1AN  ♦  CL ALF ♦DIR) 

ALPHA  =  ALPHA  -  DALF 
NIT  *  NIT  ♦  I 
IF  (N1T.LE.20)  GOTO  tOO 
IERR  «  102 
GOTO  900 
(10  CONTINUE 
C 

C  NOU  FIND  DRAG 
C 

C  1=1  FINDS  CLMAX  AND  CHECKS  CLT 
C  1=2  FINDS  CDMIN 
C  1=3  FINDS  CDL 
C 

NJ  =  2 

DO  400  1=1,3 

GOTO  (120,(20,140)  I 
120  X  *  XH 

Y  =  0. 

GOTO  140 

C 

C  IF  IKMHIN,  USE  NHIN  CURVE 
C 

140  X  •  CLT 

Y  *  AHAX1  (E(NJ»2),XM) 

140  CALL  SURF  (E(HJ),  X,  Y,  Z,  01,  92,  03,  IERR) 

IF  (IERR.NE.0)  GOTO  900 
L  *  E(NJ-l)  ♦  0.001 
HJ  *  NJ  ♦  L 
GOTO  (200,220,240)  I 
200  CLMAX  >  Z 

IF  (CLT. LE. CLMAX)  GOTO  400 
IERR  >  103 
GOTO  900 

220  CDMIN  •  Z 

60TG  400 

240  CDL  *  Z 

400  CONTINUE 

C 

C  FIND  STORE  DRAG 
CDS  *  0. 

IF  (NDCDS.EO.O)  GOTO  550 
C 

IF  (XH.LE.XHDCDSt I ) )  CDS  «  DELCDS(I) 

IF  ( XN.GE.XHDCDS(NICDS) )  CDS  •  DELCDS(NDCDS) 

IF  (XR.LE.XMBCDSt 1 )  .OR.  XM.GE.XNDCDS(NDCDS) >  GOTO  540 
DO  430  1«2,NBCBS 

430  IF  (XM.LE.XHDCDS(D)  GOTO  300 
IERR  >  104 
GOTO  900 
300  IN  •  1-1 

CDS  >  F1NTRP  (XH,  XMDCDS(IN),  BELCBSf IN),  XHDCISU),  DELCDS(D) 
340  CDS  >  CDS  •  1.0E-04 
330  A  =  (CDNIN  ♦  CDL  *  CDS)  »  DIN 
900  CONTINUE 
RETUNN 
END 


ooo  ono  o  o  o  Anonoo 


SUBROUTINE  TRIHCL 


ROUTINE  FOR  CLTRIN  AND  DCL/DALPHA ,  GIVEN  ALPHA  IN  DEGREES,  USING 
LINEAR  APPROXIMATIONS  TO  CLTRIN  VS  ALPHA.  (C  ONLY) 

G'JK  NAR  76 

CONNON  /»/  B <  200 ) 

DIMENSION  CLTAB<4,4>,  AAT<4,4),  NT(4),  XMT<4),  CLDUM(4),  CLAPUNt 

I 

EQUIVALENCE  (  XN,B(  3)),<  ALPHA, D<  43)),  <  CLT,B(  46)), 

♦  <  CLALF,B<  47)), «  IERR,B(  75)) 

FINTRP  ( X , X 1 1 YT , X 2 , Y 2 )  =  Y!  +  (X-Xl)  *  (Y2-Y1)  /  <X2-XI ) 

DATA  CLTAB  /. 05, 1 .47, I .80,0. , .02, .67, 1 .42, 1 .52,0. , .97,3*0. , 

+  .39,2*0./ 

DATA  AAT  /O. , 15.2,30.0,2*0. ,6.6,23.0,27.0,0.2, 12.6,2*0. , 1 .0, 

+  9. 2, 2*0./ 

DATA  NT  /3,4,2,2/ 

DATA  XNT  /. 2, .8, 1.2, 2.0/ 

DATA  NCURV  /4/ 

CHECK  ALPHA  LIMITS 

XMTR  =  XN 

IF  (XH.LT.XNTII))  XNTR  =  XHT(1) 

DO  5C  11=2, NCURV 
I  =  II 

50  IF  ( XM.LT .XNT ( I )  >  GOTO  100 
100  IM  =  I-l 

ALFLIH  =  FINTRP  (XNTR,  XNT(IN),  AAT (NT( IH) , IH) ,  XMT(I), 

♦  AAT (NT( I ) , I } ) 

IF  (hLPHA.GT. ALFLIH)  GOTO  900 

USE  LINEAR  INTERPOLATION  IN  2-DINENSIONS  TO  GET  CL,  CLA 

DO  300  K=IN, I 

DO  200  J=2,NT(K) 

200  IF  (ALPHA. LE.AAT(J,K))  GOTO  300 

300  CALL  IKTRP  (CLTABM ,K),  AAT ( 1 ,K) ,  CLADUN(K),  ALPHA,  CLDUN(K),  J) 
CLT  =  FINTRP  (XMTR,  XMT(IM),  CLDUN(IM),  XMT(I),  CLDUM(I)) 

CLALF  =  FINTRP  (XMTR,  XMT(IM),  CLADUMi IN) ,  XMT(I),  CLADUM(I)) 
GOTO  950 

EITHER  MACH  NO  OR  ALPHA  EXCEEDS  DATA  LIMITS 

900  IERR  =  101 
950  RETURN 
END 
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APPENDIX  6 

Stonge  Allocation  of  Labelled  Common  Arena 

The  following  paragraphs  define  storage  allocation  in  the  various  labelled  common  areas  in 
programs  AIRCRAFT,  P2  and  P4. 

In  general,  an  attempt  has  been  made  to  restrict  data  communication  to  a  COMMON 
area  B  of  dimension  200,  and  to  use  this  area  for  all  three  programs ;  however,  this  was  not  always 
possible,  and  other  COMMON  areas  are  used. 

In  particular,  COMMON  C  has  been  used  for  data  constants  set  by  the  BLOCK  DATA 
subroutine,  and  areas  TABLES  and  TABLET  (AIRCRAFT),  and  ESLIST,  RPTLST  and 
TABLES  (P2)  have  been  used  to  store  grid  data  prior  to  calling  output  routines. 

COMMON  areas  CLIMBS  and  ZEROPS  are  used  with  the  optimised  grids  when  producing 
climb  profiles  and  optimum  sustained  turn  rate  data.  Areas  EDATA  and  FDATA,  both  of 
dimension  2000,  are  reserved  for  aerodynamic  and  propulsion  data,  and  may  be  enlarged  or 
reduced  as  required.  COMMON  area  E  is  used  by  programs  P2  and  P4  in  transferring  aircraft 
identification  data  front  input  to  output  files. 

The  units  of  variables  in  program  AIRCRAFT  are  those  associated  with  interna!  storage; 
in  general  Imperial  units  are  used  with  slight  variations  to  allow  for  standard  aeronautical 
practice.  Output  units  of  this  program  are  chosen  by  the  user.  The  units  of  variables  in  programs 
P2  and  P4  are  Imperial  or  SI,  depending  on  the  output  of  program  AIRCRAFT. 
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A6.1  Program  AIRCRAFT 
A6.1.I  Labelled  Common  B 


Word 

Variable 

Description 

1 

..  ES 

Energy  height  (ft) 

2 

..  GN 

. .  Normal  load  factor  (g) 

3 

..  XM 

M  ich  number 

4 

..  T 

Ambient  temperature  (K) 

5 

. .  TS 

. .  iCAO  standard  temperuture  (K) 

6 

. .  HP 

. .  Pressure  altitude 

7 

..  P 

. .  Ambient  pressure  (lb/ft2) 

8 

..  A 

. .  Speed  of  sound  (ft/s) 

9 

..  V 

True  airspeed  (ft/s; 

10 

..  H 

Geopotential  height  (ft) 

11 

. .  1  UN  ITS  . 

Unit  flag 

12-13  .. 

. .  — 

. .  Unallocated 

14 

..  SV/P 

. .  Wing  sweep  angle  (deg) 

15 

. .  PLA 

Throttle  setting  code 

16 

..  SREF 

Wing  reference  area  (ft2) 

17 

..  W . 

. .  Combat  weight  (lb) 

18 

..  DTEMP  . . 

. .  Deviation  from  reference  atmosphere  (K) 

19 

. .  MT 

Atmosphere  type 

20-35  . . 

. .  — 

Local  working  storage 

36 

. .  ALAPSE 

. .  Ambient  temperature  lapse  rate  (K/ft) 

37 

..  LEV 

ICAO  atmosphere  layer  number 

38 

. .  ALPSTD 

. .  ICAO  atmosphere  temperature  lapse  rate  (K/ft) 

39 

..  Cl 

. .  Constant  -  yRiW2/2g  (ft/K) 

40 

. .  C2 

Constant  =  R/g  (ft/K) 

41 

..  FN 

. .  Net  thrust  (lbf) 

42 

..  WF 

Fuel  flowrate  (Ib/hr) 

43 

..  ALPHA.. 

. .  Trimmed  angle  of  attack  (deg) 

44 

..  DIM 

. .  Dimensionalising  force  =  O-SyPM^S  (lbf) 

45 

..  D 

.  Drag  (lb) 

46 

.  •  CLT 

. .  Trimmed  lift  coefficient 

47 

..  CLALF  .. 

. .  Lift  curve  slope  =  dCiJda.  (deg  ') 

48 

. .  — 

..  Unallocated 

49 

. .  IGN 

. .  Load  factor  loop  index 

50 

. .  CDT 

. .  Trimmed  drag  coefficient 

51 

..  DELES  .. 

. .  Height  variable  increment  (ft) 

52. 

. .  NES 

. .  Number  of  points  on  height  grid 

53 

. .  DELXM  . . 

. .  Mach  number  increment 

54 

. .  NXM 

. .  Number  of  points  on  Mach  number  grid 

55 

. .  DELGN 

. .  Load  factor  increment  (g) 

56 

. .  NGN 

. .  Number  of  points  on  load  factot  grid 

57 

. .  ES0 

. .  Initial  height 

58 

XM0 

. ,  Initial  Mach  number 

59 

..  GN0 

. .  Initial  load  factor 

60 

. .  10PT  . . 

. .  Optimised  grid  switch 

61 

..  XMLIM.. 

..  Upper  Mach  number  limit  on  data 

62-63  . . 

. .  FRAME 

Aircraft  name  (double  precision  text) 

64 

. .  IES 

. .  Height  loop  index 

65 

..  IXM 

Mach  number  loop  index 

66-67  .. 

. .  DAY(2)  . . 

. .  Date 

68-69  . . 

•  •  CLOCK(2) 

Time  of  day 

70-71  .. 

. .  — 

Unallocated 

72 

..  SWPM1N 

Minimum  wing  sweep  angle  (deg) 

Word 

Variable 

73 

. .  SWPMAX 

74 

.  ,  - - 

75 

IBRR  .. 

76 

..  10UT  .. 

77 

..  IGRID  .. 

78-79  , 

. .  — 

SO  .. 

. .  NDCDS 

81-100 

. .  XMDCDS(20)  . 

101-120 

. .  DEL.CDS(20)  . . 

121-134 

.  .  ROLE(14) 

135  .. 

.  IPSTYP 

136  . 

. .  CGPERC 

137 

..AM 

138  .. 

. .  ALB 

139  .. 

..  AWF  .. 

140-200 

A6.1.2  Labelled  Comma  C 

Word 

Variable 

1-4  .. 

. .  HICAO(4) 

5-8  .. 

. .  TICA014) 

9-14  .. 

. .  HARDU(6) 

15-20  .. 

. .  T4RDU(6) 

21-30  .. 

. .  N  10(10)  . . 

31-33  ... 

. .  PTAB(3) 

34-36  . . 

. .  AM£T(3) 

37 

. .  GAMMA 

38 

..  R . 

39 

. .  G . 

40 

..  RADIAN 

41 

..  HOUR  .. 

42-120 

— 

121-122 

..  GRIDK2) 

123-128 

. .  GR1S2(6) 

129-200 

. .  — 

46.1.3  Labelled  Common  CLIMBS 

Word 

Variable 

1 

.  TI 

2 

.  FU 

3 

.  RA 

A6.1.4  Labelled  Common  EDATA 

Word 

Variable 

1-2000 

.  E(2000)  .. 

A6.I.5  Labelled  Common  FPATA 

Word 

Variable 

1-2000 

.  F(2000)  .. 

Description 

Maximum  wing  sweep  angle  (deg) 

Unallocated 

Error  Flag 

Output  file  switch 

Height  grkt  switch 

Unallocated 

Numbct  of  points  in  drag  count  table 
Mach  number  list  for  drag  count  tabic 
Drag  count  table 
Description  of  current  role 
F.nei  gy  parameter  switch 
e  g.  position  (%MAC) 

Conversion  factor  (metres  per  foot) 
Conversion  factor  (pounds  per  kilogram) 
Conversion  factor  (kg/sec  per  lb/hv) 
Unallocated 


Description 


Description 

Storage  reserved  for  propulsion  data 


ICAO  atmosphere  pressure  altitude  table  (ft) 
ICAO  atmosphere  temperature  table  (K) 
ARDU  atmosphete  pressure  altitude  table  (ft) 
ARDU  atmosphere  temperature  table  (K) 
Integer  constants,  1-10 
Base  pressures  of  ICAO  temperature  profile 
(lbf/ft2) 

Atmosphere  text 
Specific  heat  ratio  (1  -4) 

Gas  constant  for  air  (3089  -78  ft2/Ksz) 
Acceleration  due  to  gravity  (32- 17405  ft/sa) 
Degrees/radian  (57  •  2957795  deg/rad) 
Seconds/hour  (3600  s/hr) 

Unallocated 

Height  variable  input  text 
Height  variable  output  text 
Unallocated 


Description 

Time  to  climb  (s) 
Fuel  used  (lb) 
Range  (n.m.) 


Description 

Storage  reserved  for  aerodynamic  data 


f  l'T 


A6.1.6  Labelled  ComwoM  TABLES 

#W 

Variable 

Description 

1-50 

. .  PSTAIX50) 

,  , 

. .  Energy  parameter  output  table 

51-100 

. .  OMTAB(SO) 

.  . 

. .  Turn  rate  output  table  (deg/s) 

101-150 

..  GNTAB(50) 

. .  Load  factor  grid  (g) 

A6 1.7  Labelled  Common  TABLET 

Word 

Variable 

Description 

1-50 

. .  XMTAB(50) 

,  , 

. .  Mach  number  output  table 

51100 

. .  VTAB(50) 

. .  True  airspeed  output  table  (kn) 

10t-150 

. .  HPTAB(50) 

,  , 

. .  Pressure  height  output  table  (ft) 

151-200 

. .  WFTAB(50) 

. .  Fuel  flowrot?  output  table  (Ib/hr) 

A6.1.8  labelled  Common  ZEROFS 

Word 

Variable 

Description 

1 

..  GNPS0  .. 

. .  Optimum  sustained  load  factor  ( g ) 

2 

..  OMPS0  .. 

. .  Optimum  sustained  turn  rate  (deg/s) 

> 

%(  ,  , 

. .  XMPX0  . . 

. .  Mach  number  at  optimum  condition 

4 

. .  VPS0 

. .  True  airspeed  at  optimum  condition  (kn) 

5 

..  HPPS0  .. 

. .  Pressure  altitude  at  optimum  condition  (ft) 

6 

. .  WFPS0  . . 

Fuel  flowrate  at  optimum  condition  (slug/s) 

A6.2  Program  P2 

A6.2.1  Labelled  Common  B 

Word 

Variable 

Description 

1 

. .  ES 

. .  Energy  height  (ft  or  m) 

2 

..  GN 

. .  Load  factor  (g) 

3 

.  XM 

. .  Mach  number 

4-10  .. 

. .  — 

. .  Not  used 

11 

. .  IUMITS  . . 

. .  Unit  f>ag 

12-13  .. 

. .  — 

. .  Not  used 

14 

..  SWP 

. .  Wing  sweep  (deg) 

15 

. .  PLA 

. .  Throttle  setting  (%  max) 

16 

, .  — 

. .  Not  used 

17 

. .  W 

Combat  weight  (lb  or  kg) 

18 

. .  DTEMP . . 

. .  Deviation  from  reference  atmosphere  (K.) 

19-24  , . 

. .  — 

. .  Not  used 

25 

. .  AOPT 

. .  Terminal  reply  (local  to  P2IN) 

26 

. .  ICREAT. . 

. .  Flag  to  indicate  on-line  data  generation 

27  .. 

. .  IONDSK 

. .  Flag  to  indicate  both  differential  MMD  files  on 

disk 

28-47  . . 

. .  — 

. .  Not  used 

48 

..  MSECS0. . 

. .  CPU  time  used  (ms) 

49-50  . . 

. .  — 

. .  Not  used 

51 

. .  DELES  . . 

. .  Height  variable  increment  (ft  or  m) 

52 

. .  NES 

. .  Number  of  points  on  height  grid 

/ DELXM . . 

Mach  number  increment 

JJ  *  * 

\  DELPS  .. 

. .  Energy  rate  increment  (MMD) 

<A 

/  NXM 

Number  of  points  on  Mach  number  *;rid 

JH  .  . 

• '  \NPS 

. .  Number  of  points  on  energy  late  grid  (MMD) 

Description 

. .  Load  factor  increment  (g) 

. .  Number  of  points  on  load  factor  grid 
. .  Initial  height  (ft  or  m) 

Initial  Mach  number 
. .  Initial  energy  rate  (MMD) 

. .  Initial  load  factor 
. .  Not  used 

. .  Aircraft  name  (double  precision) 

. .  Not  used 
. .  Date 
. .  Time  of  day 
. .  Not  used 

. ,  Description  of  current  role 
. .  Energy  parameter  flag 
. .  Not  used 

Description 
. .  Not  used 

, .  Units  text  (double  precision) 

Energy  variable  text  (double  precision) 
. .  Energy  unit  text  (double  precision) 

. .  Energy  variable  scaling  factor  text 
. .  Length  unit  text 
. .  Weight  unit  text 
Not  Used 


Description 
.  Not  used 

Maximum  value  on  height  grid  (ft  or  m) 
Maximum  value  on  Mach  number  grid 
Maximum  value  on  load  factor  grid 
, .  Power  regime 
Not  used 

. .  Atmosphere  description 
Output  filename 

. .  Output  filename  (double  precision) 

. .  Not  used 

Description 
. .  Energy  state  grid 
. .  Turn  rate  grid  (MMD) 

Description 

, .  Table  of  load  factor  request  replies 
Description 

. ,  Mach  number  grid 

Energy  rate  grid  (MMD) 

. .  Load  factor  grid 


A6.3  Program  P4 


A6.3.1  Labelled  Common  B 


Word 

Variable 

Description 

1-10 

_ 

. .  Not  used 

11  .. 

. .  IUNITS 

..  Unit  flag 

12,  13 

. .  Not  used 

14 

..  SWP 

. .  Wing  sweep  (deg) 

15  .. 

. .  PLA 

. .  Throttle  setting  (%  max) 

16 

. .  — 

. .  Not  used 

17  .. 

..  W . 

. .  Combat  weight  (lb  or  kg) 

18  .. 

. .  DTEMP 

. .  Deviation  from  reference  atmosphere  (K) 

19-24 

. ,  — 

. .  Not  used 

25  .. 

..  YLMIN  , . 

. .  Lowest  energy  rate  for  zero  turn  rate  boundary 

26-27 

. .  NAM(2) . . 

. .  Input  filename 

28-29 

. .  — 

. .  Not  used 

30  .. 

..  1DATA 

. .  Data  set  counter 

31  .. 

. .  NDATA 

. .  Number  of  data  sets 

32  .. 

. .  IC 

. .  Running  count  of  points  on  zero  turn  rate  contour 

33-34 

. .  1C0UNT(2) 

. .  Number  of  points  on  two  zero  turn  rate  contours 

35 

. .  NSAV 

. .  Default  pen  command  if  above  zero  turn  rate 
boundaries 

36-50 

. .  — 

Not  used 

51  .. 

..  DELES  .. 

. .  y-axis  increment  (ft  or  m) 

52  .. 

. .  NES 

. .  Number  of  points  on  >>-axis  grid 

53 

..  DELXM . . 

. .  x-axis  increment 

54  .. 

. .  NXM 

. .  Number  of  points  on  x-axis  grid 

55  .. 

. .  DELGN  . . 

. .  Load  factor  increment 

56  .. 

..  NGN  .. 

Number  of  points  on  load  factor  grid 

57  .. 

. .  ES0 

Initial  j'-axis  value 

58  .. 

. .  XM0 

. .  Initial  x-axis  value 

59  .. 

..  GN0 

. .  Initial  load  factor 

60-61 

. .  — 

Not  used 

62-63 

. .  FRAME . . 

. .  Aircraft  name  (double  precision) 

64-65 

. .  — 

. .  Not  used 

66-67 

..  DAY(2)  .. 

. .  Date 

68-69 

. .  CLOCK(2) 

. .  Time  of  day 

70-120 

. .  — 

. .  Not  used 

121-134 

..  ROLE(14) 

. .  Description  of  current  role 

135 

. .  IPSTYP  . . 

. .  Energy  parameter  flag 

136-200 

. .  — 

. .  Not  used 

A6.3.2  Labelled  Common  C 

Word 

Variable 

Description 

1-43  .. 

.  ,  - 

. .  Not  used 

44-45  .. 

. .  ENERGY 

. .  Energy  variable  text  (double  precision) 

46-47  .. 

. .  PSUN1T 

. .  Energy  unit  text  (double  precision) 

48 

. .  FACT  . . 

Energy  variable  scaling  factor  text 

49 

. .  ALUNIT 

. .  Length  unit  text 

50-200 

. .  — 

Not  used 

A6.3.3  Labelled  Common  E 


Word 

Variable 

Description 

1-9  . 

f  ' — 

Not  used 

10 

. .  RATING 

Power  regime 

11-14  .. 

.  t  — 

Not  used 

IS 

..  AMET . 

Atmosphere  description 

16-2C0 

. .  — 

Not  us«l 

A6.3.4  Labelled  Common  CONPLT 

Word 

Variable 

Description 

1 

..  CONLAB 

Flag  to  request  contour  identification 

2 

..  XSCALE 

x-axis  scale  in  units  per  inch  of  plot 

3 

. .  YSCALE 

y-axis  scale  in  units  per  inch  of  plot 

4 

. .  SMTH . 

Flag  to  request  smoothing  of  contours 

5 

. .  I  DOT . 

Integer  indicating  contour  texture 

A6.35  Labelled  Common  GDMDOT 

Word 

Variable 

Description 

1 

. .  DINT . 

Contour  smoothing  interval  (in.) 

A6.3.6  Labelled  Common  GDMSWT 

Word 

Variable 

Description 

1 

. .  DGS  . 

Flag  to  bypass  biquadratic  diagonal  interpolation 

2 

. .  NLEV . 

Levei  of  contour  being  plotted 

A6.3.7  Labelled  Common  LIMS 

Word 

Variable 

Description 

1-400 

. .  XTAB{400) 

x -co-ordinates  of  differential  MMD  boundaries 

401-800 

..  YTAB(400) 

^-co-ordinates  of  differential  MMD  boundaries 

A6.3.8  Labelled  Common  NLEV0 
Word  Variable  Description 

1  . .  . .  NLEV0 . I-evcl  number  of  zero  contour 

2  . .  . .  IDOT0 . Integer  giving  texture  of  zero  contour 

3  . .  IDEF . Integer  giving  texture  of  nonzero  contours 
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APPENDIX  7 
P2  Program  Library 

A7.I  P2  Program 
Purpose 

The  program  P2  processes  data  files  produced  by  program  AIRCRAFT  and  produces 
plotter  files  and  output  files  for  input  to  program  P4. 

Plotter  files  (named  P2.PLT)  produced  are  of  two  types:  (a)  energy  parameter  versus  turn 
rate  for  a  specified  altitude,  and  (6)  turn  rate  and  energy  parameter  versus  Mach  number  for 
various  load  factors,  for  a  specified  height  (energy  height  or  pressure  height). 

P4  input  files  produced  are  P2.CON,  P2A.CON  and  P2DIFF.CON  (energy  parameter 
contours  on  a  Mach  number  versus  pressure  altitude  grid)  and  P2,OFT,  P2A.OPT  and 
P2DIFF.OP7  (optimum  turn  rate  contours  on  an  energy  state  versus  energy  parameter  grid). 

In  addition  a  conversational  routine  is  included  which  assists  in  calculating  Maximum 
Manoeuvre  Persistence. 

The  main  program  calls  routine  P21N  to  determine  parameters  defining  the  type  of  run, 
and  then  calls  a  control  subroutine  for  that  run,  after  which  execution  ceases. 

Use 

RUN  P2  (Monitor  command) 

User  Routines  Called 

P2IN  PSCON  RATE1  RATE2  PSDIFF  R2DIFF  MMP  GRID 


A7.2  GRID  Subroutine 

Purpose 

GRID  reads  from  disk  storage  data  files  produced  by  program  AIRCRAFT,  expressed  on  a 
(Mach  number,  load  factor)  grid  for  each  height  requested,  and  prepares  a  disk  file  P2.PLT  for 
plotting.  Turn  rate  and  energy  parameter  are  plotted  against  Mach  number  for  each  load  factor 
requested.  A  separate  page  of  plots  is  produced  for  each  value  of  the  height  variable  (energy  height 
or  pressure  height). 

COMMON  areas  B,  C  and  TABLES  provide  variable  communication. 


Use 

CALL  GRID  (WORK) 

where  WORK  is  an  8000-word  vector  of  working  storage. 

Input 

COMMON  B:  NES  NXM  NGN  1UN1TS  IPSTYP 
COMMON  C:  ENERGY  PSUNIT  FACT  ALUNIT 
COMMON  TABLES:  GNLIST 

Output 

COMMON  3.  ES  GN 

[for  Routines  Called 

PLTLAB  AXIS  LINE 
System  Routines  Called 
PLOT  TEXT 
Calling  Routine 
P2  program 
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A7.3  INLAB  Subroutine 

Purpose 

INLAB  reads  an  alphanumeric  identifying  header  from  an  input  disk  file  produced  as  output 
by  program  AIRCRAFT. 

COMMON  B  and  COMMON  E  provide  data  communications. 


Use 

CALL  INLAB  (LU) 

Output 

LU  The  logical  unit  number  for  the  disk  file 
COMMON  Output 

COMMON  B:  SWP  PLA  W  DTEMP  FRAME  DAY  CLOCK 

ROLE 

COMMON  E.  RATING  A  MET 

Calling  Routine 
PSDIFF 


A7.4  INMMD  Subroutine 

Purpose 

INMMD  performs  a  conversational  dialogue  with  the  user  to  accept  text,  grid  data  and 
configuration  parameters  for  on-line  provision  of  maximum  manoeuvre  diagram  (MMD)  data. 
COMMON  B  and  COMMON  E  provide  data  communication. 

ILc 

Call  INMMD 

COMMON  Output 

COMMON  B:  SWP  PLA  W  DTEMP  FRAME  DAY  CLOCK 

ROLE  DELES  NES  ES0  DELPS  NPS  PS0 
ROLE 

COMMON  E:  RATING  A  MET 

User  Routines  Colied 

BADINP  UNITS 

System  Routines  Called 
DATE  TIME 

Calling  Routine 

INMMD  R2DIFF 


A7.5  MMP  Subroutine 
Purpose 

MMP  is  a  conversational  routine  to  assist  in  calculating  data  for  the  Maximum  Manoeuvre 
Persistance  (MMP)  Diagram.  Fuel  available  for  optimum  manoeuvres  at  discrete  range  intervals, 
together  with  energy  state,  maximum  sustained  turn  rate  and  fuel  flow,  are  requested  as  terminal 
input,  and  the  number  of  sustained  turns  at  each  range  is  provided  as  terminal  output. 

Use 

CALL  MMP 

User  Routine  Called 
BADINP 

Catling  Routine 
P2  program 
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A7.6  FLTLAB  Subroutine 

Purpose 

PLTLAB  outputs  the  identifying  header  read  by  INLAB  to  a  plotter  file.  A  header  for  a 
comparison  aircraft  is  output  with  a  second  call  to  PLTLAB  with  input  argument  N  =-  2. 
Data  communication  is  via  COMMON  areas  B,  C  and  E. 

Use 

CALL  PLTLAB(LU.N) 

Input 

LU  The  logical  unit  number  of  the  plotter  file 
N  The  aircraft  number 

N  -  1  base  aircraft 
N  =  2  comparison  aircraft 
COMMON  Input 

COMMON  B:  SWP  PLA  W  DTEMP  FRAME  DAY  CLOCK 

ROLE 

COMMON  C:  AM  UNIT 

COMMON  E:  RATING  AMET 

Systems  Routine  Called 

WHERE  PLOT  TEXT 

Calling  Routines 

RATE1  GRID 

A7.7  PSCON  Subroutine 

Purpose 

PSCON  reads  from  disk  storage  data  files  produced  by  program  AIRCRAFT,  expressed  on 
a  (load  factor,  Mach  number)  grid  for  each  pressure  height  requested,  and  prepares  a  disk  file 
P2.CON  for  input  to  program  P4.  Data  is  ovitput  on  a  (Mach  number,  altitude)  grid  for  each 
load  factor. 

Variable  coramunicttion  is  via  COMMON  areas  B,  E,  TABLES,  ESLIST  and  RPTLST. 


Use 

CALL  PSCON)  WORK) 

where  WORK  is  an  8000-work  vector  of  working  storage. 

Input 

COMMON  B:  DELES  NES  DELXM  NXM  DELGN  NGN 

ES0  XM0  GN0 
COMMON  E:  NAMOUT 

COMMON  TABLES:  GNTAB 

COMMON  RPTLST:  GNLST 

Output 

COMMON  TABLES:  GNTAB 
COMMON  ESLIST:  ESTAB 
COMMON  RFTLST:  GNLST 

User  Routine  Collet 
WRLAB 

Calling  Routines 

PSDIFF  P2  program 

A7.8  PSDIFF  Suf routine 
Purpose 

PSDIFF  performs  the  functions  of  PSCON  with  two  separate  disk  files  representing  com¬ 
parison  aircraft  and  piepares  output  files  P2.CON  and  P2A.CON  for  input  to  program  P4. 


While  doing  this  the  grid  specifications  are  checked  for  equality,  and  a  third,  differential  file, 
P2D1FF.CON  is  produced  for  input  to  program  5*4.  An  error  halt  occurs  if  the  data  grids  differ 
in  any  way. 

Data  communication  is  via  COMMON  areas  B,  E  and  ESL1ST. 


Use 

CALL  PSDIFF(WORK) 

where  WORK  is  an  8000-word  vector  of  working  storage. 

Input 

COMMON  B:  DELES  NES  DELXM  NXM  DELGN  NGN 

ES0  XM0  GN0  1PSTYP  1UNITS 
COMMON  ESLIST:  ESTAB 

Output 

COMMON  B:  SWP  PLA  W'  DTEMP  FRAME  DAY  CLOCK 

ROLE  IPSTYP  1UN1TS 
COMMON  E:  RATING  AMET  NAMOUT 

COMMON  ESLIST.  ESTAB 

User  Routine  Called 

PSCON  INLAB  WRLAB 

Calling  Routine 
P2  program 


A7.9  P2IN  Subroutine 

Purpose 

P21N  is  called  by  the  main  program  to  process  the  aircraft  identification  header  on  the 
input  disk  file  and  requests  an  option  code  for  the  type  of  run.  An  error  halt  is  forced  if  the  8000 
words  allocated  for  vector  WORK  is  too  small  to  accommodate  the  data  required  by  the  input 
grid.  It  is  the  user’s  responsibility  to  ensure  that  the  option  requested  is  consistent  with  the  pur¬ 
pose  for  which  program  AIRCRAFT  produced  the  data. 

Data  communication  is  via  COMMON  areas  B,  E  and  TABLES. 


Use 

CALL  P2IN(!OPT) 

Output 

10PT  The  option  control  as  follows: 

=  1  Energy  parameter  contour  plot 

=  2  Plot  of  energy  parameters  versus  turn  rates 

•=  3  Maximum  Manoeuvre  Diagram  (MMD)  contour  plot 

=  4  Differential  energy  parameter  contour  plot 

=  5  Differential  MMD  contour  plot 

=  6  Maximum  manoeuvre  persistence  (MMP)  calculation 


=; 

7  Plot  of  turn  rate  and 

energy  parameter  versus  Mach  number 

COMMON  Output 

COMMON  B: 

SWP  PLA  W 

DTEMP 

MSECS0  DELES 

NES 

DELXM  NXM 

FRAME  DAY 

COMMON  TABLES:  GNTAB 

DF.LGN 

CLOCK 

NGN  ES0  XM0 

ROLE 

GN0 

COMMON  E: 

ESP  IN  XMF1N 

NAMOUT 

GNF1N 

RATING  AMET 

User  Routine  Called 

UNITS  1NMMD 


System  Routines  Called 
TIMES 

Calling  Routine 
P2  program 


A7.10  RATE!  Subroutine 

Purpose 

RATEl  reads  data  produced  by  program  AIRCRAFT  from  disk  expressed  as  energy 
parameter  versus  turn  rate,  for  a  range  of  altitudes  (pressure  height  or  energy  height)  and  Mach 
numbers,  and  prepares  a  plotter  file  of  this  data.  A  separate  plotter  page  is  produced  for  each 
value  of  the  altitude  variable. 

COMMON  B  and  COMMON  C  provide  data  communication 
Use 

CALL  RATEl(WORK) 

where  WORK  is  an  8000-word  vector  of  working  storage 
Input 

COMMON  B:  NES  DELXM  NXM  DELGN  NGN  XM0 

COMMON  C:  ENERGY  PSUN1T  FACT  ALUN1T 

User  Routines  Called 

PLTLAB  AXIS  LINE  SYMBOL 

System  Routines  Called 
TEXT  PLOT 

Calling  Routine 
P2  program 


A7.ll  RATE2  Subroutine 

Purpose 

RATE2  reads  data  produced  by  program  AIRCRAFT  from  disk  representing  the  maximum 
Manoeuvre  Diagram  (,MMD).  This  provides  optimum  energy  parameter  and  turn  rate  versus 
load  factor  for  a  range  of  energy  states.  A  linear  interpolation  is  performed  and  a  disk  file  P2.0PT 
for  input  to  program  P4  is  prepared.  This  data  provides  turn  rate  in  an  (energy  state,  energy  para¬ 
meter)  grid. 

Data  communication  is  via  COMMON  areas  B,  E,  TABLES  and  ESL1ST. 

Use 

CALL  RATE2(WORK) 

where  WORK  is  an  8000-word  vector  of  working  storage. 

Input 

COMMON  B:  DELES  NES  DELGN  NGN  ES0  GN0 

COMMON  E:  NAMOUT 

User  Routines  Called 
WRLAB  1NTRP 

System  Routine  Called 
ABS 

Calling  Routine 

R2DIFF  P2  program 


A7.12  R2DIFF  Subroutine 


Purpose 

R2DIFF  performs  the  functions  of  RATE2  with  two  separate  disk  files  representing  com¬ 
parison  aircraft  and  prepares  output  files  P2.0PT  and  P2A.OPT  for  input  to  program  P4. 
While  doing  this  the  grid  specifications  are  checked  for  equality,  and  a  third  differential  file, 
P2D1FF.OPT  is  produced  for  input  tc  program  P4.  An  error  halt  occurs  if  the  data  grids  differ 
in  any  way. 

Data  communication  is  via  COMMON  areas  B,  E  and  ESL1ST. 

Use 

CALL  R2DlFF(WORK) 
where  WORK  is  an  8000-word  vector  of  working  storage. 

Input 

COMMON  B:  MSECS0  DELES  NES  DELGN  NGN  ES0 

GN0  IPSTYP  RJNITS 

Output 

COMMON  B:  SWP  PLA  W  DTEMP  FRAME  DAY  CLOCK 

ROLE  IPSTYP  IUN1TS 
COMMON  £:  RATING  AMET  NAMOUT 

COMMON  ESLIST.  OMTAB 

User  Routines  Called 
RATE2  INMMD 

System  Routines  Called 

CSTR1NG  TIMES  RUNPRG 

Calling  Routine 
P2  program 


A7.I3  UNITS  Subroutine 
Purpose 

UNITS  allocates  text  strings  to  variables  for  use  in  terminal  dialogue  and  plotter  axis 
labelling,  dependent  on  the  values  of  variables  IUNITS  and  IPSTYP  read  from  input  data  files. 
Data  communication  is  via  COMMON  areas  B  and  C. 

Use 

CALL  UNITS 
Input 

COMMON  B:  IPSTYP  IUNITS 

Output 

COMMON  C.  TYPE  ENERGY  PSUNIT  FACT  ALUNIT 

AM  UNIT 


Calling  Routine 
P1IN 


APPENDIX  S 


P24LIB  Subroutine  Library 

A8.1  Nature  of  Library 

This  relocatable  binary  library,  used  when  loading  programs  P2  and  P4,  is  derived  from 
several  sources.  It  consists  of  routines  for  numerical  interpolation,  common  input/output 
routines,  plotting  routines  and  a  routine  to  allocate  texts  for  the  two  unit  systems.  The  routines 
are  described  briefly  below;  fuller  details  may  be  found  in  the  sources  indicated. 

A8.2  Routines  from  Library  PI 

Routine  INTRP  (numerical  interpolation)  and  routine  BAD1NP  (input  error  text)  are 
obtained  from  subroutine  library  PI  which  is  fully  described  in  Appendix  3, 


A8.3  Routines  from  Library  P2 

Input/output  routines  INLAB  andPLTLAB,  and  text  allocation  routine  UNITS  are  obtained 
from  program  library  P2  which  is  fully  described  in  Appendix  6. 


A8.4  Plotting  Subroutines 

Routines  AXIS,  LINE,  NUMBER  and  SYMBOL  are  obtained  from  the  CALCOMP 
plotter  software  library  in  use  at  ARL’s  Computer  Centre.  Full  machine-readable  documentation 
may  be  obtained  by  running  program  ALLPL3.FOR  available  on  DECtape  17  at  the  Computer 
Centre.  A  brief  statement  of  the  purpose  of  each  routine  is  given  below. 


A8.1.I  AXIS  Subroutine 

This  routine  draws  labelled  axes  suitable  fer  plotting  graphs.  The  axes  may  be  drawn  et  any 
specified  angle  from  an  arbitrary  origin.  It  calls  P24LIB  routines  SYMBOL  and  NUMBER  and 
system  routine  PLOT. 


A8.4.2  LINE  Subroutine 

This  routine  plots  a  line  on  the  plotter  using  a  set  of  co-ordinates  stored  in  input  vectors.  It 
calls  P24LIB  routine  SYMBOL  and  system  routines  PLOT  and  WHERE. 


A8.4.3  NUMBER  Subroutine 

This  routine  plots  the  value  of  a  real  number  on  the  plotter,  calling  P24LIB  routine  SYMBOL 
and  system  routine  PLOT 

A8.4.4  SYMBOL  Subroutine 

This  routine  plots  ASCII  characters  and  special  symbols  on  the  plotter,  calling  system 
routines  PLOT,  SIN  and  COS. 
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APPENDIX  9 


A9.1  P4  Program 


P4  Program  library 


Purpose 

The  program  P4  processes  data  files  produced  by  program  P2  and  produces  contour  plots 
for  off-line  plotting. 

The  contour  plots  on  file  P4.PLT  are  of  four  types  (input  filenames  shown  in  parentheses): 

(a)  energy  rate  contour  plots  (P2.CON,  P2A.CON); 

( b )  differential  energy  rate  contour  plots  (P2DIFF.CON); 

(c)  maximum  manoeuvre  diagrams  (P2.0PT,  P2A.OPT); 

(</)  differential  maximum  manoeuvre  diagrams  (P2DIFF.OPT). 


The  type  of  plot  is  determined  by  the  input  filename. 

A  conversational  dialogue  allows  the  user  to  vary  several  features  of  the  plots: 


(i)  smoothness  of  the  plotted  output; 

(ii)  texture  of  critical  contours; 

(iii)  number  and  level  of  contours; 

(iv)  selection  of  load  factors 

(v)  inclusion  of  energy  state  contours 


jenetgy  rate  plo's  only. 


Energy  rate  contours  arc  plotted  on  a  Mach  number  versus  altitude  grid  and  maximum 
manoeuvre  diagrams  present  turn  rate  contours  on  an  energy  state  versus  energy  parameter  grid. 

The  brief  main  program  performs  preliminary  input  operations  and  calls  routine  P4MA1N 
to  perform  the  major  control  operations. 

Plotted  output  is  written  on  logical  unit  I.  COMMON  areas  B  and  C  are  used  for  variable 
communication. 


Use 

RUN  P4  (monitor  command) 

User  Routines  Called 

P4MAIN  UNITS  INLAB 


A9.2  OU TXT  Subroutine 

Purpose 

OUTXT  is  called  by  routine  P  at  the  start  of  a  contour  line  to  output  the  contour  level  on 
the  plotter. 

Use 

CALL  OUTXT(N) 
with  input 

N  The  contour  level  to  be  plotted 

System  Routines  Called 
PLOT  TEXT 

Calling  Routine 
CONT 


49,3  PSateMrtfce 

Purpose 

Subroutine  P  controls  pen-up  and  pen-down  increment?,  taking  into  account  zero  turn 
rate  contour  boundaries  when  plotting  differential  maximum  manoeuvre  diagrams. 

Plotting  is  achieved  by  calls  to  PLOTD  or  SMOOTH,  for  straight  line  or  smoothed  plots, 
respectively.  The  pen  is  lifted  (N  reset  to  5)  if  the  requested  point  is  outside  either  turn  rate 
boundary  (differential  MMD  only). 

Data  communication  is  via  COMMON  areas  B,  UMS,  CONPLT,  NLEVfl,  and 
ODMSWT. 

Use 

CALL  PKX,  Y,  N) 

Input 

X.Y  Co-ordinates  of  the  contour  point  in  units  of  X  and  Y  grid  intervals. 

N  Integer  supplying  pen  position  requested  by  the  contouring  routine: 

3  lift  pen,  move  to  (X,Y),  drop  pen; 

2  drop  pen,  move  to  (X,Y); 

5  lift  j*?a,  more  to  (X,Y). 

COMMON  Input 

COMMON  B:  IC  YLMIN  1DATA  NDATA  ICOUNT 
COMMON  L1MS:  XTAB  YTAB 

COMMON  CONPLT:  CNOLAB  XSCALE  YSCALE  SMTH  1DOT 
COMMON  NLEV0:  NLEV0  IDOT0  IDEP 
COMMON  ODMSWT:  NLEV 

COMMON  Output 

COMMON  B  :  IC  ICOUNT  NSAV 

COMMON  CONPLT:  I  DOT 

User  Routines  Called 

PLOTD  SMOOTH  1NTRP  OUTX1 

System  Routine  Called 
AMIN1 

Calling  Routine 

CONT  DIAG 


A9.4  PLOTD  Subroutine 

Purpose 

PLOTD  is  a  dummy  interface  between  the  plot  calls  in  routines  P  and  SMOOTH  and  the 
CALCOMP  software  routine  PLOT.  Use  of  this  interface  renders  the  contouring  routines 
independent  of  conventions  for  plotter  pen  commands.  In  this  application  it  is  needed  tc  change 
the  pen-down  command  from  N  *  2  (COSNT)  to  N  =  4  (PLOT). 

Use 

C4LL  PLOTD  (X,  Y.  N) 
with  inputs  X,  Y  and  N  defined  as  for  routine  P. 

System  Routine  Calhd 
PLOT 

Calling  Routines 
P  SMOOTH 
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A9JS  P4MAIN  Ssferowta* 

Purpose 

P4MAIN  is  the  control  routine  for  the  contour  plotting  program  P4.  It  uses  the  input  data 
filename  to  differentiate  between  requests  for  single  aircraft  or  comparison  plots,  and  the  file¬ 
name  extension  to  differentiate  between  requests  for  maximum  manoeuvre  diagrams  (turn  rate 
contour  plots)  and  energy  rate  contour  plots.  In  the  latter  case,  a  preliminary  plot  of  energy 
state  contours  may  also  be  requested.  F4MA1N  controls  the  setting  uo  and  plotting  of  identifying 
text  and  rxes,  the  reading  of  grid  data,  and  the  calling  of  routine  CONT  to  plot  the  contours. 

COMMON  areas  B,  C,  CONPLT,  NI  EVB,  ODMDOT,  and  CDMSWT  are  used  for  data 
communication. 

Use 

CALL  P4MAIN  (Z,  ZG,  NXM,  HES) 

Input 

Z  The  (empty)  data  buffers  for  turn  rate  or  energy  rate  data. 

ZG  The  (empty)  data  buffer  for  energy  state  data 

NXM  The  number  of  points  of  the  r-axes  grid  (energy  state  or  Mach  numbers). 

NES  The  number  of  points  on  the  y-axes  grid  (energy  rate  or  pressure  altitude). 

COMMON  Input 

COMMON  B:  NAM  DELES  DELXM  NGN  ES0  XM0 
COMMON  C:  ENERGY  PSUNIT  FACT  ALUNIT 

User  Routines  Called 

SMOOTH  BADINP  SYMBOL  AXIS  CONT  INLAB  PLTLAB 
System  Routines  Called 

TEXT  AMAX1  FLOAT  ARS  WHERE  PLOT 

Calling  Routine 
P4  program 
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APPENDIX  10 

Routines  from  Libraries  GRAFIC  and  EXTRAS 
A10.1  Nature  of  Libraries 

GRAFIC1*  is  a  system  for  graphical  presentation  of  three-dimensional  fluid  flows,  and 
EXTRAS  is  a  library  used  by  GRAFIC.  Both  of  these  libraries  are  available  on  disk  area 
[1033, 1022]  for  use  by  ARL  staff. 

The  contouring  routine  CONT  and  associated  routines  DIAG,  SMOOTH,  REAL1N  and 
PROMPT  form  but  a  small  subset  of  the  GRAFIC  system  required  by  program  P4.  These 
routines  are  outlined  brefly  below. 


A10.2  CONT  Subroutine 


Purpose 

CONT  is  a  routine  in  GRAFIC  for  drawing  a  contour  map  based  on  data  supplied  for  a 
regular  rectangular  grid. 

A  feature  of  CONT  is  that  it  is  independent  of  plotter  conventions,  calling  a  user-written 
interface  routine  P,  which  then  processes  the  plot  request. 

COMMON  area  GDMSWT  is  used  for  additional  data  communication. 

Use 

CALL  CONT  (T,  M,  N,  CONLEV,  NC,  LABFLG) 


Input 


1 

M 

N 

CONLEV 

NC 

LABFLG 


The  array  containing  the  grid  values. 

The  first  dimension  of  T. 

The  second  dimension  of  T. 

The  vector  containing  the  requested  contour  levels. 

The  number  of  contour  levels. 

A  logical  switch  set  to 

.TRUE,  if  the  contour  level  is  *o  be  output  at  the  start  of  each  contour, 
.FALSE,  otherwise. 


Output 

Data  is  plotted  via  calls  to  P(X,  Y,  NN)  with 
X  The  x-co-ordinate  in  the  range  (I,  FLOAT(M)). 

Y  The  y-co-ordinate  in  the  range  (1,  FLOAT  (N)). 

NN  The  pen  code : 

NN  —  3  plot  with  pen  up, 

NN  =  2  plot  with  pen  down. 

COMMON  Input 

COMMON  GDMSWT:  DGS,  a  logical  variable  set  to  .TRUE,  in  P4MA1N  to  bypass 

calls  to  DIAG. 


COMMON  Output 

COMMON  GDMSWT :  NLEV,  the  level  number  of  the  contours  currently  being 

plotted. 


User  Routines  Called 
DIAO  P 


Calling  Routine 
P4MAIN 

A 103  D1AG  S«brovdM 
Purpose 

DIAO  is  a  routine  in  GRAF1C  called  by  CONT  to  interpolate  along  a  mesh  diagonal 
using  a  simple  biquadratic  representation  of  the  mesh.  As  used  in  program  P4,  diagonal  interpo¬ 
lation  is  bypassed  by  setting  switch  DGS  to  'TRUE'. 

COMMON  area  GDMSWT  is  used  for  additional  data  communication. 

Use 

CALL  DIAG  (T,  M,  Ml,  Nl,  I,  J,  Cl) 

Input 

T  The  array  of  grid  values. 

M  The  first  dimension  of  T. 

Ml,  Nl  Direction  along  M  and  N  direction,  respectively: 

+ 1  denotes  positive, 

—  1  denotes  negative. 

(I,  J)  Index  of  current  grid  reference. 

C  Current  contour  level. 

COMMON  Input 

COMMON/GDMSWT/DGS 

User  Routine  Called 
P 

System  Routine  Called 
SQRT 

Calling  Routine 
CONT 

A10.4  SMOOTH  Subroutine 
Purpose 

SMOOTH  is  a  routine  in  GRAF1C  which  smooths  the  contour  plots  by  joining  successive 
points  on  the  contour  with  cubic  arcs  having  tangential  coincidence  at  their  common  points 
(knots).  The  tangential  slope  is  set  equal  to  the  mean  of  the  slopes  of  the  straight  lines  that  would 
otherwise  join  points  on  either  side  of  the  knot.  When  smoothing  is  selected,  an  additional  call 
to  SMOOTH  is  required  to  complete  the  current  contour  before  proceeding  to  the  next  contour 
request. 

Common  area  GDMDOT  is  used  for  additional  data  transmission. 

Use 

CALL  SMOOTH  (X,  Y,  N,  IDOT) 

Input 

X,  Y  The  co-ordinates  of  the  point  to  be  plotted. 

N  The  pen-up  command : 

N  -  2  put  pen  down  and  draw  to  (X,  Y), 

N  «  3  move  to  (X,  Y)  with  pen  up  and  drop  pen, 

N  =  4  finish  smoothing  the  current  contour. 
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I  DOT  An  integer  specifying  the  mark:  space  ratio  for  broken  lines: 
IDOT  -  0  continuous  line, 

1DOT  >  1  space  is  larger  than  mark, 

IDOT  <  I  mark  is  greater  than  space, 

IDOT  -  1  mark  »  space. 


COMMON  Input 

COMMON  GDMDOT:  DINT,  the  smoothing  interval,  set  to  0  05  in,  »n  routine 

P4MAIN. 

User  Routines  Called 

PROMPT  REAUN 

System  Routines  Celled 
SQRT 

Calling  Routines 
P4MAIN  P 

AI0.5  REAUN  Subroutine 
Purpose 

REAUN  is  a  routine  in  GRAFIC,  which  accepts  real  numerical  data  from  the  user’s 
terminal.  It  is  not  cailed  during  the  execution  of  P4,  but  is  required  to  complete  the  loading 
process. 

User  Routines  Called 
PROMPT 

Calling  Routine 
SMOOTH 

A10.6  PROMPT  Routine 

Purpose 

PROMPT  is  a  routine  in  library  EXTRAS,  which  writes  a  text  on  the  user's  terminal.  It  is 
not  called  during  the  execution  of  P4,  but  is  required  to  complete  the  loading  process. 

Calling  Routines 

SMOOTH  REAUN 
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APPENDIX  11 


System  Routine  Requirements 


Routine 

Required  by 

Brief  description 

E 

Ki 

2 

ca 

3 

a. 

P24LIB 

u 

E 

< 

K 

O 

ABS 

X 

X 

X 

Absolute  value  of  real  argument 

FLOAT 

X 

X 

Integer  to  floating  point  conversion 

AMIN1 

X 

X 

Mihimum  of  real  arguments 

AMAX1 

X 

Maximum  of  real  arguments 

MIN0 

X 

Integer  minimum  of  integer  arguments 

ALOCi 

X 

Natural  logarithm 

EXP 

X 

Exponential 

SQRT 

X 

X 

X 

Square  root 

SIN 

X 

Sine 

ASIN 

X 

Inverse  sine 

COS 

X 

X 

Cos 

COSD 

X 

Cos,  argument  in  degrees 

DATE 

X 

X 

Current  date 

TIME 

X 

X 

Current  daytime 

TIMES 

X 

Current  C  PU  time  used 

RUNPROG 

V 

•  ■ 

Transfer  control  to  another  program 

CSTRlNG 

X 

Generate  CUSP  command  string 

PLOT 

X 

X 

X 

Plot  a  point  on  the  plotter 

TEXT 

X 

X 

Write  a  text  on  toe  plotter 

WHERE 

X 

X 

1 

Provide  current  plotter  co-ordinates 
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